Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,914,5230CT23.5% 122.0 / 11.2 68V226V (GTG→GTAglmUfused N‑acetylglucosamine‑1‑phosphate uridyltransferase and glucosamine‑1‑phosphate acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (14/38);  new base T (16/0);  total (30/38)
Fisher's exact test for biased strand distribution p-value = 9.90e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GATTTCGCAATCATCGCCAATCACGCTGTTTTTAATCACGCAACCGGTGCCAATTTTCACGCGATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCACGG  >  NC_000913/3914274‑3914772
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                  
gATTTCGCAATCATCGCCACTCACGCTGTTTTTAATCACGCAACCGGTGCCAATTTTCACGCGATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGGTATTcc                                                                                                                                                                                                                                                                   >  1:95325/1‑250 (MQ=255)
 aTTTCGCAATCATCGCCAATCACGCTGTTTTTAATCACGCAACCGGTGCCAATTTTCACGCGATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTcc                                                                                                                                                                                                                                                                   >  1:130090/1‑249 (MQ=255)
         atcatcGCCAATCACGCTGTTTTTAATCACGCAACCGGTGCCAATTTTCACGCGATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTa                                                                                                                                                                                                                                                        <  1:107330/251‑1 (MQ=255)
            atcGCCAATCACGCTGTTTTTAATCACGCCACCGGTGCCAATTTTCACGCGATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACtt                                                                                                                                                                                                                                                    <  2:347044/251‑1 (MQ=255)
            atcGCCAATCACGCTGTTTTTAATCACGCAACCGGTGCCAATTTTAACGCGATGACCGAGAGTCACGTTGCCCTCGATGATACCGTTAGACTCAATTTCAACATCGCGCCCGCGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGCTAAACACGCTCCAGACGGGAGCGTTGCAGGCGGTTATTCACCCCTTCTACtt                                                                                                                                                                                                                                                    >  1:403266/1‑251 (MQ=255)
               gCCAATCACGCTGTTTTTAATCACGCAACCGGTGCCCATTTTCACGCGATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTTTTTACGCCTTcctccttgc                                                                                                                                                                                                                                                  >  1:265122/1‑243 (MQ=255)
               gCCAATCACGCTGTTTTTAATCACGCAACCGGTGCCAATTTTCACGCGATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGGTATTTAAGCcttcttctttgc                                                                                                                                                                                                                                                 >  1:238293/1‑248 (MQ=255)
                     cACGCTGTTTTTAATCACGCAACCGGTGCCAATTTTCACGCGATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTTGAGCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCGTCTACTTCGCTTAAAc                                                                                                                                                                                                                                           <  2:361538/251‑1 (MQ=255)
                      aCGCTGTTTTTAATCACGCAACCGGTGCCAATTTTCACGCGATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCAAGCCTTCTACTTCGCTTTAACg                                                                                                                                                                                                                                          >  1:98914/1‑251 (MQ=255)
                       cGCTGTTTTTAATCACGCAACCGGTGCCAATTTTCACGCGATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCTCGCCTTCTACTTCGCTTAAACGt                                                                                                                                                                                                                                         >  1:326917/1‑251 (MQ=255)
                              tttAATCACGCAACCGGTGCCAATTTTCACGCGATGACCGAGAGTCAAGTTGCCCTCGATGATAACGTTAGTATCCATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGCTCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCCCTTAAACGTTGccg                                                                                                                                                                                                                                    >  1:378264/1‑247 (MQ=255)
                                tAATCACGCAACCGGTGCCAATTTTCACGCGATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGGTTTTTCCCCCTTCTTCTTTCCTTAAACGTTggcggtgg                                                                                                                                                                                                                                >  2:272643/1‑244 (MQ=255)
                                       gCAACCGGTGCCAATTTTCACGCGATGACCGAGAGTCACGTTGACCTCGATGATAACGCTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTTACCGCCTTTACTTcccccaaaacttgaggaagaacgggg                                                                                                                                                                                                                          >  1:50374/1‑225 (MQ=255)
                                         aaCCGGTGCCAATTTTCACGCGATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCCCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGcgccg                                                                                                                                                                                                                      >  1:327699/1‑248 (MQ=255)
                                                    aTTTTCACGCGATGACCGAGAGTCACGTTGCCCTCGATCATAACGTTAGTCTCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCCTAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTTTTCCACCCTTCTACCTTCCTTAAACGGTGCGGGTGAACGGGGGCGATTTTACGccc                                                                                                                                                                                                           >  1:291595/1‑251 (MQ=255)
                                                         cACGCGATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTTACGCCCTCTACTTAGCTTAAAAGTTGCGGATGAAAGGCGACGAATTTACGCCCtttcc                                                                                                                                                                                                      >  1:235506/1‑248 (MQ=255)
                                                              gATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACTCTCCAGACGGGAGAGTTGCAGGCTGTTATTCACTCCTTCTACTTCGCTTAATCGTGGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATAc                                                                                                                                                                                                 <  1:269331/251‑1 (MQ=255)
                                                                     gagaGTGACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGCGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACCCGCTCCAGACGGGAGAGTTGCAGGCGGATATTTAAGCCCTCTACTTCTCTTTAACGTTGTGGGTGGACGGCGGCGATTTCaacgccttctctatttgccagcg                                                                                                                                                                                          >  2:328257/1‑229 (MQ=255)
                                                                      agagTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAgcgc                                                                                                                                                                                         <  1:73321/251‑1 (MQ=255)
                                                                       gagTCTCGTTGCCCTCGATGATAACGTTAGTCTCAATTTCAACATCGCTCCCGTGAGTTAGCGTACCACGCAGATCAAAGCGCGCTGTATCGCGCAGCATAACGCCTGCAAACAGCAGGTTTTCAGCCGTTTCGGACTGATAAACACGCTCTAGCCTGGAGAGTTGCAGGCGTTTATTCTCGCCTTCTACTTCCCTTAATCGTTGCGGAAGAACGTCCACGATTTCACGCCCACCCTGAAATGCCAGCGCa                                                                                                                                                                                        <  2:292076/251‑1 (MQ=255)
                                                                                 gCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTTAAAGGTGCGGATGCACGGCGACGATTTTACGCCCCTCCTGATACGCCAGCGCCATAATGTCgg                                                                                                                                                                              >  1:259277/1‑251 (MQ=255)
                                                                                   ccTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCTAGCGCAATAATGTCGGTg                                                                                                                                                                            <  2:87697/251‑1 (MQ=255)
                                                                                      cGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTTACGCCCTCTACTTCTCTTTAAAGGTGCGGGTGAACGGGGACGATTTCACGCCCTTCCTGATATGCCCGCGCCATTATGTCgggggtg                                                                                                                                                                         >  2:290528/1‑246 (MQ=255)
                                                                                      cGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACTCCAGCGCAAGAATGTCGGTGa                                                                                                                                                                           <  1:336371/249‑1 (MQ=255)
                                                                                         tgatAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTTACGCCCTCTACTTCGCTTAAAAGTTTCGGGTGAAAGGCGAAGATTTTACGCCCTTCCTGGAAAGCCCGCGCAAAAATGTCGGTGATgg                                                                                                                                                                        >  1:238335/1‑248 (MQ=255)
                                                                                            tAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCCCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTTACTCCTCTTATTTTGCTTAAACGGTGCGGGTGAACGGCCAGGATTTCACGCCCTTCCTTATAGCCCATCGCCATAATGTCGGTgctgccttat                                                                                                                                                                   >  1:231625/1‑242 (MQ=255)
                                                                                             aaCGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTTACCCCTTCTCATTCGCTTAAAAGTTGGGGATGAAAGGCCAGGATTTTCCGACCCTCCTGTGACGCCAGGGCaataattttgggggtgtaggatt                                                                                                                                                                  >  1:253861/1‑234 (MQ=255)
                                                                                             aaCGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGGACCACGCAGAGCAAAACGCGCTGGATCGCGCAGCATAACGCTTGCTAACAGCAGTTTTTCAGCCTGTTCTGACTGATAAACACGCTCCAGACGGGAGAGTGGCAGGCGGTTATTCACGCCTTATACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCTTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTAtt                                                                                                                                                                  <  2:164997/251‑1 (MQ=255)
                                                                                              aCGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGGACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTc                                                                                                                                                                 <  2:238888/251‑1 (MQ=255)
                                                                                                    gTATCAATTTCAACAACGCGCCCGTGAGTTAGCTTACCACGCAGATCAAAAAGCGCTGGATCGCGCAGCATAACGCCTTCTAACAGCAGTTTTTCAGCCCGTGCGGACTGATGAACACGCACCAGACGGGATAGTTGCAGGCGGTTATGCGCGCCTTCTACTTCGGTTAAACGTTGCGGATGAACGGCGACGATTGCTCGTCCTTCCAGAGACGCCAGCGCAATAATGTCGGGGATGTAGTATTCGctccg                                                                                                                                                           <  2:343784/251‑5 (MQ=255)
                                                                                                         aaTTTCAACATCGCGCCCGTGAGTTAGCGTACCTCGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACGGATAAACACGCTCCAGACGGGAGAGGTGCAGGCGGTTATTCACGCCTTCTACTGCGCGTAAACGTTGCGGAGGAACGGCGACGATGTCACGCCCTGCCTGATACGCCAGCGCAATAAAGTCGGCGATGTAGGATTCGCCCTGAGCat                                                                                                                                                     <  2:177853/251‑1 (MQ=255)
                                                                                                          aTTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCatt                                                                                                                                                    <  1:53588/251‑1 (MQ=255)
                                                                                                          aTTTCAACATCGCGCCCGTGACTGAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCTCTTAAACGTTGCGGATGAACGGCGACGAGTTCACGCCCTGCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCatt                                                                                                                                                    <  2:244130/251‑1 (MQ=255)
                                                                                                          aTTTCAACAGCGCGCCCGTGAGTTAGCGTACCACGCAGAGCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGAGAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCatt                                                                                                                                                    <  2:327471/251‑1 (MQ=255)
                                                                                                              cAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAACACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACCCCTTCTAATTCGCTTTAACGTTGCCGATGAACCGCGACGGTTTCACGCCCTTCCTGCTACGCCCGCGCCATAATGTCGGTGGTGGTGTATTCGCCCTGAGCATTAttg                                                                                                                                                >  1:65854/1‑251 (MQ=255)
                                                                                                                aCATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTGCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTAttgtt                                                                                                                                              <  2:382054/251‑1 (MQ=255)
                                                                                                                  aTCGCGCCCGTGAGGTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGACTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTg                                                                                                                                                          <  1:117925/237‑1 (MQ=255)
                                                                                                                             gAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTTACGCCTTCTACTTCGCTTTAACGTTGGGGGTGGACGGCGACGATTTCACGCCCTTCCTGATACGCCCGCGCAATAATGTCGGGGGTGTAGGATTCGCCCCGAGGATTATTTTTTGGTAGCTTCgcc                                                                                                                              >  1:168590/1‑251 (MQ=255)
                                                                                                                              aGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTgcgcaa                                                                                                                             <  2:82213/251‑3 (MQ=255)
                                                                                                                                 tAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGcc                                                                                                                                                                                              <  1:423264/187‑1 (MQ=255)
                                                                                                                                    cGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTTACGCCTTCTAATTCGCTTAAACGTTGCGGGTGTACGGCGACGATTTTACGCCCTTCCTTATACGCCCGCGCAATAATGTCGGTGATGTAGGATTCGCCCCGAGCATTTTTGTTGGGCAGCttttcccaccacc                                                                                                                        >  1:322618/1‑240 (MQ=255)
                                                                                                                                     gTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTTTTTACGCCTTTTACTTCGCTTTAACGGTGCGGGTGAACGGCGACGCATTCACGCCCTTTCTGATACGCCCGCGCCATAATGTCGGGGATGTAGTATTTCCCCTGAGCATTTTTTTTGGTCAGCTTCGCCCGCCcgcg                                                                                                                      >  2:329281/1‑251 (MQ=255)
                                                                                                                                       aCCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGt                                                                                                                     <  1:362513/250‑1 (MQ=255)
                                                                                                                                        ccACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGttt                                                                                                                   <  1:272643/251‑1 (MQ=255)
                                                                                                                                            gCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTTACACCTTCTACTTCGCTTAAACGTTGCGGGTGAACGGCGACGATTTCACGCCCTTCCTGATACGCCCGCGCCAGAATTTCGGGGATGTAGGATTCTCCCTGAGCATTTTTTTTTGGTAGCTTCGCCAGCCAGCGTTTCatt                                                                                                               >  1:104375/1‑250 (MQ=255)
                                                                                                                                                  cAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTgc                                                                                                          <  1:370592/250‑1 (MQ=255)
                                                                                                                                                               aTCGCGCAGCATAACGCCTGCTAACAGCAGTTTGTCGGCCTGTTCGGACTGATAAACCCGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCACTTAAACGTTGCGGATGAACGGCTACGATTTCACGCCCGTCCTTATACGCCAGCTCAATAAGGTCGGTGATGTAGTATTCGCCCTGAGCATTATTTTTTATCACCTTCGCCAGCCAGCGTTTCATATCTTCGCCGTTGGCAATCa                                                                                           <  2:390647/251‑1 (MQ=255)
                                                                                                                                                                        caTAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCgg                                                                                  <  1:416863/251‑1 (MQ=255)
                                                                                                                                                                           aaCGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTGATTCCTCCCCTCTACTTCCCTTTAAAGTTGTGGGTGGACGGCGGCGATTTCTCGCCCCTCCTGATTCGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCCTGGCATTATTTTTGGTCAGCTTTGCCCGACCGCGTTTCCTATCTGCGCCCTTTGCGATCAGaaacgccgggg                                                                               >  2:71357/1‑242 (MQ=255)
                                                                                                                                                                                  gCTAACAGCAGTTTTTCAGCCTGTTCGGCCTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGCTATTCACGCCCTCTCCTTCGCTTAAACGTTGCGGGTGAACGGCGACGATTTCACCCCCCTCCTGATCCGCCAGCGCAATAATGTCGGTGGTGTAGTATTTCCCCTGCGCATTTTTTTTGGTCCGCTTCGCCCGCCCGCGTTTCCTATTTGCGCCGTTGGCAATCAGAATGCCCGTGTTGAtctc                                                                        >  1:348113/1‑251 (MQ=255)
                                                                                                                                                                                  gCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCCCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGAtctc                                                                        <  2:130090/251‑1 (MQ=255)
                                                                                                                                                                                  gCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTACAGACGGGAGAGTTGCAGGCGGTTTTTTACGCCCTCTACTTCCCTTAAACGTCGCCGGTGCAACGCGACCGTTTAACCCCCCTCCCGATACCCCAGAGCACAAATGTCGGGGATGTGGTATTCTCCCTGAGCGTTATTTTTGGTCACCTTCGTCCACCCGCGGTTTTTATCTGCGCCCCTGGGGATCAGAATGCCGGGGTTGctctt                                                                        >  1:46557/1‑250 (MQ=255)
                                                                                                                                                                                  gCTAACAGCAGGTTTTCAGCCTGTGCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGAtc                                                                          <  2:22060/249‑1 (MQ=255)
                                                                                                                                                                                        agcagTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTTACGCCTTCTACTTTGCTTTAAAGTTGCGGGTGAAAGGCGCCGCTTTCACCCCCTTCCTTATACGCCAGCGCAATAATGTCCGTGATGTAGTATTCGCCCCTAGCATTTTTGTTGGGTCGCCTCGCCAGCCAGCGTTTCATATCTTCGCCCGTTGGAATAAGAATGGCGGGGTTTATTTCCTggat                                                                  >  1:412524/1‑248 (MQ=255)
                                                                                                                                                                                         gcagTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGTGAGAGTTGCAGGCGGTTATTCCCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCTCGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCTCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATc                                                                 <  2:153109/251‑1 (MQ=255)
                                                                                                                                                                                          cagTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTTACACCTTTTTCTTTCCTTTAACGTTGCGGGTGAACGGCGACGATTTCACCCCCTTCCTGATACCCCCGCGCAATAATGTTGGGGGTGTAGTATTCGCCCTGCGCATTTTTGTTGGGCTGCCTTGCCCGCCAGGGGTTTTTATCTGCGCCGTTGGCAATTAGAATTCCGGGGGTGGTCTCCTGGATCt                                                                >  1:62084/1‑251 (MQ=255)
                                                                                                                                                                                              ttttCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGGTGGCAATCAGAATGCCGGTGTTGAGCTCCTGAATCTGAc                                                             <  2:62084/250‑1 (MQ=255)
                                                                                                                                                                                                        gTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCtcgt                                                    <  1:392040/249‑1 (MQ=255)
                                                                                                                                                                                                        gTTCGGACGTATAAACACGCTCCTGACGGGTGAGTTTCAGGCGGTTAGTCTCGCTTTCTACTTCGCTTAAACGTTGCGGATGAACGTCGTCGATTTCACGCCCTTCCTGATACGCCAGCGCAATAAGGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGTGCAGCTTCTCCAGCCAGCGTTTCATATCTGCGCAGTTGGCTATCAGAAGGCCTGTGTTGATCTCCTGAATCTGACGCTGCtcgtcg                                                  <  2:265122/251‑1 (MQ=255)
                                                                                                                                                                                                                 gATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATc                                          <  1:290528/250‑1 (MQ=255)
                                                                                                                                                                                                                  aTAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCCCGTCGGTGGCATCtt                                        <  1:15943/251‑1 (MQ=255)
                                                                                                                                                                                                                   tAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCTCCCTGAGCATTATTTTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCtcgtc                                                   <  1:29933/239‑1 (MQ=255)
                                                                                                                                                                                                                            tCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCCTCTACTTCGCTTAAAAGTTGCGGGTGAACGGCGACGATTTCACGCCCCTCCCGATACGCCCGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTTTTTGTTAGCTTTTCCAGCCCGCGCTTTATATCTGCGCCGCTGGCAATCAGAATGCCGGTGTTGTTCTCCTGAATCTGCCGCCGCTCGTCGGGGGCATCTTTGTGCTccac                              >  1:19697/1‑248 (MQ=255)
                                                                                                                                                                                                                                   gggAGAGTTGCAGTCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGGAGTATTCGCCCTGAGCATTAGTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGAGCTCCTGAAGCTGACGCTGCtc                                                      <  1:132499/220‑1 (MQ=255)
                                                                                                                                                                                                                                         gTTGCAGGCGGTTATTCACTCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTGCACGCCCTTCCTGATACGCCAGCGCAATAAGGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCATCTTCGCCAGCCAGCGTTTCATATCGGCGCCGTTGGCAATCAGAATGCTGGTGTTGATCTCCTGGATCTGACGCTGCTCTTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACt                 <  2:269975/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 cGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGtt         <  1:8126/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      tATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTc      <  2:322618/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      tATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTc      >  1:370082/1‑250 (MQ=255)
                                                                                                                                                                                                                                                      tATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTc      <  1:162520/250‑1 (MQ=255)
                                                                                                                                                                                                                                                         tCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCACgg  >  2:275017/1‑251 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                  
GATTTCGCAATCATCGCCAATCACGCTGTTTTTAATCACGCAACCGGTGCCAATTTTCACGCGATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCACGG  >  NC_000913/3914274‑3914772

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: