Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,326,088 T→G 5.6% intergenic (‑347/+311) yjdN ← / ← yjdM metalloprotein superfamily protein/zinc‑ribbon family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,326,0880TG 5.6% 386.3 / 11.8 132intergenic (‑347/+311)yjdN/yjdMmetalloprotein superfamily protein/zinc‑ribbon family protein
Reads supporting (aligned to +/- strand):  ref base T (63/61);  new base G (4/4);  total (67/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 = 1.00e+00

Only 100 of 136 total aligned reads displayed.

GTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGC  >  NC_000913/4325939‑4326237
                                                                                                                                                     |                                                                                                                                                     
gTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAt                                                                                                                                                       <  2:1598597/150‑1 (MQ=38)
gTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAt                                                                                                                                                       <  2:1214962/150‑1 (MQ=38)
      aCCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCg                                                                                                                                                 >  1:1866143/1‑150 (MQ=35)
          tggCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGaa                                                                                                                                             <  1:429429/148‑1 (MQ=17)
             ccGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAAttt                                                                                                                                          <  1:1841087/150‑1 (MQ=32)
             ccGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAAttt                                                                                                                                          <  2:369969/150‑1 (MQ=32)
                gATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTa                                                                                                                                       >  2:2298876/1‑150 (MQ=32)
                gATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTa                                                                                                                                       <  1:2502008/150‑1 (MQ=32)
                 aTGCGGCGTAACCGCCTTATCCGGCCTACGGAGGGTGCGGGACTTTGGAGGCCTGCTAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTACACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAg                                                                                                                                      <  1:2619197/150‑1 (MQ=11)
                 aTGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCTGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAg                                                                                                                                      <  1:2614586/150‑1 (MQ=17)
                        gTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGAt                                                                                                                               <  1:2114623/150‑1 (MQ=32)
                        gTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGAt                                                                                                                               <  1:2114629/150‑1 (MQ=32)
                              gCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGAcg                                                                                                                         <  1:958404/150‑1 (MQ=32)
                               ccTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGAcgc                                                                                                                        <  2:2616578/150‑1 (MQ=32)
                                cTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGAcgcg                                                                                                                       >  1:381644/1‑150 (MQ=32)
                                    tCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAg                                                                                                                   >  1:49685/1‑150 (MQ=32)
                                       ggCCTTCGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGt                                                                                                                <  2:276533/150‑1 (MQ=17)
                                        gCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTc                                                                                                               >  1:1143679/1‑150 (MQ=32)
                                        gCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTc                                                                                                               >  1:1143399/1‑150 (MQ=32)
                                        gCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTc                                                                                                               <  1:668103/150‑1 (MQ=32)
                                             cGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATc                                                                                                          <  1:1912392/150‑1 (MQ=32)
                                             cGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATc                                                                                                          <  1:2040015/150‑1 (MQ=32)
                                                   gTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCat                                                                                                    >  1:1077335/1‑149 (MQ=14)
                                                   gTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCat                                                                                                    >  1:1077353/1‑149 (MQ=14)
                                                    tGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGt                                                                                                   >  2:1107470/1‑150 (MQ=33)
                                                     gCGGGAATTTGTAGGCCTGATAAGACGCGCTAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGtt                                                                                                  <  1:1596381/150‑1 (MQ=21)
                                                     gCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGtt                                                                                                  <  1:714536/150‑1 (MQ=34)
                                                      cGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACCGATGGACAGGGAATTTTTAGTCCTGATAAGACGCGCTATAGTCGCATCAggaatatg                                                                                                 >  2:851083/1‑144 (MQ=11)
                                                       gggAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTgg                                                                                                >  2:75081/1‑150 (MQ=35)
                                                       gggAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTga                                                                                                >  2:2090324/1‑149 (MQ=25)
                                                        ggAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGc                                                                                               >  2:1273351/1‑150 (MQ=37)
                                                         gAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCa                                                                                              >  1:1554652/1‑150 (MQ=37)
                                                         gAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGGGAGCGCCTTTTACGGCCTAGCGATGGTGGGGTAAAGTGTATGCCTGATAAATCGCGCAAGCGGCGCATCGGGCAGATGGcc                                                                                              >  1:510964/1‑149 (MQ=11)
                                                             ttGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGt                                                                                          <  1:1470781/150‑1 (MQ=39)
                                                                 aGGCCCGATACGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTCAACGCCTGACCCGGCCACCGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGcc                                                                                      <  1:2230391/150‑1 (MQ=14)
                                                                  ggCCTGATAAGACGCGCAAGCGTCGCCTCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGCTAAGACGCGCAAGCGTCGCATCGGGCAGTTGGCACTGGTGCCg                                                                                     <  1:72402/150‑1 (MQ=255)
                                                                   gCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCg                                                                                                                                                 <  2:1222514/89‑1 (MQ=18)
                                                                    ccTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCCGTTGGCACTGTTGCCGGa                                                                                   >  2:551183/1‑150 (MQ=255)
                                                                             aCGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATgcggc                                                                             >  2:2590869/1‑147 (MQ=39)
                                                                               gcgcAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGGTGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTaaa                                                                        >  1:305278/1‑150 (MQ=255)
                                                                               gcgcAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAAAGCCTTATCCGGCCGACGGATGGTGCGGGAATTTGTTGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACAGTTGCCGGGTGCGGCGTtaa                                                                        >  1:2045228/1‑150 (MQ=11)
                                                                               gcgcAAGCGTCGCATCAGGCAGTCGGCACCGTGTCCGGATGCGGCGTAAACGCCATATCCGGCCTACGGATGGTGCGGGAATTTGTTGGCCAGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTTCCGGATGCGGCGTtaa                                                                        >  1:305230/1‑150 (MQ=14)
                                                                               gcgcAAGCGTCGCATCAGGCAGTCGGCACCGTGGCCGGATGCGGCGTAAACGCCTTCTCCGGCCTACGGATGGTGCGGGAATTTGTCGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTGGGCACTGTTGCCGGATGCGGCGTaaa                                                                        <  2:1229755/150‑1 (MQ=255)
                                                                                cgcAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGATGGCACAGTTGCCGGATGCGGCGTAAAc                                                                       >  2:2602628/1‑150 (MQ=255)
                                                                                cgcAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGCTGCGTCAATATGTAGGGCTGATAAGACGCGCAAGGGTCGCAGCAGGCAGATGGCACAGATGCCGGAGGCGGCGTacac                                                                       >  2:2556256/1‑150 (MQ=11)
                                                                                  cAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGc                                                                     <  1:510209/150‑1 (MQ=39)
                                                                                    aGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACGGTTGCCGGATGCGGCGTAAACGCCt                                                                   >  2:483356/1‑150 (MQ=37)
                                                                                      cGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTa                                                                 <  1:556642/150‑1 (MQ=39)
                                                                                         cGCATCAGGCAGTCGGCACCGTTGCCGGAGGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGGATTTGTAGGCCTGATCAGACGCGCAAGCGTCGCATCAGGCAGGTGGCACGGTTGCCGGATGCGGCGGAAGCGCCTTATcc                                                              >  2:888950/1‑150 (MQ=11)
                                                                                            aTCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                           <  1:558348/150‑1 (MQ=255)
                                                                                            aTCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                           <  1:2264580/150‑1 (MQ=21)
                                                                                            aTCAGGCAGTCGGCACCGGTGCCGGATGCGGCGTAACCGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                           <  1:2264626/150‑1 (MQ=255)
                                                                                             tCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                          >  2:1566159/1‑150 (MQ=39)
                                                                                                 gCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCGTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGGTGGCACTGATGCCGGAGGCGGCGTAAACGCCTTATCCGGCCTACg                                                      >  2:1384933/1‑150 (MQ=255)
                                                                                                   aGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGa                                                    <  2:1693931/150‑1 (MQ=39)
                                                                                                         cACCGTGGCCGGATGCGGCGTAACCGCCTTATCCGGCCTACGGATGGTGCGGGTATTTGGAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGc                                              <  1:1661449/150‑1 (MQ=255)
                                                                                                           ctgTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGAcgc                                                                                                                        >  1:2616574/3‑74 (MQ=14)
                                                                                                           ctgTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGAcgc                                                                                                                        <  2:2616574/72‑1 (MQ=14)
                                                                                                             gTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGa                                          >  1:2717805/1‑150 (MQ=37)
                                                                                                             gTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGa                                          >  1:2717789/1‑150 (MQ=37)
                                                                                                              gtGCCGGATGCGGCGTAACCGCCTTCTCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGaa                                         <  2:1876776/149‑1 (MQ=255)
                                                                                                               tGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAt                                        <  1:2298876/150‑1 (MQ=37)
                                                                                                               tGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACGGTTGCCGGATGCGGCGAAAACGCCTTATCCGGCCTACGGAGGGGGCGGGAAt                                        >  2:1152213/1‑150 (MQ=255)
                                                                                                                 ccGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAttt                                      <  1:75081/150‑1 (MQ=37)
                                                                                                                 ccGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCg                                             <  1:586177/143‑1 (MQ=17)
                                                                                                                 ccGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCg                                             >  2:586177/1‑143 (MQ=17)
                                                                                                                   ggAGGCGGCGTACCCGCCGTCGCCGGCCTTGGGCTGGGGCGGGAATTGGGCGGCCTGATACGCCGCGCACGCGGCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGt                                    <  1:1854896/150‑1 (MQ=11)
                                                                                                                    gATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTa                                   <  1:2523435/150‑1 (MQ=37)
                                                                                                                    gATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTa                                   <  1:1130049/150‑1 (MQ=37)
                                                                                                                       gcggcgTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTTGGc                                >  2:909795/1‑150 (MQ=37)
                                                                                                                         ggcgTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCt                              <  2:1143679/150‑1 (MQ=37)
                                                                                                                         ggcgTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                           <  1:558238/121‑1 (MQ=17)
                                                                                                                          gcgTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATgcggc                                                                             <  2:683378/102‑1 (MQ=40)
                                                                                                                           cgTAAACGCCTTATCCGGCCTACGGATGGTGCTGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGa                            >  2:474853/1‑150 (MQ=37)
                                                                                                                           cgTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCa                                                                                                         <  2:1037766/73‑1 (MQ=34)
                                                                                                                           cgTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCa                                                                                                         >  1:1037766/1‑73 (MQ=34)
                                                                                                                           cgTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGa                            >  1:327743/1‑150 (MQ=37)
                                                                                                                           cgTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATATGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGATGGCACTGTTGCCGGATGCGGCGTAAAAGCCTTATCCGGCCGACGGAGGGGGCTGGAATTTGTAGGCCTGa                            >  1:2488274/1‑150 (MQ=11)
                                                                                                                             tAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCTGTTGGCACTGGTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATa                          <  2:1890077/150‑1 (MQ=255)
                                                                                                                              aaaCGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATaa                         <  2:381644/150‑1 (MQ=37)
                                                                                                                               aaCGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAg                                                                                                        <  1:1219962/70‑1 (MQ=35)
                                                                                                                               aaCGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTg                                                                                                 <  1:2336903/77‑1 (MQ=38)
                                                                                                                                  gCCTTATCCGGCCTACGGATGGTGCGGGAGTTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGAcg                     >  1:1931566/1‑150 (MQ=37)
                                                                                                                                  gCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTg                                               <  1:223012/124‑1 (MQ=38)
                                                                                                                                  gCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGAcg                     >  1:2051521/1‑150 (MQ=37)
                                                                                                                                  gCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACAGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGAcg                     >  2:243747/1‑150 (MQ=21)
                                                                                                                                  gCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTg                                               >  2:223012/1‑124 (MQ=17)
                                                                                                                                   ccTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGAcgc                    >  2:1734121/1‑150 (MQ=37)
                                                                                                                                   ccTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAGTTGTAGGCCTGATCAGAcgc                    >  2:1838599/1‑150 (MQ=255)
                                                                                                                                          cGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCg             >  2:954012/1‑150 (MQ=37)
                                                                                                                                          cGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGAGGGAATTTGTAGGCCTGATAAGACGCGCAAGCg             >  2:2063041/1‑150 (MQ=37)
                                                                                                                                            gCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGAcgc                    >  1:1289216/1‑141 (MQ=37)
                                                                                                                                            gCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGAcgc                    <  2:1289216/141‑1 (MQ=37)
                                                                                                                                                aCGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCAt       <  2:2341086/150‑1 (MQ=37)
                                                                                                                                                 cGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATc      >  2:681635/1‑150 (MQ=37)
                                                                                                                                                  ggATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCa     <  1:434654/150‑1 (MQ=37)
                                                                                                                                                  ggAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTTCGGAGGGTGCGGGAATTTGTa                                   >  1:1758236/1‑120 (MQ=12)
                                                                                                                                                  ggAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTa                                   <  2:1758236/120‑1 (MQ=17)
                                                                                                                                                    aTGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAgg   <  1:1566159/150‑1 (MQ=35)
                                                                                                                                                     tGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGc  >  1:1919709/1‑150 (MQ=35)
                                                                                                                                                     |                                                                                                                                                     
GTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGC  >  NC_000913/4325939‑4326237

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: