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.

GGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTT  >  NC_000913/915066‑915364
                                                                                                                                                       |                                                                                                                                                     
ggggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAtt                                                                                                                                                       <  1:816219/150‑1 (MQ=32)
 gggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTg                                                                                                                                                      >  2:2726425/1‑150 (MQ=32)
 gggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTATTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTg                                                                                                                                                      <  1:3150627/150‑1 (MQ=17)
 gggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTATTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTg                                                                                                                                                      <  1:3150700/150‑1 (MQ=17)
    cAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgcc                                                                                                                                                   <  1:2690812/150‑1 (MQ=32)
      gCCCACATCCATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgccgc                                                                                                                                                 >  1:2549950/1‑150 (MQ=17)
       cccACATCCATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgccgcc                                                                                                                                                >  2:3116808/1‑150 (MQ=17)
        ccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGccc                                                                                                                                               <  1:975049/150‑1 (MQ=32)
        ccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGccc                                                                                                                                               >  2:645505/1‑150 (MQ=32)
        ccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTGACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGccc                                                                                                                                               <  2:1620137/150‑1 (MQ=17)
            aTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTg                                                                                                                                           >  1:3038868/1‑150 (MQ=32)
              ggATGTGGGCTGAGGGCGCGTTTTACAGGGTTGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTgtg                                                                                                                                         >  2:2909736/1‑148 (MQ=17)
              ggATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTgtg                                                                                                                                         >  2:519689/1‑148 (MQ=17)
              ggATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTgtg                                                                                                                                         >  2:2909732/1‑148 (MQ=17)
              ccaTGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTgcg                                                                                                                                         >  1:1851452/3‑150 (MQ=17)
                     ggCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCt                                                                                                                                  >  2:1400793/1‑150 (MQ=32)
                     ggCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGGTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCt                                                                                                                                  >  2:451947/1‑150 (MQ=17)
                             gcgcGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAg                                                                                                                          >  2:1574939/1‑150 (MQ=32)
                              cggtttttACGGGTTTGTTACCCCGCGCCCGA‑CCGGTTCGCGTTATGGTTAATTTGGCGGCCCCCGCGCCGCGGAT‑TTTTTGTTGGCCAGAGCCGGGGGGTGATGGGGGCGGGGCCGTTTGGCCCCCCCCCGCGCGCTCCTTGCGCCAGt                                                                                                                         <  2:1974729/146‑1 (MQ=2)
                              cgcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGt                                                                                                                         >  1:1652930/1‑150 (MQ=32)
                              cgcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCTTGCGCCAGt                                                                                                                         >  2:2855305/1‑150 (MQ=17)
                              cgcgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCATGCGCCAGt                                                                                                                         >  1:2611496/1‑150 (MQ=17)
                                 gTTTTACTGGGTTGTTACCTCGCGCCCGCCCGGGTAGCCGTATGGGATAC‑GTGGTGGGCACCGCGCCGCGGCG‑ATTTTGTTCGCCTGCGCCCGGGGGTGCTGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGc                                                                                                                      <  2:2986580/150‑1 (MQ=11)
                                  ttttACTGGTTTGTTACCTCGCGCCCGACCGGGTAGCCGTTAGGTATAC‑GTCGTGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCTGAGCCGGGGGGTGCTGGGGGCGGCGGCGTTGGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCa                                                                                                                     <  2:2679756/150‑1 (MQ=11)
                                    ttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAat                                                                                                                   >  1:198849/1‑150 (MQ=35)
                                        aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGt                                                                                                               >  2:2448849/1‑150 (MQ=39)
                                        aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGt                                                                                                               >  2:1676758/1‑150 (MQ=39)
                                        aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGt                                                                                                               >  2:1676738/1‑150 (MQ=39)
                                        aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGt                                                                                                               >  2:1297345/1‑150 (MQ=39)
                                         gggATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGtt                                                                                                              >  2:2294866/1‑150 (MQ=40)
                                         gggATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAATGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGtt                                                                                                              >  2:1200112/1‑150 (MQ=255)
                                               ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑CTTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCGCCCCCTGCGCGCTCCGTGCGCCAGTGGCAAGATGTTGCTTAg                                                                                                        >  2:2967226/1‑150 (MQ=255)
                                               ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAg                                                                                                        >  2:804777/1‑150 (MQ=255)
                                               ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAAAATGTTGCTTAg                                                                                                        >  1:3056633/1‑150 (MQ=255)
                                                    tCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCAt                                                                                                   >  2:699917/1‑150 (MQ=255)
                                                    ccgcgCCCGACCGGGTAGCCGTAGGGGATAT‑GTGGTGGGCCCCGCGCCGGGGCGTTTTTTGTCGGCCAGCGCCGGGGGGTGCTGGGGGGGGGGGGGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCa                                                                                                    <  2:245425/149‑1 (MQ=12)
                                                          ccGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCGTGCGCCAGGGGCAATATGTTGCTTAGCTCATGAAAgg                                                                                             >  1:1136070/1‑150 (MQ=255)
                                                          ccGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCGCCTTGCGCCCGTGGCAATATGTTGCTTAGCTCATgccagg                                                                                             >  1:18490/1‑145 (MQ=255)
                                                          ccGACACGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGACAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAgg                                                                                             >  2:2119505/1‑150 (MQ=255)
                                                            gACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAg                                                                                           >  2:2627003/1‑150 (MQ=255)
                                                            gACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAg                                                                                           >  2:387306/1‑150 (MQ=255)
                                                            gACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAg                                                                                           >  2:2437037/1‑150 (MQ=255)
                                                            gACCCGGTAGCCGTAAGGGAGAA‑TTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCAAGAGCCCGGGCGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCACTGGCCATATGTTGCTTAGCTCATGTAAGGAg                                                                                           >  2:1927330/1‑150 (MQ=17)
                                                             aCCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGATCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTAATGAAAGGAgc                                                                                          >  2:2851335/1‑150 (MQ=255)
                                                             aCCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCGCCCCCAGCGCGCGCCTTGCGCCCGTGGCAATATGTTGCGTAGCTCATGAAAGGAgc                                                                                          >  1:1250011/1‑150 (MQ=17)
                                                             aCCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTTCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAgc                                                                                          >  2:974465/1‑150 (MQ=255)
                                                             aCCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAgc                                                                                          >  2:1769879/1‑150 (MQ=255)
                                                             aCCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAgc                                                                                          >  2:2623255/1‑150 (MQ=255)
                                                             aCCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAgc                                                                                          >  2:2135203/1‑150 (MQ=255)
                                                             aCCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAgc                                                                                          >  2:732043/1‑150 (MQ=255)
                                                             aCCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAgc                                                                                          >  1:389065/1‑150 (MQ=255)
                                                             aCCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCACGGGGCGGCGGCGAGTGGCCGCCCCCTGCGCGCGCCATGCGCCAGTGGCAATATGTTTCTTAGCTCATGTAAGGTgc                                                                                          >  1:1099641/1‑150 (MQ=17)
                                                              cccGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑AGTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCCGCGGCCAGGGGCCGCCCCCGGCGCCCGCCCTGCGCCCGGGGCCAAATGTTGCGTAGCGCCTGCAAGGAgcg                                                                                         >  1:274241/1‑150 (MQ=14)
                                                                       cGTAGGGTATAC‑GTCGGGGGCACCGCGCTGGGGGT‑TTTTTGTTCGCCAGCGCCGGGGGGTGCTGGGGGGGGGGGGTTGGGGCCGCCCCGTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgat                                                                                <  2:2021200/150‑1 (MQ=16)
                                                                          aaGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCATGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCTCCTTTCGCCAGTGGCAAAATGTTGCTTAGCTCATGAAAGGAGCGCAACAAtatgat                                                                             >  2:49714/1‑150 (MQ=18)
                                                                          aaGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgatgat                                                                             >  2:2887306/1‑150 (MQ=255)
                                                                          aaGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑AATTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATGGGCCGCCCCCTGCGCGCTCCATGCGCCAGTGGCAAAATGTTGATTAGCTCAAGAAAGGAGCGCAACAAgatgat                                                                             >  2:1950171/1‑150 (MQ=18)
                                                                              gagaa‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGATTGGCCGCCCCCAGCGCGCGCCATGCGCCAGTGGCAATATGTTTCTTAGCTCATGAAAGGAGCGCACCAAGATGATGAAt                                                                         >  1:2811850/4‑150 (MQ=18)
                                                                              gATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCCTGAAAGGAGCGCAACAAGATGATGAAt                                                                         >  2:585477/1‑150 (MQ=255)
                                                                                  ttgtggtgggCACCGCGCTGCGGGT‑TTTTTTTTCGCCTGAGCCGGGGGGTGCTGGGGGGGGCGGCTTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCaa                                                                      <  2:534541/143‑1 (MQ=18)
                                                                                  acgTCGCGGGCACCGCGCAGGGGCG‑ATTTTGTTCGCCTGAGCCGGGGGGTGCGGGGGGCGGCGGCGTTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCaa                                                                      <  2:534519/148‑1 (MQ=18)
                                                                                     tCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACat                                                                   >  2:2197083/1‑150 (MQ=255)
                                                                                         gggCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCGTGGCCGCCCCCCGCGCGCGCCGTGCGCCCCGGGCCAAAATTTTCTTAGCTCATGAAAGGCGCCCAACAAAAAGATGAATCaaacaaaaa                                                               >  1:1728673/1‑143 (MQ=9)
                                                                                         gggCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATaa                                                               >  2:1900839/1‑150 (MQ=255)
                                                                                            cACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCCCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATaacaa                                                            <  1:761286/150‑1 (MQ=255)
                                                                                                    tggggcg‑TTTTTGTGCGCCAGAGCCGGGGGGTGCCGGGGGCGGCGGCGTTTGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaa                                                    <  2:1851452/147‑1 (MQ=25)
                                                                                                          g‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCACTATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAATAAAAg                                              >  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)
                                                                                                                                                       |                                                                                                                                                     
GGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCG‑ATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTT  >  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: ‑