Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 915,226 T→G 15.6% intergenic (‑369/+126) ybjE ← / ← aqpZ putative transporter/aquaporin Z

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913915,2260TG15.6% 347.4 / 14.8 195intergenic (‑369/+126)ybjE/aqpZputative transporter/aquaporin Z
Reads supporting (aligned to +/- strand):  ref base T (101/42);  new base G (13/18);  total (128/67)
Fisher's exact test for biased strand distribution p-value = 3.38e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 196 total aligned reads displayed.

ATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGT  >  NC_000913/915078‑915375
                                                                                                                                                     |                                                                                                                                                     
   gATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTgcgc                                                                                                                                                   >  1:5325296/1‑150 (MQ=32)
     tgtgggctgtgggcGAGTTTTACTGGGTTGTAACCTCGCGCCCTACGGGGTAGCCTTACGGGTTAAGTCGAGGGCACCGCGCAGCGGCTATTTTTTTGGCCTGTGCCGGGGGGTGCTGGGGGGGGCGGCGATGGGCCGCCCCGTgcgcgc                                                                                                                                                 <  2:3634154/150‑1 (MQ=255)
     tgtgGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTgcgcgc                                                                                                                                                 >  2:4594862/1‑150 (MQ=32)
     ggtgGGGTGTGGGGGGTTTTTACTGGTATTTTACCTCGCGCCCGTCCGGTTAGCCTTATGGGTTATGTCGTGGGCACCGCGCTGGGGCTTTTTTGTTCGCCAGAGCCCGGGGGTGCTGGGGGGGGCGGGGTTTGGCCGCCCCCGgcgcgc                                                                                                                                                 <  2:3633403/149‑1 (MQ=2)
         ggCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCt                                                                                                                                             >  2:5585235/1‑150 (MQ=32)
           cTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTg                                                                                                                                           >  1:1478767/1‑150 (MQ=32)
           cTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTg                                                                                                                                           >  1:6656521/1‑150 (MQ=32)
                ggCGCTTTTTACCGGGATGTTCCCGCGCGCCCGACCGGGTCGCCTTACGGGATATGTGGCGGGCACCGCGCAGCGGCTTTTTTGTTCGCCAGCGCCGGGGGGTGCTGGGGGGGGCGGCGTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCa                                                                                                                                      <  2:5535570/150‑1 (MQ=11)
                ggCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCa                                                                                                                                      >  2:3062891/1‑150 (MQ=32)
                ggCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCa                                                                                                                                      >  1:208613/1‑150 (MQ=32)
                 gcgcGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAg                                                                                                                                     >  2:557280/1‑150 (MQ=32)
                 gcgcGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCGTGCGCCAg                                                                                                                                     >  1:3036550/1‑150 (MQ=17)
                  cgcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAgg                                                                                                                                    >  1:4528875/1‑149 (MQ=17)
                     gTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGc                                                                                                                                 >  1:4888087/1‑150 (MQ=33)
                      ttttACTGGTTTTTTCCCCCGCGCCCGCCCGGGTTGCGTTTCGGGATAAGGCGCGGGCACCGGGCCGGGGGGTTTTTGTTCACCCGCGCCGGGGGGTGCTGGGGGGGGGGGCGTGGGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCa                                                                                                                                <  2:6615461/150‑1 (MQ=2)
                      ttttACTGGTATGTTCCCTGGCGCCCGCCCGGGTAGCCTTATGGGTTAGGTCGTGGGCACCGCGCCGCGGCGTTTTTTTCCGCCAGCGCCGGGGGGTGCTGGGGGGGGGGGCGTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCa                                                                                                                                <  2:3757494/150‑1 (MQ=2)
                      ttttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCa                                                                                                                                >  1:3368214/1‑150 (MQ=34)
                      ttttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCa                                                                                                                                >  1:3368269/1‑150 (MQ=34)
                       tttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCGGCGCGCTCCGTGCGCCAGTGGCaa                                                                                                                               >  1:1959769/1‑150 (MQ=25)
                            aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCGGCGCGCTCCATGCGCCAGTGGCAATATGt                                                                                                                          >  1:2796747/1‑150 (MQ=255)
                             gggATGTTACCTCGCGCCCGACCCGGTAGCCGTATGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGtt                                                                                                                         <  2:4975803/150‑1 (MQ=255)
                              ggATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTg                                                                                                                        >  2:4390228/1‑150 (MQ=40)
                                  gTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTa                                                                                                                    >  2:2615146/1‑150 (MQ=255)
                                   ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAg                                                                                                                   >  2:2536617/1‑150 (MQ=255)
                                     aCCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCt                                                                                                                 >  1:3595478/1‑150 (MQ=255)
                                      ccTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGGTGGGCCGCCCCCAGCGCGCTCCGTGCGCCAGTGGCACTATTTTTGTTAtctc                                                                                                                >  1:1706353/1‑150 (MQ=11)
                                      ccTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTc                                                                                                                >  2:1896718/1‑150 (MQ=255)
                                       cTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCa                                                                                                               >  2:5278311/1‑150 (MQ=255)
                                       cTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCa                                                                                                               >  2:1940879/1‑150 (MQ=255)
                                        tCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCGTGGCCCCCCCCAGCGCGCGCCATGCGCCCGGGGCAAGAAGTTGCGTAGcgccg                                                                                                              >  1:5603254/1‑146 (MQ=11)
                                        tCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCAt                                                                                                              >  2:6350489/1‑150 (MQ=255)
                                         cgcgCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCACCTTGCGCCCGTGGCAATATGTTGCTTAGCTCATg                                                                                                             >  1:6348202/1‑150 (MQ=255)
                                         cgcgCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCACCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATg                                                                                                             >  1:1382848/1‑150 (MQ=255)
                                              ccGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCCGCGCTTGGCCGCCCCCAGCGCGCTCCGTGCGCCACTGGCAAAATGTTGCTTAGCTCATGaaatg                                                                                                        >  1:1177191/1‑148 (MQ=17)
                                              ccGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTATGTTCGCCAGAGCCCGGGCGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCGGCGCGCACCTTGCGCCAGTGGCAAAATGTTGCGTAGCTCATGACAgg                                                                                                        >  2:682495/1‑150 (MQ=14)
                                               cGACCGGGTAGCGTTATGGGATATGTGGCGGGCACCGCGCTGGGGCTATTTTGTGCGCCAGCGCCGGGGGGGGCTGGGGGCGGCGGCTTTGGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGa                                                                                                       <  2:1584649/150‑1 (MQ=17)
                                               cGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAggg                                                                                                       >  1:4028745/1‑149 (MQ=255)
                                                  ccggttAGCCGTATGGTATATGTCGAGGGCACCGCGCTGCGGCTTTTTTGTTCGCCTGCGCCGGGGGGTGCCGGGGGCGGCGGCGTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAgcg                                                                                                    <  2:6070022/145‑1 (MQ=17)
                                                   cgggTAGCCTTATGGTATATGTCGTGGCCACCGCGCAGCGGCTTTTTTGTTCGCCACCGCCGGGGGTTGCTGGGGGCGGCGCCTTTGGCCCCCCCCCCCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAgcgc                                                                                                   <  2:4266455/148‑1 (MQ=12)
                                                   ccGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGACCCCGGGGGTGCACGGCGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAgcgc                                                                                                   >  2:5840533/1‑150 (MQ=255)
                                                          ccGTATGGGATAAGTCGAGGGCCCCGCGCCGGGGCGATTTTGTGCGCCAGAGCCGGGGGGTGCAGGGGGCGGGGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAga                                                                                            <  2:3866307/150‑1 (MQ=18)
                                                          ccGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCCCTCCATGCGCCCGGGGCAAAAAGTTGCTTAGCTCATGAAAGGAGCGCAACAAga                                                                                            >  1:5804057/1‑150 (MQ=17)
                                                           cGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAAAATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgat                                                                                           >  2:5424080/1‑150 (MQ=255)
                                                           cGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGGCCCCCCCCGGCGCGCGCCATGCGCCCGGGGCCAAATGTTGCTTAGCGCAAGAAAGGAGCGCAACAAcat                                                                                           >  1:3696585/1‑150 (MQ=17)
                                                            gTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTTGCCGCCCCCTGCGCCCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgatg                                                                                          >  1:4985133/1‑150 (MQ=255)
                                                             tAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCCGCGATTGGCCCCCCCCATCGCGCACCAATCGCCAATGGCAAAAATTTTCTTATCTCAAGAAACGAGCGCAACAagaaga                                                                                         >  2:6170324/1‑150 (MQ=14)
                                                                 ggATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCGGCGCGCTCCGTGCGCCAGTGGCCATATTTTTCGTAGCTCAGGGAAGGGGCGCAACAAgatgtggta                                                                                     >  1:1063346/1‑145 (MQ=17)
                                                                   aTAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGGTTGGCCGCCCCCAGCGCGCTCCGTGCGCCAGTGGCAATATTTTTGTTATCTCATGAAAGGAGCGCAACAAGATGTTGTAta                                                                                   >  1:104591/1‑149 (MQ=17)
                                                                        tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGGGGCGATTTGGCCCCCCCGGCGCGCTCCATGCGCCAGTGGCAAAATGTTGCTTATCTCATGAAAGGAGCGCAAAAAGAAGATGaaaaaacca                                                                              >  1:3488705/1‑143 (MQ=18)
                                                                        tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCACGGGGCGGCGGCGATTTGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCAAGAAAGGAGCGCAACAAGATGATGAATCAACat                                                                              >  2:2309858/1‑150 (MQ=255)
                                                                        tCGAGGGCACCGCGCAGCGGCGATTGTGTTCGCCCGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGGCCGCCCCCTGCGCGCGCCTTGCGCCAGGGGCCATATGTTGCTTAGCGCATGGAAGGAGCGCAACAAGATGATGAATCacccg                                                                              >  2:79399/1‑146 (MQ=18)
                                                                          gAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCCGGGGTGCACGGGGCGGCGGCGAGTTGCCCCCCCCCGCGCGCTCCGTGCGCCCCTGGCAAAATGTTGCTTAGCGCGAGAAAGGGGCGCCACAAGATTATTGATCAACataa                                                                            >  1:871020/1‑149 (MQ=17)
                                                                           aGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGCTTGGCCGCCCCCGGCGCGCTCCATGCGCCAGGGGCCATATGTTGCTTAGCTCATGAAATGTGCGCAACAAGATTATGAATCAACAaata                                                                           >  1:4159082/1‑150 (MQ=18)
                                                                            gggCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCAGGAAAGGAGCGCAACAAGATGATGAATCAACATATaa                                                                          >  2:5900689/1‑150 (MQ=255)
                                                                               cACGGCGCTGCGGGTTTTTTGTTGGCCTGAGCCGGGGGGTGCTGGGGGGGGCGGGGTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATaacaa                                                                       <  2:4458032/150‑1 (MQ=21)
                                                                                 ccGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAacaaca                                                                     >  2:5989212/1‑150 (MQ=255)
                                                                                          ggcgATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCACAAGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaaa                                                            >  2:219695/1‑150 (MQ=255)
                                                                                           gcgATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGTGCAACAAGATGATGAATCAACATATAACAACATCTTaaaaaa                                                           >  2:3370168/1‑150 (MQ=255)
                                                                                           gcgATTTTGTTCGCCAGAGCCCGGGGGTGCACGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCATGCGCCCGTGGCAATATGTTGCTTAGCTCAGGAACGGAGCGCCACAAGATGATGACTCAACATATAACAACAGCTGaaaaaa                                                           >  2:2334079/1‑150 (MQ=21)
                                                                                             ttttttGTTCGCCCGGGGGGGGGGGGGCGGGGGGGGGGGGGGGTGGGCCCCCCCCGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAg                                                         <  1:237937/148‑1 (MQ=255)
                                                                                             ttttttGTGCGCCTGCGCGGGGGGGGGGGGGGGGGGGGGGGTATGGCCCGCCCCCGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAg                                                         <  1:982682/148‑1 (MQ=255)
                                                                                             taTTTGGTGCGCCTGCGCCGGGGGGGGCTGGGGGCGGGGGCGTTGGCCCCCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAg                                                         <  2:3711597/149‑1 (MQ=255)
                                                                                                tttGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCCCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCAAGAAAGGAGCGCAAAAAGATGATGAATCAACATAAAACAACATCTTaaaaaaacgac                                                      >  2:4707894/1‑146 (MQ=25)
                                                                                                   tttCGCCCGCGCCGGGGGGTGCTGGGGGGGGCGGCTTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGa                                                   <  2:365113/149‑1 (MQ=37)
                                                                                                   gTTCGCCAGCGCGGGGGGGTGCTGGGGGCGGGGGCTTTGGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGa                                                   <  2:1803982/150‑1 (MQ=37)
                                                                                                    ttCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGAc                                                  >  2:4393984/1‑150 (MQ=255)
                                                                                                          agagCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCGCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACg                                            >  2:3139002/1‑150 (MQ=255)
                                                                                                             gcggggggggGCTGGGGGGGGGGGCGTGGGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCa                                         <  2:5871347/146‑1 (MQ=255)
                                                                                                             ccgggggggTGCCGGGGGCGGCGGGTTTGGCCCCCCCCCTCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCa                                         <  1:3254207/146‑1 (MQ=255)
                                                                                                              cccGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCATGCGCCCGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTAAGCCAg                                        >  2:1145482/1‑150 (MQ=255)
                                                                                                               cgggggggGCTGGGGGGGGCGCCGTTGGGCCCCCCCCGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAgg                                       <  1:5968637/148‑1 (MQ=38)
                                                                                                                 gggggTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGcc                                     >  2:5669585/1‑150 (MQ=255)
                                                                                                                         cgggggCGGCGGCGATGGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGt                             <  2:5035422/149‑1 (MQ=255)
                                                                                                                         aGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGTCCTGACATTACGCCAGGCCTTCTGCGt                             >  1:2561840/1‑150 (MQ=255)
                                                                                                                           ggggCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtta                           >  2:6211384/1‑150 (MQ=255)
                                                                                                                           ggggCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtta                           >  2:4093966/1‑150 (MQ=255)
                                                                                                                                 gcggcgATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTtaattaat                     >  2:884092/1‑150 (MQ=255)
                                                                                                                                   ggcgATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCa                   >  1:4414812/1‑150 (MQ=255)
                                                                                                                                    gcgATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCAc                  >  1:2583096/1‑150 (MQ=255)
                                                                                                                                        tgggCCGCCCCGTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGACAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCtt              <  2:326809/148‑1 (MQ=255)
                                                                                                                                        tgggCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCtt              <  1:4437108/148‑1 (MQ=255)
                                                                                                                                         tGGCCGCCCGATGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCttt             <  2:4929043/150‑1 (MQ=255)
                                                                                                                                           gccgccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTc           >  2:3459399/1‑150 (MQ=255)
                                                                                                                                           gccgccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTc           <  1:521095/150‑1 (MQ=255)
                                                                                                                                           gccgccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTc           >  1:3913982/1‑150 (MQ=255)
                                                                                                                                             ccgccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTcc          >  1:4761125/1‑150 (MQ=255)
                                                                                                                                             ccgccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTcc          >  2:3208477/1‑150 (MQ=255)
                                                                                                                                             ccgccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACAGATAACAACATCTTAAAAAATGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTcc          >  2:599710/1‑150 (MQ=255)
                                                                                                                                              cgccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCca         >  1:461264/1‑150 (MQ=255)
                                                                                                                                              cgccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCca         <  1:1594070/150‑1 (MQ=255)
                                                                                                                                              cgccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCca         >  1:4296465/1‑150 (MQ=255)
                                                                                                                                               gccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCcag        >  2:1411879/1‑150 (MQ=255)
                                                                                                                                                cccccTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCcagc       <  2:6085322/150‑1 (MQ=255)
                                                                                                                                                 cgccTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCcagca      <  2:6524532/148‑1 (MQ=255)
                                                                                                                                                 ccccTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCcagca      <  1:3946442/150‑1 (MQ=255)
                                                                                                                                                 ccccTCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCcagca      <  2:5061620/150‑1 (MQ=255)
                                                                                                                                                 ccccCGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCcagca      <  1:6350489/150‑1 (MQ=255)
                                                                                                                                                   ccTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAgg    <  2:4435028/150‑1 (MQ=255)
                                                                                                                                                    cTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGGTAATTAATCACGCTTTTCCAGCAggg   >  1:4794032/1‑150 (MQ=255)
                                                                                                                                                     tGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGt  >  2:48078/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGT  >  NC_000913/915078‑915375

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