Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 915,227 G→C 5.7% intergenic (‑370/+125) ybjE ← / ← aqpZ putative transporter/aquaporin Z

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913915,2270GC 5.7% 1012.0 / 10.2 331intergenic (‑370/+125)ybjE/aqpZputative transporter/aquaporin Z
Reads supporting (aligned to +/- strand):  ref base G (217/92);  new base C (3/19);  total (220/111)
Fisher's exact test for biased strand distribution p-value = 1.83e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 356 total aligned reads displayed.

ATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTC  >  NC_000913/915078‑915376
                                                                                                                                                     |                                                                                                                                                      
aTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTg                                                                                                                                                        <  2:9622822/150‑1 (MQ=32)
 tGGTTGTGGGCTGAGGGCGCGTTTTTCTGGGATGTTACCTCGCGCCCGACCCGTTAGCCTTAAGGTATAAGTCGAGGGCACCGCTCAGCGGCTTTTTTGTTCGCCAGAGCCCGGGGGTGCTGGGGGCGGCGGCGATTGGCCGCCCCCTgc                                                                                                                                                      <  2:5149336/150‑1 (MQ=11)
  ccaTGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAgcg                                                                                                                                                     >  1:7814652/3‑150 (MQ=17)
      gtgGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGATTGGCCGCCCCCGgcgcgcg                                                                                                                                                 >  1:4891324/1‑149 (MQ=17)
         ggCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCt                                                                                                                                              >  2:7509676/1‑150 (MQ=32)
          gCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCtt                                                                                                                                             >  2:92308/1‑150 (MQ=17)
              aGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTgcgc                                                                                                                                         >  2:562307/1‑150 (MQ=32)
                 gcgcGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAg                                                                                                                                      >  2:8133277/1‑150 (MQ=32)
                 gcgcGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAg                                                                                                                                      >  2:150304/1‑150 (MQ=32)
                 gcgcGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCATGCGCCAg                                                                                                                                      >  1:9366437/1‑150 (MQ=17)
                  cgcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGATTTGCCGCCCCCAGCGCGCTCCATGCGCCAGt                                                                                                                                     >  1:9483342/1‑150 (MQ=255)
                  cgcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGt                                                                                                                                     >  2:10736102/1‑150 (MQ=32)
                     gTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGc                                                                                                                                  >  1:4413896/1‑150 (MQ=33)
                     gTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGc                                                                                                                                  >  1:5389230/1‑150 (MQ=33)
                     gTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGGGGc                                                                                                                                  >  2:9169954/1‑150 (MQ=17)
                      ttttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCACTGGCa                                                                                                                                 >  2:3090216/1‑150 (MQ=21)
                        ttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAat                                                                                                                               >  2:525169/1‑150 (MQ=35)
                            aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGt                                                                                                                           >  1:7604759/1‑150 (MQ=39)
                            aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGt                                                                                                                           >  2:7740592/1‑150 (MQ=39)
                              ggttGTTCCCGCGCGCCCGCCCGGGTCGCCTTATGGGATACGTCGCGGGCCCCGCGCCGGGGCTTTTTGGGGCGCCCGCGCCGGGGGGGGCCGGGGGGGGCGGCGTTTGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTg                                                                                                                         <  2:8885620/147‑1 (MQ=9)
                              ggATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTg                                                                                                                         >  2:8747944/1‑150 (MQ=40)
                               gATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGc                                                                                                                        >  2:9832778/1‑150 (MQ=255)
                                aTGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCt                                                                                                                       >  1:2454705/1‑150 (MQ=255)
                                   ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCGTGCGCCAGTGGCAATATGTTGCTTAg                                                                                                                    >  1:10422008/1‑150 (MQ=255)
                                   ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCGGCGCGCGCCATGCGCCAGTGGCAATATGTTGCTTAg                                                                                                                    >  1:19131/1‑150 (MQ=255)
                                       cTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCGCCCCCAGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCa                                                                                                                >  1:3557027/1‑150 (MQ=255)
                                       cTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCGTGCGCCAGTGGCAATATGTTGCTTAGCTCa                                                                                                                >  1:4573068/1‑150 (MQ=255)
                                       cTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCGGCGCGCTCCGTGCGCCAGTGGCAATATGTTGCTTAGCTCa                                                                                                                >  1:902808/1‑150 (MQ=255)
                                       cTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCGTGCGCCCGTGGCCATATGTTGCTTAtctcc                                                                                                                >  1:3992528/1‑149 (MQ=14)
                                        tCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCCCCCCCAGCGCGCTCCATGCGCCAGTGGCAAAATGTTGCTTATCTCAt                                                                                                               >  1:10214906/1‑150 (MQ=14)
                                            gcCCGACCCGGTAGCCGTATGGGTTACGTCGTGGGCACGGCGCCGCGGCTTTTTTGTTGGCCAGAGCCGGGGGGTGCTGGGGGCGGGGGCTTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGaaa                                                                                                           <  2:1854774/150‑1 (MQ=17)
                                             cccGACCCCGTAGCCGTAAGGGATAAGTCGAGGGCACCCCGCAGCGGCGATTTTGTTCGCCAGAGCCCCGGGGTGCAGGGGGCGGCGGCGATTGGCCCCCCCCTGCGCGCTCCTTGCGCCACTGGCAATATGTTTCATAGCTCATGAAAg                                                                                                          >  2:2073245/1‑150 (MQ=17)
                                                   ccGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAgcgc                                                                                                    >  2:4113454/1‑150 (MQ=255)
                                                     ggTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGGAGGGGGCGGCGGCGATTGGGCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGGTTTTGTCATGAAAGGAGCGcaa                                                                                                  >  1:2921288/1‑150 (MQ=17)
                                                            gTATGGGATAAGTCGAGGGCACCGCGCCGCGGCGATTTTGTTCGCCAGCGCCGGGGGGTGCTGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgatg                                                                                           <  2:6338828/150‑1 (MQ=18)
                                                            gTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCGGCGCGCGCCATGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGGGCGCAACAAcatg                                                                                           >  1:2166496/1‑150 (MQ=18)
                                                            gTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAATGGCCGCCCCCAGCGCGCTCCTTTCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAAAAAgaag                                                                                           >  2:8291080/1‑148 (MQ=255)
                                                             tAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAATGGCCGCCCCCACCGCCCTCCATGCGCCACTGGCAAAAAGTTGCTTAGCTCAAGAAAGGAGCCcaacaacaaga                                                                                          >  2:2825580/1‑145 (MQ=17)
                                                                      aGTCGAGGGCACCGCGCCGCGGCTATTTTGTTCGCCAGAGCCGGGGGGTGCTGGGGGCGGGGGCTATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAAc                                                                                 <  2:7333182/150‑1 (MQ=18)
                                                                       gTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGATCCCGGGGGAGCAGGGGGCGGCGGCGATTGGCAGCCCCCAGCGCGCTCCTTGCGCCAGGGGCAATAAGTTGCTAAGCTCATGAAAGGAGCGCAACAACATGATGAATCAACa                                                                                >  2:1577174/1‑150 (MQ=18)
                                                                       gTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCCCCCCCACCGCGCGCCCTGCGCCCCGGGCAATATGTTGCTTTTCTCAGGGAAGGGGCGCCACAAGATGATGGATAAACa                                                                                >  1:5551359/1‑150 (MQ=17)
                                                                        tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACat                                                                               >  2:10067517/1‑150 (MQ=255)
                                                                        tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACat                                                                               >  2:3606542/1‑150 (MQ=255)
                                                                        tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACat                                                                               >  2:479245/1‑150 (MQ=255)
                                                                        tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACaa                                                                               >  2:5515217/1‑149 (MQ=255)
                                                                         cGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACata                                                                              >  2:7867340/1‑150 (MQ=255)
                                                                          gAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCCGCGAGTGGCCCCCCCCGGCGCGCTCCTTGCGCCCCGGGCAATATTTTTCTTATCGCATGAAAGGAGCGCAACAAGATGTTGAATCAAAatat                                                                             >  1:7197946/1‑150 (MQ=17)
                                                                           aGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACAtata                                                                            >  2:1346958/1‑150 (MQ=255)
                                                                           aGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACAtata                                                                            >  2:6772391/1‑150 (MQ=255)
                                                                           aGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGGGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACAtata                                                                            >  2:4362213/1‑150 (MQ=255)
                                                                            gggCACCGCGCTGGGGCTTTTTTTTTCGCCTGTGCGGGGGGGTGCTGGGGGGGGGGGGTTTTGGCGGCCCCTTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATaa                                                                           <  2:5902556/150‑1 (MQ=255)
                                                                             ggCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCCCGCCCCATGCCCCCGGGGCAAAAAGTTGCGTATCTCAAGGAAGGAGCGCAACAACATGTTGAATCAACAAATaac                                                                          >  1:1691194/1‑150 (MQ=17)
                                                                                  cgcgCAGCGGCGAGGTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGGCCGCCCCCTGCGCGCTCCGTGCGCCAGTGGCAATATGGTGCTTAGCTCATGACAGGAGCGCAACAAGATGATGAATCAACATATAACAACAt                                                                     >  2:7012264/1‑150 (MQ=18)
                                                                                           gcgATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTTGCCCCCCCCAGCGCGCGCCGTGCGCCAGTGGCAAAATTTTGCTTAGCTAATGAAAGGAGCGCAACAATATTGTGAAAAAACAAATAAAAAAAAATTaaaaaa                                                            >  1:6149324/1‑150 (MQ=17)
                                                                                            cttttttGTGCGCCTGCGCCGGGGGGTGGTGGGGGGGGCGGCGCTGGCCCCCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaaaaa                                                           <  2:3404944/147‑1 (MQ=255)
                                                                                             gATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAg                                                          >  2:10168248/1‑150 (MQ=255)
                                                                                                tttGTGCGCCAGCGCCGGGGGGGGGCGGGGGGGGGGGCGAGTGCCCCCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGcc                                                       <  2:9382805/150‑1 (MQ=25)
                                                                                                    ttCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCCTGAACGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTGAAAAAAAGGCCTGAc                                                   >  2:7067246/1‑150 (MQ=255)
                                                                                                       gCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACAtt                                                >  2:10944860/1‑150 (MQ=255)
                                                                                                           gcgCCGGGGGGTGCCGGGGGCGGCGGGGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGc                                            <  2:4617519/148‑1 (MQ=255)
                                                                                                            tgcgggggggTGATGGGGGCGGCGGGGTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGcc                                           <  2:819926/145‑1 (MQ=38)
                                                                                                             gCCGGGGGTTGCTGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCa                                          <  2:625384/150‑1 (MQ=255)
                                                                                                             gCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCa                                          >  2:747305/1‑150 (MQ=255)
                                                                                                              gggggggtTGCGGGGGGGGGGGGCGATTGGCCCCCCCCTGCGCGCTCCTTGCGCCAGTGGAAATATGTTGCTTTGCTCATGAAAGGAGCGCAACAAGATTATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAg                                         <  2:7763144/147‑1 (MQ=37)
                                                                                                                     gggCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACAGTACGCCAGGCCTTCt                                  <  2:10063819/148‑1 (MQ=255)
                                                                                                                     gTGCTGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCt                                  <  2:648000/150‑1 (MQ=255)
                                                                                                                      ggggggggggggggGCGGTGGCCCCCCCCGGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTg                                 <  2:3149768/146‑1 (MQ=255)
                                                                                                                       ggggggggcgggggcGTGGGCCCCCCCCCCGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGc                                <  1:6683051/147‑1 (MQ=255)
                                                                                                                        gtggggggggggGCTTTTGCCCCCCCCCCCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCg                               <  2:10422008/148‑1 (MQ=255)
                                                                                                                          ggggggggCGGCTTTGGGCCCCCCCCCGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtt                             <  1:3664450/150‑1 (MQ=39)
                                                                                                                          gggggCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtt                             <  2:9938922/150‑1 (MQ=255)
                                                                                                                           ggggggggggCGTTGGCCCGCCCCCGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtta                            <  2:3256820/150‑1 (MQ=255)
                                                                                                                           ggggcgggggcgCTGGGCCGCCCCGCGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGACAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtta                            <  2:8722109/150‑1 (MQ=38)
                                                                                                                           ggggCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtta                            >  2:2311564/1‑150 (MQ=255)
                                                                                                                            gggCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaa                           >  2:2954243/1‑150 (MQ=255)
                                                                                                                            agggggcggccAATCGCCGCCGCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaa                           >  1:9938266/5‑150 (MQ=39)
                                                                                                                             gggggcggctTTGGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaat                          <  2:4517079/147‑1 (MQ=255)
                                                                                                                                    ggttttGGCCCGCCCCTGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCa                    <  2:10214906/146‑1 (MQ=255)
                                                                                                                                     atcgccgccgccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACg                  >  1:10033189/7‑150 (MQ=255)
                                                                                                                                      gtttGGCCCCCCCCCGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACg                  <  2:4920775/148‑1 (MQ=255)
                                                                                                                                      gATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGc                 >  1:8615959/1‑150 (MQ=255)
                                                                                                                                        ttGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCtt               >  2:9296867/1‑150 (MQ=255)
                                                                                                                                           gccgccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTc            >  2:10511765/1‑150 (MQ=255)
                                                                                                                                           gccgccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTc            <  2:10021353/150‑1 (MQ=255)
                                                                                                                                            cggccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTcc           <  2:6952580/148‑1 (MQ=255)
                                                                                                                                            ccgccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGc                                                                                                                   <  2:3843137/46‑1 (MQ=255)
                                                                                                                                            ccgccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGc                                                                                                                   >  1:10645130/1‑46 (MQ=255)
                                                                                                                                            ccgccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTcc           >  2:1989527/1‑150 (MQ=255)
                                                                                                                                            ccgccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTcc           >  1:625675/1‑150 (MQ=255)
                                                                                                                                            ccgccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTcc           >  2:2548551/1‑150 (MQ=255)
                                                                                                                                            ccgccCCCGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTcc           <  2:4827646/150‑1 (MQ=255)
                                                                                                                                            ccccccccGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTcc           <  1:8792049/147‑1 (MQ=255)
                                                                                                                                             cgccCCCCCCGCGCTCCTTGCGCCAGTCGCAATATGTAGCTTAGCACATGAAAGGACCGCAACAAGATCATGAAACAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCca          <  1:10067517/150‑1 (MQ=255)
                                                                                                                                             cccccccGGCGCGCTCCTTGCGCCAGTGGCACGATGTTGCTGAGCTCATGAAAGGAGAGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCca          <  1:10882805/148‑1 (MQ=255)
                                                                                                                                             cccccccGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCca          <  1:5126052/148‑1 (MQ=255)
                                                                                                                                                 cccTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCcagcag      <  2:10033189/150‑1 (MQ=255)
                                                                                                                                                  ccTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAAAAAGATGATGAATCAACATATAACAACATCTTAAAAAAATGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAgg     >  1:2634531/1‑150 (MQ=255)
                                                                                                                                                    tGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGt   >  2:8061787/1‑150 (MQ=255)
                                                                                                                                                     gcgcgcTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTc  <  1:1989527/150‑1 (MQ=255)
                                                                                                                                                     gcgcgcTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTc  >  2:3872199/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
ATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTC  >  NC_000913/915078‑915376

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: