Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 915,213 A→T 12.7% intergenic (‑356/+139) ybjE ← / ← aqpZ putative transporter/aquaporin Z

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913915,2130AT12.7% 266.5 / 11.3 123intergenic (‑356/+139)ybjE/aqpZputative transporter/aquaporin Z
Reads supporting (aligned to +/- strand):  ref base A (75/26);  new base T (2/15);  total (81/42)
Fisher's exact test for biased strand distribution p-value = 1.39e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 123 total aligned reads displayed.

CGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATT  >  NC_000913/915065‑915353
                                                                                                                                                     |                                                                                                                                              
cGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAt                                                                                                                                               <  2:3771773/150‑1 (MQ=32)
 ggggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATACGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAtt                                                                                                                                             <  2:1887485/150‑1 (MQ=17)
 ggggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAtt                                                                                                                                             <  2:112981/150‑1 (MQ=32)
  gggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTg                                                                                                                                            <  2:3234682/150‑1 (MQ=32)
  gggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTg                                                                                                                                            <  1:110212/150‑1 (MQ=32)
  gggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTg                                                                                                                                            <  2:1910374/150‑1 (MQ=32)
  gggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTg                                                                                                                                            <  2:1845883/150‑1 (MQ=32)
    gCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgc                                                                                                                                          <  1:217068/150‑1 (MQ=32)
    gCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTCACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgc                                                                                                                                          <  2:3577259/150‑1 (MQ=17)
       gCCCACCTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgccgc                                                                                                                                       <  2:185200/150‑1 (MQ=17)
       gCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgccgc                                                                                                                                       <  1:3249302/150‑1 (MQ=32)
       gCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgccgc                                                                                                                                       >  1:372529/1‑150 (MQ=32)
       gCCCACATCCATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgccgc                                                                                                                                       >  2:3322037/1‑150 (MQ=17)
       gCCCACAGGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGGTACCGCGCGCCCGACCCGGTAGCCGTAAGGGATACGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgccgc                                                                                                                                       <  1:1247452/150‑1 (MQ=255)
        cccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgccgcc                                                                                                                                      >  1:1598147/1‑150 (MQ=32)
        cccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgccgcc                                                                                                                                      >  1:2887995/1‑150 (MQ=32)
        cccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgccgcc                                                                                                                                      >  1:2888309/1‑150 (MQ=32)
          cacaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGcccc                                                                                                                                    <  2:3216345/150‑1 (MQ=32)
          cacaGGGTTGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGcccc                                                                                                                                    <  1:2092602/150‑1 (MQ=17)
               ggATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTgcg                                                                                                                              >  2:53672/1‑150 (MQ=32)
               ggATGTGGGCTGAGGGCGCGTTTTAAAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAATGGCCGCCCCCTgcg                                                                                                                              >  2:53595/1‑150 (MQ=17)
                       gCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTccat                                                                                                                      >  2:3012395/1‑148 (MQ=17)
                       gCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCtt                                                                                                                      >  1:2059347/1‑150 (MQ=32)
                       gCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTccgt                                                                                                                      >  1:3756487/1‑148 (MQ=17)
                       gCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCtt                                                                                                                      >  2:1455653/1‑150 (MQ=17)
                         tGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTgc                                                                                                                    >  1:3156206/1‑150 (MQ=32)
                                  tttttACTGGGTTGTTACCTCGCGCCCGACCGGGTAGCGTTATGGTATATGTCGTGGCCACCGCGCTGCGGCTTTTTTTTACGCAAGAGCCCGGGGGTGCGGGGGGCGGCGGCTTTGGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGc                                                                                                           <  2:579695/149‑1 (MQ=2)
                                  gTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGc                                                                                                           >  2:2579567/1‑150 (MQ=33)
                                   ttttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCGCCCCCTGCGCGCGCCATGCGCCCGTGGCa                                                                                                          >  1:338800/1‑150 (MQ=255)
                                   ttttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCa                                                                                                          >  1:2446567/1‑150 (MQ=34)
                                   ttttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCATGCGCCAGTGGCa                                                                                                          >  1:2739283/1‑150 (MQ=21)
                                    tttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCaa                                                                                                         >  2:775605/1‑150 (MQ=35)
                                        ctggtttGTTCCCTGGCGCCCGACCCGGTCGCCGTTGGGTTAACGTGGTGGGCACCGCGCAGCGGCGTTTTTGTTCGCCCGCGCCGGGGGGTGCCGGGGGCGGCGGCGTTGGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATg                                                                                                     <  2:2526672/144‑1 (MQ=2)
                                         aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGt                                                                                                    >  1:253428/1‑150 (MQ=39)
                                         aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTACATGCGCCAGTGGCAATATGt                                                                                                    >  1:1897433/1‑150 (MQ=255)
                                           ggATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGGGGCGTTTTGCCGCCCCCTGGGCGCTCCTTGCGCCAGTGGCAATATGTTg                                                                                                  >  2:1373214/1‑150 (MQ=255)
                                                ttACCTCGCGCCCGACCGGTTCGCGTTATGGTCAAAGTCGGGGGCACGGCGCCGCGGCTTTTTTTTTCGCCAGGGCCGGGGGGTGCCGGGGGGGTGGGCGTTTGGCCGCCCCCTCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAg                                                                                             <  2:1647926/150‑1 (MQ=9)
                                                ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTTTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGGGATTGGCCGCCCCCAGCGCGCTCCTTGCGCCAGTGGCAATAAGTTGCATAg                                                                                             >  1:3402687/1‑150 (MQ=14)
                                                ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTAGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTTGCCGCCCCCAGCCCGCACCCTGCGCCCGTGGCAATAATTTGCTTAg                                                                                             >  1:3402647/1‑150 (MQ=14)
                                                 tACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGc                                                                                            >  2:3735356/1‑150 (MQ=255)
                                                   ccTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTc                                                                                          >  2:728090/1‑150 (MQ=255)
                                                    cTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCa                                                                                         >  2:371296/1‑150 (MQ=255)
                                                    cTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCGTGCGCCAGTGGCAATATGTTGCTTAGCTCa                                                                                         >  2:371253/1‑150 (MQ=255)
                                                     tCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCag                                                                                        >  1:3879567/1‑149 (MQ=255)
                                                     tCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCCGTGGCAATATGTTGCTTAGCTCAt                                                                                        >  1:3879605/1‑150 (MQ=255)
                                                     tCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCGTGCGCCCGGGGCAATATGTTGCTTAGCTCAt                                                                                        >  1:1812638/1‑150 (MQ=255)
                                                     ccgcgCCCCACCCGGTTGCGGTACGGTTTTCGTGGGGGCCACCGCGCTGGGGGTTTTTTTTTCGCCAGCGCCGGGGGGGGCTGGGGGCGGGGGGTTTGGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCAt                                                                                        <  2:3805790/149‑1 (MQ=12)
                                                      cgcgCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATg                                                                                       >  2:3549649/1‑150 (MQ=255)
                                                      cgcgCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATg                                                                                       >  2:776079/1‑150 (MQ=255)
                                                          cccGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGGCCGCCCCCAGCGCGCTCCTTGCGCCCGTGGCAATATGTTGCTTTGCGCCTGAAAg                                                                                   >  1:3679929/1‑150 (MQ=14)
                                                            cGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGa                                                                                 >  1:2284225/1‑150 (MQ=255)
                                                             gACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAg                                                                                >  2:3050960/1‑150 (MQ=255)
                                                             gACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAg                                                                                >  2:3817314/1‑150 (MQ=255)
                                                             gACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCACCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAg                                                                                >  1:3775347/1‑150 (MQ=255)
                                                                  ggTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATGGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCACCGGAAAGGAGCGcaa                                                                           >  2:3697319/1‑150 (MQ=255)
                                                                   gTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGcaac                                                                          >  2:3581451/1‑150 (MQ=255)
                                                                       ccGTACGGTATACGTGGAGGGCACCGCGACGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAga                                                                      <  2:3604670/150‑1 (MQ=18)
                                                                         gTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCATGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgatg                                                                    >  2:1924740/1‑150 (MQ=255)
                                                                         gTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCCGCCATTGGCCCCCCCCTGCGCGCGCCTTTCGCCACTGGCAATAATTTTCTTATCTCATTCAAGGAGCGCAACAAgaat                                                                    >  1:146865/1‑148 (MQ=17)
                                                                          tAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGGAGGGGGCGGCGGCGTTTGGGCGCCCCCTGCGCGCTCCTTGCGGCAGTGGCAATTTGTTGGTTTGCTGATGTTAGGTGCGCAACAAgattt                                                                   >  2:105784/1‑148 (MQ=14)
                                                                          tAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCGCCTTGCGCCACTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgatga                                                                   >  1:3816009/1‑150 (MQ=255)
                                                                              ggATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGaa                                                               >  2:3529613/1‑150 (MQ=255)
                                                                              ggATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCGCCATGCGCCCCTGGCAAAAATTTGCTTAGCTAATGAAAGGAGCGCAACAAAATTATGaa                                                               >  1:2695989/1‑150 (MQ=17)
                                                                              ggATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCTGCGCGCTCCTTGCCCCAGTGGCAAAAAGTTTCTTATCTCATGAAAGGAGCGCAAAAAGAATATGAAt                                                              >  2:2835342/1‑150 (MQ=17)
                                                                                   acgccgtgggCACCGCGCCGCGGCTTTTTTTTTCGCCAGAGCCGGGGGGTGCTGGGGGCGGCGGCGTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCaa                                                           <  2:1997124/143‑1 (MQ=18)
                                                                                      tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCGCCCCCAGCGCCCGCCCTGCGCCCCTGGCAAAATGTTTCATATCTCAGGGAAGGAGCGCAACAAGAAGAATAATAAACaa                                                        >  1:3529964/1‑149 (MQ=14)
                                                                                      tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCCCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGAAGATGAATCAACat                                                        >  2:3015331/1‑150 (MQ=255)
                                                                                              aCCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATaacaac                                                >  2:813808/1‑150 (MQ=255)
                                                                                              aCCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCATGAGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCCACAAGATGATGAATCAACAAATaacaac                                                >  2:813836/1‑150 (MQ=255)
                                                                                                cgcgCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACAt                                              >  2:1735789/1‑150 (MQ=255)
                                                                                                     cgcggcgTTTTTGTTCGCCAGCGCCGGGGGGTGCTGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaa                                         <  2:3956706/149‑1 (MQ=255)
                                                                                                      ggggcttttttGTTCGCCTGAGCGGGGGGGGGGGGGGGGGGGCGGCTTTGGCCCGCCCCCCGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaa                                        <  1:3289126/143‑1 (MQ=255)
                                                                                                       gggcttttttttTCGCCAGAGCGGGGGGGTGCTGGGGGCGGCGGCTATGGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaa                                       <  2:3868980/144‑1 (MQ=25)
                                                                                                         gcgATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGACTCAACATATAACAACATCTTaaaaaa                                     >  2:928315/1‑150 (MQ=255)
                                                                                                          tttttttGTTCGCCAGTGCGGGGGGGTGCGGGGGGGGGGGGCTTTGGCCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaaaaa                                    <  2:2759414/147‑1 (MQ=255)
                                                                                                          cttttttGTTCGCCAGCGCCGGGGGGTGCTGGGGGGGGGGGCGATTGGCGGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaaaaa                                    <  2:2130842/147‑1 (MQ=25)
                                                                                                           gtttttGTTCGCCAGAGCCCGGGGGTGCTGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAg                                   <  2:2972112/148‑1 (MQ=255)
                                                                                                           gATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAg                                   >  2:122623/1‑150 (MQ=255)
                                                                                                           gATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAg                                   >  2:464054/1‑150 (MQ=255)
                                                                                                           gATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAg                                   >  2:122581/1‑150 (MQ=255)
                                                                                                             ttttttTCGCCTGCGCCGGGGGGTGCGGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGc                                 <  2:3920946/150‑1 (MQ=25)
                                                                                                             ttttGTTCGCCAGAGCCCGGGGGTGCAGGGGGTGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGc                                 <  2:883642/150‑1 (MQ=255)
                                                                                                             ttttGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACAGATAACAACATCTTAAAAAAAGGc                                 >  2:2743181/1‑150 (MQ=255)
                                                                                                                     gCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACAtt                         >  2:698874/1‑150 (MQ=255)
                                                                                                                      ccAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCATGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAAAATCTTAAAAAAAGGCCTGAAATTa                        >  2:1469379/1‑150 (MQ=255)
                                                                                                                       cAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATGGGCCGCCCCCTGCGCGCTCCGTGCGCCAGGGGCAATATGTTGCGTAGCTCATGACAGGAGCGCAACAAGATGAAGCAGCAACATATCACAACATCTTCCAAAAAGGCCTGAccgtcc                       >  2:3519878/1‑145 (MQ=18)
                                                                                                                       cAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGGCCGCCCCCTGCGCGCGCCTTGCGCCCGGGGCAAGATGGTGCGTCGCGCCGGAAAGGAGCGCACCAAGATGAGGACTCAACCTAGAACACCCGCTTACAAAAAGGCCTGACagtcc                       >  2:3519876/1‑146 (MQ=12)
                                                                                                                         gggggggggggTGGTGGGGGGGGTGGTTTTTGGCCGCCGACGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGATCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGc                     <  2:3402687/144‑1 (MQ=255)
                                                                                                                          ttgcggggggTGTGGGGGGGGGGGGGGTTGTGGCCCCCCCGGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGc                     <  2:3402647/145‑1 (MQ=255)
                                                                                                                          cgcgggggggTGGTGGGGGGGGGGGGTTTGGGCCCCCCCCTCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGcc                    <  2:3720223/145‑1 (MQ=255)
                                                                                                                             ccGGGTGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAgg                 >  2:2610340/1‑150 (MQ=255)
                                                                                                                              cGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCTTGCTCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGAAGATGAATCAACAAATAACAACAACTTAAAAAAAGGCCTGACATTACGCCAGGc                >  2:1436878/1‑150 (MQ=37)
                                                                                                                                 gggTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCtt             >  1:2391853/1‑150 (MQ=255)
                                                                                                                                     gggggggggggggGCGGTGGGCCCCCCCCGGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATGTAACAACATCTTAAAAAAAGGCCTGGCGTAATGTCAGGCCttttt          <  2:1159784/147‑4 (MQ=255)
                                                                                                                                      cAGGGGGCGGCGGCGAGTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCg        >  1:1497719/1‑150 (MQ=255)
                                                                                                                                        ggggggggCGGCGTTGGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtt      <  2:2887995/150‑1 (MQ=255)
                                                                                                                                          gggCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaa    >  1:3331629/1‑150 (MQ=255)
                                                                                                                                            ggggcggccAATCGCCGCCGCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaatt  >  2:251787/3‑150 (MQ=255)
                                                                                                                                            ggggcggccAATCGCCGCCGCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaatt  >  1:3530897/3‑150 (MQ=255)
                                                                                                                                            gcgggggcgATTGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaatt  <  2:1288777/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
CGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATT  >  NC_000913/915065‑915353

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: