Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,184,318 T→G F1025C (TTC→TGC)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi4,184,3180TG97.3% 502.6 / 2.9 147F1025C (TTC→TGC) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base G (80/63);  total (82/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.57e-01

Only 100 of 147 total aligned reads displayed.

TGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGG  >  NC_000913_3_bme_pgi/4184172‑4184464
                                                                                                                                                  |                                                                                                                                                  
tgctgtgctgGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCga                                                                                                                                                 >  2:1653140/1‑150 (MQ=255)
 gctgtgctgGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCgag                                                                                                                                                <  1:374325/150‑1 (MQ=255)
   tgtgctgGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCgagat                                                                                                                                              >  1:348436/1‑149 (MQ=255)
     tgctgGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGagaaga                                                                                                                                            >  2:1207604/1‑150 (MQ=255)
       ctgGTAGCCGGTGGCGTTGAAGCTGAGAAGCTAGACAAACTGCCGCGCGATCGCTGGCTGGAGCGGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTNAGCAGTATGACGNACTGAAACACGAGTGCGAGAAGaaa                                                                                                                                          <  3:996772/150‑1 (MQ=255)
        tgGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAAc                                                                                                                                         >  1:1645500/1‑150 (MQ=255)
          gTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCCGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTc                                                                                                                                       >  2:537500/1‑150 (MQ=255)
            aGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGa                                                                                                                                     >  1:1217357/1‑150 (MQ=255)
            aGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGa                                                                                                                                     >  3:701638/1‑150 (MQ=255)
              gcGGTGGCGTTGAAGGTGAGGAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGGTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAATCTCGAAg                                                                                                                                   <  2:155539/149‑1 (MQ=255)
              ccGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAACTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAg                                                                                                                                   <  4:677153/150‑1 (MQ=255)
                      gTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGAGAACAGAAACTCGAAGCGAAAcgc                                                                                                                           >  2:241775/1‑150 (MQ=255)
                         gAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAAcgccgc                                                                                                                        >  1:195941/1‑150 (MQ=255)
                         gAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAAcgccgc                                                                                                                        <  1:994040/150‑1 (MQ=255)
                          aaGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCa                                                                                                                       >  2:559526/1‑150 (MQ=255)
                            gCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGATTTGCGCCGAACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCaaa                                                                                                                     <  1:241775/150‑1 (MQ=255)
                                 gaAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTc                                                                                                                                       >  1:7210/1‑127 (MQ=255)
                                      tCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAAcgcc                                                                                                                          >  1:1450275/1‑135 (MQ=255)
                                       cGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAggg                                                                                                          <  1:1653140/150‑1 (MQ=255)
                                        gACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGc                                                                                                         >  3:845809/1‑150 (MQ=255)
                                         aCAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGcg                                                                                                        >  1:1513608/1‑150 (MQ=255)
                                           aaaCTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAAt                                                                                                                   <  1:1469345/137‑1 (MQ=255)
                                           aaaCTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAAt                                                                                                                   >  2:1469345/1‑137 (MQ=255)
                                              cTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAGGCGaaa                                                                                                                              >  3:602270/1‑123 (MQ=255)
                                                gCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCt                                                                                                 <  1:1126125/150‑1 (MQ=255)
                                                  cgcgcgATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTgg                                                                                               <  3:1002510/150‑1 (MQ=255)
                                                  cgcgcgATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTgg                                                                                               >  4:140457/1‑150 (MQ=255)
                                                   gcgcgATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGc                                                                                              <  3:20485/150‑1 (MQ=255)
                                                     gcgATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCAc                                                                                            >  4:587668/1‑150 (MQ=255)
                                                         tCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCggg                                                                                        >  3:14577/1‑150 (MQ=255)
                                                            ctggctggAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGt                                                                                     <  2:806056/150‑1 (MQ=255)
                                                              ggctggAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGc                                                                                   <  2:1281590/150‑1 (MQ=255)
                                                                 tggAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGa                                                                                >  2:292881/1‑150 (MQ=255)
                                                                  ggAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGaa                                                                               <  2:1180058/150‑1 (MQ=255)
                                                                     gCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGAt                                                                            >  1:1145462/1‑150 (MQ=255)
                                                                     gCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGAt                                                                            <  4:520432/150‑1 (MQ=255)
                                                                      cTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTg                                                                                    <  2:529679/141‑1 (MQ=255)
                                                                      cTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTg                                                                                    >  1:529679/1‑141 (MQ=255)
                                                                        gggCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGt                                                                         >  3:38448/1‑150 (MQ=255)
                                                                          gCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTa                                                                       <  2:1645500/150‑1 (MQ=255)
                                                                             tGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGt                                                                                     >  2:1253487/1‑133 (MQ=255)
                                                                             tGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGt                                                                                     <  1:1253487/133‑1 (MQ=255)
                                                                             tGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGAtt                                                                           >  1:121537/1‑143 (MQ=255)
                                                                             tGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGAtt                                                                           <  2:121537/143‑1 (MQ=255)
                                                                              gACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGCGCAGTATGACGAACTGAAACACGAGTGCGAGACGAAACTCGAAGCGAACCGACGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGATGAAGATTGTTAAGGt                                                                   >  3:628753/1‑150 (MQ=255)
                                                                              gACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGt                                                                   <  2:1322921/150‑1 (MQ=255)
                                                                                 aGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGtata                                                                >  2:1069861/1‑150 (MQ=255)
                                                                                  gACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTatat                                                               >  2:760618/1‑150 (MQ=255)
                                                                                   aCGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGcga                                                                                                       <  1:220382/109‑1 (MQ=255)
                                                                                   aCGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGcga                                                                                                       >  2:220382/1‑109 (MQ=255)
                                                                                      aaGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAAc                                                                                                                                         >  1:7437/1‑72 (MQ=255)
                                                                                      aaGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAAc                                                                                                                                         <  2:7437/72‑1 (MQ=255)
                                                                                      aaGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTgg                                                           <  1:1311129/150‑1 (MQ=255)
                                                                                         agaAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCgg                                                        <  1:138156/150‑1 (MQ=255)
                                                                                           aaaCAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGaaa                                                                                                                                          >  1:101221/1‑66 (MQ=255)
                                                                                           aaaCAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGtt                                                      >  2:686307/1‑150 (MQ=255)
                                                                                            aaCAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTa                                                     >  2:1419798/1‑150 (MQ=255)
                                                                                            aaCAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCAcc                                                                                                               <  2:577540/92‑1 (MQ=255)
                                                                                            aaCAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCAcc                                                                                                               >  1:577540/1‑92 (MQ=255)
                                                                                              cAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCGGAGCGGCGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTaaa                                                   >  1:1192538/1‑150 (MQ=255)
                                                                                              cAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCgg                                                                                         >  1:619150/1‑112 (MQ=255)
                                                                                              cAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCgg                                                                                         <  2:619150/112‑1 (MQ=255)
                                                                                              cAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTaaa                                                   <  1:1446804/150‑1 (MQ=255)
                                                                                              cAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTaaa                                                   >  4:392802/1‑150 (MQ=255)
                                                                                              cAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTaaa                                                   <  4:244292/150‑1 (MQ=255)
                                                                                                  aTCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGcc                                               <  3:257661/150‑1 (MQ=255)
                                                                                                       cTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATc                                          <  1:805519/150‑1 (MQ=255)
                                                                                                        tGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTaa                                                                      >  1:650568/1‑121 (MQ=255)
                                                                                                        tGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTaa                                                                      <  2:650568/121‑1 (MQ=255)
                                                                                                         ggAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGAtt                                                                           <  1:1280213/115‑1 (MQ=255)
                                                                                                              aGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTg                                   >  3:598493/1‑150 (MQ=255)
                                                                                                              aGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTg                                   <  2:232977/150‑1 (MQ=255)
                                                                                                              aGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTg                                   >  1:765605/1‑150 (MQ=255)
                                                                                                                ctggctgAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAg                                                                     <  2:760949/114‑1 (MQ=255)
                                                                                                                ctggctgAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGt                                 <  1:559526/150‑1 (MQ=255)
                                                                                                                  ggctgAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGAAGCAAAATCACACAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGa                               <  3:1032974/150‑1 (MQ=255)
                                                                                                                   gctgAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCgg                                                                                         >  1:539519/1‑91 (MQ=255)
                                                                                                                   gctgAGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGAc                              >  2:1258279/1‑150 (MQ=255)
                                                                                                                       aGCAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGa                          >  4:103637/1‑150 (MQ=255)
                                                                                                                         cAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATg                        >  3:1031336/1‑150 (MQ=255)
                                                                                                                         cAGTATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATg                        >  2:711955/1‑150 (MQ=255)
                                                                                                                            tATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCa                     <  3:392802/150‑1 (MQ=255)
                                                                                                                            tATGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCa                     >  3:564064/1‑150 (MQ=255)
                                                                                                                              tGACGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAgg                   >  2:372998/1‑150 (MQ=255)
                                                                                                                                aCGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGCGAAGATGGCAGgtc                 >  2:437038/1‑150 (MQ=255)
                                                                                                                                aCGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGgtc                 <  1:292881/150‑1 (MQ=255)
                                                                                                                                aCGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCCGGTGACAAGATGGCAGgtc                 <  4:14577/150‑1 (MQ=255)
                                                                                                                                 cGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTg                                                                                 <  4:751939/85‑1 (MQ=255)
                                                                                                                                 cGAACTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTg                                                                                 >  3:751939/1‑85 (MQ=255)
                                                                                                                                   aaCTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGgtcgtc              >  2:331207/1‑150 (MQ=255)
                                                                                                                                   aaCTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTAGTGACAAGATGGCAGgtcgtc              <  3:512102/150‑1 (MQ=255)
                                                                                                                                    aCTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGta                                                                  <  4:27759/97‑1 (MQ=255)
                                                                                                                                    aCTGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGta                                                                  >  3:27759/1‑97 (MQ=255)
                                                                                                                                      tGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACg           >  2:370898/1‑150 (MQ=255)
                                                                                                                                      tGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACg           <  4:439499/150‑1 (MQ=255)
                                                                                                                                      tGAAACACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACg           >  2:480975/1‑150 (MQ=255)
                                                                                                                                        aaaCACGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGt         >  3:720554/1‑150 (MQ=255)
                                                                                                                                            acGAGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACa     >  2:1356049/1‑150 (MQ=255)
                                                                                                                                              gAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGgtcgt               <  1:269313/138‑1 (MQ=255)
                                                                                                                                               aGTGCGAGAAGAAACTCGAAGCGAAACGCCGCAACATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAgg  <  1:1204073/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGG  >  NC_000913_3_bme_pgi/4184172‑4184464

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑