Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 915,215 T→G 7.3% intergenic (‑358/+137) ybjE ← / ← aqpZ putative transporter/aquaporin Z

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913915,2150TG 7.3% 807.5 / 10.6 289intergenic (‑358/+137)ybjE/aqpZputative transporter/aquaporin Z
Reads supporting (aligned to +/- strand):  ref base T (200/60);  new base G (3/24);  total (203/86)
Fisher's exact test for biased strand distribution p-value = 1.53e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 298 total aligned reads displayed.

GGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTT  >  NC_000913/915066‑915364
                                                                                                                                                     |                                                                                                                                                      
ggggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTGACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAtt                                                                                                                                                        <  2:9655554/150‑1 (MQ=17)
 gggCAGCCCACCTGGCTGTGGGCTGAGGGCGCGTTTTACAGGGAGGTGCCCGCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTg                                                                                                                                                       <  2:5157311/150‑1 (MQ=11)
 gggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTg                                                                                                                                                       >  2:4599475/1‑150 (MQ=32)
  ggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTgg                                                                                                                                                      <  2:2138961/150‑1 (MQ=32)
  ggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTgg                                                                                                                                                      <  1:7270112/150‑1 (MQ=32)
  ggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTgg                                                                                                                                                      <  1:3893639/150‑1 (MQ=32)
  ggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTgg                                                                                                                                                      <  1:5337895/150‑1 (MQ=32)
   tcAGCCCACATCCATGTGGGCTGAGGGCGCCTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCGAGAGCCCGGGGGTGCAGGGGGCGGCGGCGGTTGgc                                                                                                                                                     >  2:7386399/2‑150 (MQ=11)
   gCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgc                                                                                                                                                     <  2:7749111/150‑1 (MQ=32)
   gCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgc                                                                                                                                                     >  1:1060053/1‑150 (MQ=32)
   gCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgc                                                                                                                                                     >  2:3665192/1‑150 (MQ=32)
     aGCCCACCTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgccg                                                                                                                                                   <  2:871427/150‑1 (MQ=17)
     aGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgccg                                                                                                                                                   <  1:4358003/150‑1 (MQ=32)
      gCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgccgc                                                                                                                                                  <  1:2330746/150‑1 (MQ=32)
       cccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgccgcc                                                                                                                                                 >  1:9495098/1‑150 (MQ=32)
        ccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGccc                                                                                                                                                >  2:3831570/1‑150 (MQ=32)
        ccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGccc                                                                                                                                                >  1:9565666/1‑150 (MQ=32)
         cacaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGcccc                                                                                                                                               <  2:1877614/150‑1 (MQ=32)
          acaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGccccc                                                                                                                                              >  2:8037552/1‑150 (MQ=32)
          acaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGccccc                                                                                                                                              >  2:1383278/1‑150 (MQ=32)
          acaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGccccc                                                                                                                                              >  2:2804165/1‑150 (MQ=32)
           caTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACTCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCt                                                                                                                                             <  2:1014956/150‑1 (MQ=17)
           caTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCt                                                                                                                                             <  2:3513747/150‑1 (MQ=32)
              ggATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTg‑cg                                                                                                                                         <  2:6108989/150‑1 (MQ=32)
              ggATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTg‑cg                                                                                                                                         >  1:8272202/1‑150 (MQ=32)
              ggATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTg‑cg                                                                                                                                         >  2:2883037/1‑150 (MQ=32)
               gATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTAGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTg‑cgc                                                                                                                                        >  2:590300/1‑150 (MQ=17)
                aTGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTg‑cgcg                                                                                                                                       >  2:3258484/1‑150 (MQ=32)
                   tgGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGAGTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGATTGGCCGCCCCCAG‑CGCGCTc                                                                                                                                    >  1:8935483/1‑150 (MQ=12)
                    gggCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTcc                                                                                                                                   >  2:5833563/1‑150 (MQ=32)
                    gggCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTcc                                                                                                                                   >  1:6655319/1‑150 (MQ=32)
                    gggCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGGCCGCCCCCTG‑CGCGCTcc                                                                                                                                   >  1:3084735/1‑150 (MQ=17)
                     ggCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCt                                                                                                                                  >  1:3254901/1‑150 (MQ=32)
                      gCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCtt                                                                                                                                 >  1:6222103/1‑150 (MQ=32)
                      gCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAG‑CGCGCGCCtt                                                                                                                                 >  1:6575179/1‑150 (MQ=17)
                       cTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTg                                                                                                                                >  2:1297023/1‑150 (MQ=32)
                          aGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTgcgc                                                                                                                             >  1:8049704/1‑150 (MQ=32)
                             gcgcGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAg                                                                                                                          >  2:650825/1‑150 (MQ=32)
                                cgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTgg                                                                                                                       >  2:4872065/1‑150 (MQ=32)
                                  ttttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCa                                                                                                                     >  2:6143699/1‑150 (MQ=34)
                                     tACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAata                                                                                                                  >  2:239631/1‑150 (MQ=37)
                                          ggATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCGTGCGCCAGTGGCAATATGTTg                                                                                                             >  1:7878046/1‑150 (MQ=255)
                                          ggATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTg                                                                                                             >  1:9139939/1‑150 (MQ=255)
                                               ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAg                                                                                                        >  2:377693/1‑150 (MQ=255)
                                               ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGCGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAg                                                                                                        >  2:5528502/1‑150 (MQ=255)
                                                  ccTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGGCCGCCCCCTG‑CGCGCGCCGTGCGCCCGTGGCAATATGTTGCTTAGCTc                                                                                                     >  1:2396303/1‑150 (MQ=255)
                                                    tCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGTTTTGGCGCCCCCAG‑CGCGCTCCTTTCGCCAGTGGCAAAATTTTGGTTATCTCAt                                                                                                   >  1:3422613/1‑150 (MQ=14)
                                                           cGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTTGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGa                                                                                            >  1:2843708/1‑150 (MQ=255)
                                                            gACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAAAATTTTGCTTATCTCATGGAAGGAg                                                                                           >  1:3512564/1‑150 (MQ=255)
                                                             aCCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCCGAGCCCGGGGGTGGCGGGGGCGGCGGCGCTGGGCCGCCCCCAG‑CGCGCTCCCTGCGCCAGGGGCAAAAAGTTTCTTAGCTCAAGGAAGGAgc                                                                                          >  1:4134772/1‑150 (MQ=14)
                                                             aCCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAgc                                                                                          >  2:9237385/1‑150 (MQ=255)
                                                             aCCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAG‑CGCGCTCCGTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGgcgc                                                                                          >  1:696644/1‑150 (MQ=255)
                                                                cGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTAGCTTAGCTCATGAAAGGAGCGca                                                                                       >  2:8073445/1‑150 (MQ=255)
                                                                      ccGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTTGCCCCCCCCAG‑CGCCCGACAATCGCCACTGGCAATAAGTTTCGTAGCTCATGAAAGGAgcgccaacaga                                                                                 >  1:4099581/1‑143 (MQ=14)
                                                                       cGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGATTGGCCGCCCCCAG‑CGCGCTCCATGCCCCAGTGGCAAAAAGTAGCTTAGCACAGGAAAGGAGCGCAACAAgat                                                                                >  2:2553902/1‑150 (MQ=17)
                                                                       cGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGAGTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgat                                                                                >  2:6899068/1‑150 (MQ=255)
                                                                            gggATATGTCGTGGGCACCGCGCAGGGGCGTTTTTGTTCGCCTGAGCCCGGGGGTGCTGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGAtgatga                                                                           <  2:7305412/150‑1 (MQ=18)
                                                                             ggATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCCCCCCCAG‑CGCGCGCCGTGCGCCCCGGGCAATATGTTTCTTAGCTCATGAAAGGAGCGCAAAAAgaatattaa                                                                          >  1:6710089/1‑143 (MQ=17)
                                                                               aTAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCCG‑CGCGCTCCATTCGCCCGTGGCAAAATGTTGCGTAGCTCAAGAAAGGCGCGCAAAAAGatgatgtaac                                                                        >  1:8707834/1‑146 (MQ=17)
                                                                                 ttgtggtgggCACCGGGCCGGGTCTTTTTTGTTCGCCTGCGCCGGGGGGTGCTGGGGGGGGGGGGTTTGGGCCGCCCGCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCaa                                                                      <  2:1357838/143‑1 (MQ=16)
                                                                                  tgTCGTGGGCACCGCGCCGGGGCTATTTTGTGCGCCAGCGCCGGGGGGTGCTGGGGGGGGCGGCTTTTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAAc                                                                     <  2:702796/149‑1 (MQ=18)
                                                                                   gTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACa                                                                    >  2:5669050/1‑150 (MQ=255)
                                                                                    tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCCGAGCCCGGGGGTGCCGGGGGCGGCGGCGATTGGCCGCCCCCGG‑TGCGCGCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACat                                                                   >  1:2891795/1‑150 (MQ=18)
                                                                                    tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACat                                                                   >  2:768296/1‑150 (MQ=255)
                                                                                    tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACat                                                                   >  2:3013759/1‑150 (MQ=255)
                                                                                    tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACat                                                                   >  2:4335500/1‑150 (MQ=255)
                                                                                      gtgggCACGGCGCTGCGGTTTTTTTGTTCGCCTGTGCGGGGGGGTGGTGGGGGGGGGGGGTTTTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACatat                                                                 <  2:2307745/148‑1 (MQ=21)
                                                                                      gtgggCACGGCGCAGCGGCTTTTTTGTTCGCCTGAGCCCGGGGGTGCTGGGGGCGGCGGGGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACatat                                                                 <  2:6431759/148‑1 (MQ=21)
                                                                                      gAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAG‑CGCGCTCCGTGCGCCACTGGCAAAATTTTTCTTATCTCATGAAAGGAGCGCAACAAGATGATGAATCAACatat                                                                 >  1:519907/1‑150 (MQ=18)
                                                                                         ggCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGGGGCGGCGCTTTGGCGCCCCCAG‑CCCCCGCCATGCGCCCGGGGCAAAAATTTTATTAGCTCAAGAAAGGAGCGCAACAAAAATATGAAAAAACaaaaaaa                                                              >  1:5777564/1‑144 (MQ=9)
                                                                                              cgcgCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACAt                                                         >  2:6691360/1‑150 (MQ=255)
                                                                                                        tttttttGTTCGCCTGCGCCGGGGGGGGCGGGGGGGGGCGGCTTGGGCCCCCCCCCCC‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaaaaa                                               <  1:4759858/147‑1 (MQ=255)
                                                                                                        tttttttGTTCCCCCGCGCGGGGGGGTGGGGGGGGGGGGGGCTTGCGCCCCCCCCCCC‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaaaaa                                               <  1:9461568/147‑1 (MQ=255)
                                                                                                        gttttttttTCGCCTGTGCGGGGGGGGGCTGGGGGCGGGGGC‑TTTGGCCCCCCCCCGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaaaaa                                               <  1:755663/147‑1 (MQ=255)
                                                                                                         gATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTTGCCCCCCCCAG‑CGCGCTCCGTGCGCCCCTGGCAAAAAGTTTCTTATCTCATGAAAGGAGAGCAACAAGAATATTAATAAACAAATAAAAACATCTTAAAAAAAg                                              >  1:5588528/1‑150 (MQ=18)
                                                                                                         gATTATGTTCGCCAGAGCCCGGGGGAGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGAGCCAGTGGCAATATGTTGCTAAGCACATGAAACGAGCGCAACAAGATGAAGAATCAACATATAACAACAGCTAAAAAAAAg                                              >  2:6663831/1‑150 (MQ=18)
                                                                                                                    ccAGAGCCCGGGGGTGCAGGGGGCGGCGGCGTTGGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTa                                   <  2:8753807/150‑1 (MQ=255)
                                                                                                                    ccAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTa                                   >  2:1900381/1‑150 (MQ=255)
                                                                                                                     cAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTAc                                  >  2:6774952/1‑150 (MQ=255)
                                                                                                                      agagCCCGGGGGTGCAGGGGGCGGCGGCGATGGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACAGATAACAACATCTTAAAAAAAGGCCTGACATTACg                                 >  2:8844760/1‑150 (MQ=255)
                                                                                                                        cgCCGGGGGGTGCGGGGGGGGGGGGGGTGGGGCGGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGcc                               <  2:8577109/149‑1 (MQ=38)
                                                                                                                          cccGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAg                             >  2:8445842/1‑150 (MQ=255)
                                                                                                                            cGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGAAATTACGCCAGGc                           >  2:51392/1‑150 (MQ=255)
                                                                                                                                 gTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATAACGCCAGGCCTTCt                      >  2:5476520/1‑150 (MQ=255)
                                                                                                                                  tGCTGGGGGGGGGGGCGTTTGGCCGCCCCGTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTg                     <  2:6223597/150‑1 (MQ=38)
                                                                                                                                    cggggggcgggggcGATGGCCCCCCCCCGC‑CGCGCTCCGTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCg                   <  1:8240759/148‑1 (MQ=255)
                                                                                                                                     aGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAATGCCTGACATTACGCCAGGCCTTCTGCGt                  >  1:1203334/1‑150 (MQ=255)
                                                                                                                                      gggggCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtt                 >  1:7643322/1‑150 (MQ=255)
                                                                                                                                       ggggCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCGTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtta                >  1:2668144/1‑150 (MQ=255)
                                                                                                                                          ggggcggccAATCGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaatt             >  2:316293/3‑150 (MQ=255)
                                                                                                                                           cggcggcgATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaatta            >  1:709248/1‑150 (MQ=255)
                                                                                                                                             ggggcttttGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTtaattaat          <  2:4057885/143‑1 (MQ=255)
                                                                                                                                             ggggcgTTTGCCCCCCCCCGC‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTtaattaat          <  2:5958475/148‑1 (MQ=255)
                                                                                                                                             gcggcgATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCCGGCCTTCTGCGttaattact          >  2:5204186/1‑148 (MQ=255)
                                                                                                                                               ggcgATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCa        >  2:3878358/1‑150 (MQ=255)
                                                                                                                                               ggcgATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCa        >  1:5997755/1‑150 (MQ=255)
                                                                                                                                               ggcgATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGGGAAGTAATCa        >  1:1741645/1‑150 (MQ=255)
                                                                                                                                                     tGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCttt  <  2:6320003/150‑1 (MQ=255)
                                                                                                                                                     tGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCttt  <  2:2436787/150‑1 (MQ=255)
                                                                                                                                                     tGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCttt  <  2:1430387/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTG‑CGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTT  >  NC_000913/915066‑915364

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: ‑