Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 915,215 T→G 8.6% 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 8.6% 323.3 / 5.1 123intergenic (‑358/+137)ybjE/aqpZputative transporter/aquaporin Z
Reads supporting (aligned to +/- strand):  ref base T (87/23);  new base G (2/11);  total (89/34)
Fisher's exact test for biased strand distribution p-value = 9.53e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 129 total aligned reads displayed.

GGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTT  >  NC_000913/915066‑915364
                                                                                                                                                       |                                                                                                                                                     
ggggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAtt                                                                                                                                                       <  1:816219/150‑1 (MQ=32)
 gggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTg                                                                                                                                                      >  2:2726425/1‑150 (MQ=32)
 gggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTATTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTg                                                                                                                                                      <  1:3150627/150‑1 (MQ=17)
 gggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTATTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTg                                                                                                                                                      <  1:3150700/150‑1 (MQ=17)
    cAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgcc                                                                                                                                                   <  1:2690812/150‑1 (MQ=32)
      gCCCACATCCATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgccgc                                                                                                                                                 >  1:2549950/1‑150 (MQ=17)
       cccACATCCATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgccgcc                                                                                                                                                >  2:3116808/1‑150 (MQ=17)
        ccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGccc                                                                                                                                               <  1:975049/150‑1 (MQ=32)
        ccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGccc                                                                                                                                               >  2:645505/1‑150 (MQ=32)
        ccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTGACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGccc                                                                                                                                               <  2:1620137/150‑1 (MQ=17)
            aTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTg                                                                                                                                           >  1:3038868/1‑150 (MQ=32)
              ggATGTGGGCTGAGGGCGCGTTTTACAGGGTTGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTgtg                                                                                                                                         >  2:2909736/1‑148 (MQ=17)
              ggATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTgtg                                                                                                                                         >  2:519689/1‑148 (MQ=17)
              ggATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTgtg                                                                                                                                         >  2:2909732/1‑148 (MQ=17)
              ccaTGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTgcg                                                                                                                                         >  1:1851452/3‑150 (MQ=17)
                     ggCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCt                                                                                                                                  >  2:1400793/1‑150 (MQ=32)
                     ggCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGGTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCt                                                                                                                                  >  2:451947/1‑150 (MQ=17)
                             gcgcGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAg                                                                                                                          >  2:1574939/1‑150 (MQ=32)
                              cggtttttACGGGTTTGTTACCCCGCGCCCGACCGGTTCGCGTTATGGTTAATTTGGCGGCCCCCGCGCCGCGGATTTTTTGTTGGCCAGAGCCGGGGGGTGATGGGGGCGGGGCCGTTTGGCCCCCCCCCGCGCGCTCCTTGCGCCAGt                                                                                                                         <  2:1974729/146‑1 (MQ=2)
                              cgcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGt                                                                                                                         >  1:1652930/1‑150 (MQ=32)
                              cgcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCTTGCGCCAGt                                                                                                                         >  2:2855305/1‑150 (MQ=17)
                              cgcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCATGCGCCAGt                                                                                                                         >  1:2611496/1‑150 (MQ=17)
                                 gTTTTACTGGGTTGTTACCTCGCGCCCGCCCGGGTAGCCGTATGGGATACGTGGTGGGCACCGCGCCGCGGCGATTTTGTTCGCCTGCGCCCGGGGGTGCTGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGc                                                                                                                      <  2:2986580/150‑1 (MQ=11)
                                  ttttACTGGTTTGTTACCTCGCGCCCGACCGGGTAGCCGTTAGGTATACGTCGTGGGCACCGCGCAGCGGCGATTTTGTTCGCCTGAGCCGGGGGGTGCTGGGGGCGGCGGCGTTGGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCa                                                                                                                     <  2:2679756/150‑1 (MQ=11)
                                    ttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAat                                                                                                                   >  1:198849/1‑150 (MQ=35)
                                        aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGt                                                                                                               >  2:2448849/1‑150 (MQ=39)
                                        aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGt                                                                                                               >  2:1676758/1‑150 (MQ=39)
                                        aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGt                                                                                                               >  2:1676738/1‑150 (MQ=39)
                                        aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGt                                                                                                               >  2:1297345/1‑150 (MQ=39)
                                         gggATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGtt                                                                                                              >  2:2294866/1‑150 (MQ=40)
                                         gggATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAATGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGtt                                                                                                              >  2:1200112/1‑150 (MQ=255)
                                               ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGCTTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCGCCCCCTGCGCGCTCCGTGCGCCAGTGGCAAGATGTTGCTTAg                                                                                                        >  2:2967226/1‑150 (MQ=255)
                                               ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAg                                                                                                        >  2:804777/1‑150 (MQ=255)
                                               ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAAAATGTTGCTTAg                                                                                                        >  1:3056633/1‑150 (MQ=255)
                                                    tCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCAt                                                                                                   >  2:699917/1‑150 (MQ=255)
                                                    ccgcgCCCGACCGGGTAGCCGTAGGGGATATGTGGTGGGCCCCGCGCCGGGGCGTTTTTTGTCGGCCAGCGCCGGGGGGTGCTGGGGGGGGGGGGGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCa                                                                                                    <  2:245425/149‑1 (MQ=12)
                                                          ccGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCGTGCGCCAGGGGCAATATGTTGCTTAGCTCATGAAAgg                                                                                             >  1:1136070/1‑150 (MQ=255)
                                                          ccGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCGCCTTGCGCCCGTGGCAATATGTTGCTTAGCTCATgccagg                                                                                             >  1:18490/1‑145 (MQ=255)
                                                          ccGACACGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGACAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAgg                                                                                             >  2:2119505/1‑150 (MQ=255)
                                                            gACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAg                                                                                           >  2:2627003/1‑150 (MQ=255)
                                                            gACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAg                                                                                           >  2:387306/1‑150 (MQ=255)
                                                            gACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAg                                                                                           >  2:2437037/1‑150 (MQ=255)
                                                            gACCCGGTAGCCGTAAGGGAGAATTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCAAGAGCCCGGGCGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCACTGGCCATATGTTGCTTAGCTCATGTAAGGAg                                                                                           >  2:1927330/1‑150 (MQ=17)
                                                             aCCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGATCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTAATGAAAGGAgc                                                                                          >  2:2851335/1‑150 (MQ=255)
                                                             aCCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCGCCCCCAGCGCGCGCCTTGCGCCCGTGGCAATATGTTGCGTAGCTCATGAAAGGAgc                                                                                          >  1:1250011/1‑150 (MQ=17)
                                                             aCCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTTCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAgc                                                                                          >  2:974465/1‑150 (MQ=255)
                                                             aCCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAgc                                                                                          >  2:1769879/1‑150 (MQ=255)
                                                             aCCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAgc                                                                                          >  2:2623255/1‑150 (MQ=255)
                                                             aCCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAgc                                                                                          >  2:2135203/1‑150 (MQ=255)
                                                             aCCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAgc                                                                                          >  2:732043/1‑150 (MQ=255)
                                                             aCCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAgc                                                                                          >  1:389065/1‑150 (MQ=255)
                                                             aCCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCACGGGGCGGCGGCGAGTGGCCGCCCCCTGCGCGCGCCATGCGCCAGTGGCAATATGTTTCTTAGCTCATGTAAGGTgc                                                                                          >  1:1099641/1‑150 (MQ=17)
                                                              cccGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGAGTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCCGCGGCCAGGGGCCGCCCCCGGCGCCCGCCCTGCGCCCGGGGCCAAATGTTGCGTAGCGCCTGCAAGGAgcg                                                                                         >  1:274241/1‑150 (MQ=14)
                                                                       cGTAGGGTATACGTCGGGGGCACCGCGCTGGGGGTTTTTTGTTCGCCAGCGCCGGGGGGTGCTGGGGGGGGGGGGTTGGGGCCGCCCCGTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgat                                                                                <  2:2021200/150‑1 (MQ=16)
                                                                          aaGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCATGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCTTTCGCCAGTGGCAAAATGTTGCTTAGCTCATGAAAGGAGCGCAACAAtatgat                                                                             >  2:49714/1‑150 (MQ=18)
                                                                          aaGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgatgat                                                                             >  2:2887306/1‑150 (MQ=255)
                                                                          aaGGGATAAGTCGAGGGCACCGCGCAGCGGCGAATTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATGGGCCGCCCCCTGCGCGCTCCATGCGCCAGTGGCAAAATGTTGATTAGCTCAAGAAAGGAGCGCAACAAgatgat                                                                             >  2:1950171/1‑150 (MQ=18)
                                                                              gagaaGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCGCCATGCGCCAGTGGCAATATGTTTCTTAGCTCATGAAAGGAGCGCACCAAGATGATGAAt                                                                         >  1:2811850/4‑150 (MQ=18)
                                                                              gATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCCTGAAAGGAGCGCAACAAGATGATGAAt                                                                         >  2:585477/1‑150 (MQ=255)
                                                                                  ttgtggtgggCACCGCGCTGCGGGTTTTTTTTTCGCCTGAGCCGGGGGGTGCTGGGGGGGGCGGCTTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCaa                                                                      <  2:534541/143‑1 (MQ=18)
                                                                                  acgTCGCGGGCACCGCGCAGGGGCGATTTTGTTCGCCTGAGCCGGGGGGTGCGGGGGGCGGCGGCGTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCaa                                                                      <  2:534519/148‑1 (MQ=18)
                                                                                     tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACat                                                                   >  2:2197083/1‑150 (MQ=255)
                                                                                         gggCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCGTGGCCGCCCCCCGCGCGCGCCGTGCGCCCCGGGCCAAAATTTTCTTAGCTCATGAAAGGCGCCCAACAAAAAGATGAATCaaacaaaaa                                                               >  1:1728673/1‑143 (MQ=9)
                                                                                         gggCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATaa                                                               >  2:1900839/1‑150 (MQ=255)
                                                                                            cACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCCCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATaacaa                                                            <  1:761286/150‑1 (MQ=255)
                                                                                                    tggggcgTTTTTGTGCGCCAGAGCCGGGGGGTGCCGGGGGCGGCGGCGTTTGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaa                                                    <  2:1851452/147‑1 (MQ=25)
                                                                                                          gATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCACTATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAATAAAAg                                              >  2:3142636/1‑150 (MQ=255)
                                                                                                          ctaTTTTGTTCGCCAGAGCCCGGGGGTGCTGGGGGCGGCGGCGTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaaaaa                                               <  2:826997/148‑1 (MQ=255)
                                                                                                             ttttGGTCGCCTGAGCCCGGGGGTGCCGGGGGCGGCGGCGATTGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGc                                            <  2:3209538/150‑1 (MQ=255)
                                                                                                                 gtccGCCCGCGCCGGGGGGGGCGGGGGGGGGCGGCTTGGGCCCCCCCCCCCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGa                                        <  1:1754130/147‑1 (MQ=255)
                                                                                                                  ttCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATAAACATATAACAACATCTTAAAAAAAGGCCTGAc                                       >  2:3044801/1‑150 (MQ=255)
                                                                                                                     gCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACAtt                                    >  2:155389/1‑150 (MQ=255)
                                                                                                                     gCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCAGACAtt                                    >  2:155415/1‑150 (MQ=255)
                                                                                                                     gCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAAAATGTTTCTTAGCTCATGAAAGGAGCGCAACAAGATGATTAATCAACATATAACAACATCTTAAAAAAATGCCTGACAtt                                    >  2:3003996/1‑150 (MQ=255)
                                                                                                                     gCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCTTGCGCCAGTGGCAAAATTTTTCTTAGCACATGAAAGGAGCGCAAAAAGATGATGAATCAAAATATAACAAAATCTTAAAAAAAGGCCTTACAtt                                    >  2:3003994/1‑150 (MQ=25)
                                                                                                                      ccAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCGTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATAAACATATAACAACATCTTAAAAAAAGGCCTGACATTa                                   >  2:885881/1‑150 (MQ=255)
                                                                                                                       ctgcgCCGGGGGGTGGGGGGGGGGGGGGCTTTGGCCCCCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTAc                                  <  2:409194/146‑1 (MQ=37)
                                                                                                                        agagCCCGGGGGTGCAGGGGGTGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACg                                 >  1:1634536/1‑150 (MQ=255)
                                                                                                                            cccGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAg                             >  2:241667/1‑150 (MQ=255)
                                                                                                                               gggggTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGcc                          >  2:2631293/1‑150 (MQ=255)
                                                                                                                               gggggTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGcc                          >  2:2570130/1‑150 (MQ=255)
                                                                                                                                 ggttGCTGGGGGGGGCGGGTTTTGGCCCCCCCGTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCtt                        <  2:447722/147‑1 (MQ=38)
                                                                                                                                 gggCGCCGGGGCCGCCCGCGAGGGCCCCCCCCCGCCGCGCTCCGTGCGCCAGGGGCAATATGTGGCTTAGCTCCTGAAAGGAGCGCAACAAGATGATGAATCAACATATAACACCATCTTAAAAAAAGGCCTGACATTACGCCAGGCCtt                        <  2:2211261/150‑1 (MQ=255)
                                                                                                                                  ggtgagggggggGCCGGTGATGGGCCCCCCCCCCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCGTAAAAAAAGGCCTGACATTACGCCAGTCCTTc                       <  2:2191134/144‑1 (MQ=255)
                                                                                                                                    tGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTg                     >  1:2356071/1‑150 (MQ=255)
                                                                                                                                    tGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTg                     >  1:2356054/1‑150 (MQ=255)
                                                                                                                                    tGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCTTGCGCCAGTGGCAAAATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAAAATATAACAACATCTTAAAAAAATGCCTGACATTACGCCAGGCCTTCTg                     >  1:601085/1‑150 (MQ=255)
                                                                                                                                       ggggggcgggggcGTTGGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGt                  <  2:1728699/149‑1 (MQ=255)
                                                                                                                                            ggggcggcgTTGGCCCCCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaatt             <  1:236175/148‑1 (MQ=255)
                                                                                                                                             gggcggcgATGGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGTCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaatta            <  2:190875/149‑1 (MQ=255)
                                                                                                                                               gcggcgTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTtaattaat          <  2:1641161/150‑1 (MQ=255)
                                                                                                                                               gcggcgATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTtaattaat          <  2:1862193/150‑1 (MQ=255)
                                                                                                                                                   cgATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACg      >  1:1064663/1‑150 (MQ=255)
                                                                                                                                                    ttttGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGc     <  2:1706386/148‑1 (MQ=255)
                                                                                                                                                    gATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGc     >  1:2924886/1‑150 (MQ=255)
                                                                                                                                                    gATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGc     >  1:2968008/1‑150 (MQ=255)
                                                                                                                                                    gATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGc     >  1:166468/1‑150 (MQ=255)
                                                                                                                                                    gATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGc     >  1:2968018/1‑150 (MQ=255)
                                                                                                                                                       tGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCttt  <  2:2565894/150‑1 (MQ=255)
                                                                                                                                                       tGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCttt  <  2:2565866/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
GGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTT  >  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: