Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,914,5530AG28.8% 125.8 / 12.0 66H216H (CAT→CACglmUfused N‑acetylglucosamine‑1‑phosphate uridyltransferase and glucosamine‑1‑phosphate acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (16/31);  new base G (19/0);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 2.43e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CACGCAACCGGTGCCAATTTTCACGCGATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCACGGGTGATACGTCCATAACCGGTCGGATCAT  >  NC_000913/3914310‑3914800
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                           
cACGCAACCGGTGCCAATTTTCACGCGATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTTTTCACGCCTTCTACTTTGCTTAAAAGTTGCGGGTTAACgg                                                                                                                                                                                                                                                      >  2:347374/1‑251 (MQ=255)
cACGCAACCGGTGCCAATTTTCACGCGATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGGTGAACgg                                                                                                                                                                                                                                                      >  1:3997/1‑251 (MQ=255)
     aaCCGGTGCCAATTTTCACGCGATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCCCGCCCTCTACTTCGCTTAAACGTTTCGGGTGAACGGcgacg                                                                                                                                                                                                                                                 >  2:448016/1‑251 (MQ=255)
      aCCGGTGCCAATTTTCACGCGATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCCCCCCTTCTACTTCGCTTAAAAGTTGGGGGTGGACGGcgacg                                                                                                                                                                                                                                                 >  2:143532/1‑250 (MQ=255)
            gCCAATTTTCACGCGATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTAGTCTCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCa                                                                                                                                                                                                                                           <  1:68884/250‑1 (MQ=255)
               aaTTTTCACGCGATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGCCGATTTCACGcc                                                                                                                                                                                                                                       <  2:57025/251‑1 (MQ=255)
                        gcgATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGTCGACGATTCCACGCCCTTCCTGAt                                                                                                                                                                                                                              <  2:34993/251‑1 (MQ=255)
                               ccGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTTTTTACGCCTTCTACTTCGCTTAAACGTTGCGGGTGAACGGCGACGATTTCACGCCCTTCCTGATACGcc                                                                                                                                                                                                                        >  1:424374/1‑249 (MQ=255)
                               ccGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTCTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCTCTTAAACGTTGCGGGTGAACGGGGACGATTTCACGCCCTTCCTTGTACGcccg                                                                                                                                                                                                                      >  2:315916/1‑249 (MQ=255)
                                cGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCGTCCTGATACGCCAgc                                                                                                                                                                                                                     <  1:129145/251‑1 (MQ=255)
                                           ttGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTTACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGCCGGTTTCACGCCCTTCCTGATACGCCCGCGCCATAATGtt                                                                                                                                                                                                          >  1:399646/1‑250 (MQ=255)
                                               ccTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGGTGAACGGCGACGATTTCACGCCCTTCCCGGTACCCCACCGCAATTATGTCgggg                                                                                                                                                                                                      >  1:436537/1‑249 (MQ=255)
                                                cTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCCTTTACTTCGCTTAAACGTTGCGGGTGAACTGCGGAGATTTTACGCCCTTCCTCATAAGCCCGCGCAATAATGTCgggg                                                                                                                                                                                                      >  1:52544/1‑248 (MQ=255)
                                                cTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGAGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTTTTCACGCCTTCTAACTCGCTTTAACGTTGCTGATGGACGTCGGCCATTTCACACCCTTCCTGCTACACCAGCGCaataaaggcggagg                                                                                                                                                                                                     >  1:242388/1‑242 (MQ=255)
                                                      gatAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATgtagt                                                                                                                                                                                               <  1:347374/251‑1 (MQ=255)
                                                         aaCGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTAtt                                                                                                                                                                                            <  2:125687/251‑1 (MQ=255)
                                                         aaCGTTAGTATCAATTTCAACATCGCGCACGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTTACACCTTTTTCTTCGCTTAAAAGATTGGGGTGGAAGGCGACGATTTCCCGCCCTTCCTGTAAAGCCCGCGCAAGACTGTCGGGGATGTAGTAtt                                                                                                                                                                                            >  2:152853/1‑251 (MQ=255)
                                                               aGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGccc                                                                                                                                                                                       <  2:353830/250‑1 (MQ=255)
                                                                gTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTTACGCCTTCTACTTTGCTTAAACGTTGCGGGTGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGGGATGTAGTATTCTCCCTg                                                                                                                                                                                     >  1:119670/1‑251 (MQ=255)
                                                                          cAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCCGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAAAGGCGACGATTTCACCCCCTTCCTGATACGCCAGCGCAAAAATGTCGGGGATGTAGTATTCGCCCTGAGCattattt                                                                                                                                                                           >  1:467699/1‑250 (MQ=255)
                                                                            aCATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTTTTCACGCCTTCTACTTCCCTTAAACGTTGCGGGTGTACGGCGGCGATTTTACGCCCTTCCTTATACCCCCGCGCAATAAATTCGGGGGTGTAGTAtttccccctagccttcttgtt                                                                                                                                                                         >  2:386337/1‑232 (MQ=255)
                                                                             cATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTAttgttg                                                                                                                                                                        <  1:107252/251‑1 (MQ=255)
                                                                                   gcCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGGTATTCACGCCCTCTACTTCTCCTTAAAGTTGTGGATGAACGGGGACGATTTTACGCCCTTCCTGATACGCCACCGCAATAAATTAGGTGGTGTAGTATTTTCCCCGAGCATTTTTGTTGGTCatc                                                                                                                                                                 >  2:433897/1‑249 (MQ=255)
                                                                                    cccGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATGTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCt                                                                                                                                                                <  2:228224/251‑1 (MQ=255)
                                                                                     ccGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTTTTTACGCCCTCTACTTTGCTTAAAAGTTTCGGGTGGACGGCGGCGATTTTACGCCCTTCCTGATACGCCAGCGCCATAATGTTGGGGGTGTTGTATTTGTCCTTAGCCTTATTGTTGGTCAACtt                                                                                                                                                               >  2:156230/1‑251 (MQ=255)
                                                                                              aGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCgccagcca                                                                                                                                                      <  1:81865/251‑1 (MQ=255)
                                                                                                      aCGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGATAAAAGCAGTTTTTCCGCCTGTTCGGACTGATAAACACGCTCCAGCCGGGAGAGTTGCAGGCGGTTATTCACGGCTTCTACTTCGCTTAAAAGGGGCGGGTGGAAGGGGGCGGTTTCACGCCCTTCCTGATACGCCAGGGGAATAATGTCCGGGGTGTAGGATTCGCCCCGAGCATTTTTGTTGGGCAGCTTCGCCCGCCcgcggttt                                                                                                                                               >  1:20249/1‑246 (MQ=255)
                                                                                                                 aaCGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCg                                                                                                                                  <  1:78447/251‑1 (MQ=255)
                                                                                                                  aCGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGt                                                                                                                                 <  2:307301/251‑1 (MQ=255)
                                                                                                                         ggATCGCGCAGCATAACGCCTGCTAACAGCAGTTGTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAAt                                                                                                                         <  2:436537/251‑1 (MQ=255)
                                                                                                                          gATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGGTTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGGGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCCGCCCGCGGTTCATATCTGCGCCGTTGGCCATc                                                                                                                        >  1:279780/1‑251 (MQ=255)
                                                                                                                           aTCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCa                                                                                                                       <  1:143532/251‑1 (MQ=255)
                                                                                                                               cgcAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTTACGCCTTCTACTTCGCTTAAACGTTGGGGGTGAACGGCGGCGATTTTACGCCCTTCCTGATACGCCAGCGCCATAATGTCGGGGGTGTAGTATTCGCCCCGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGGTTCATATCTGCGCCGTTGGCCATCAGaa                                                                                                                    >  1:145231/1‑250 (MQ=255)
                                                                                                                                 cagcaTAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTTTTTACGCCTTCTTCTTCGCCTAAACGGTGCGGGTGGACGGCGGAGATTTTACGCCCTTCCTTATAACCCCACGCCATAATTTCGGGGATGTAGTATTCGCCCTGAGCATTATTTTTTGTCAGCTTCGCCCACCCGCGGTTCATATTTTCGCCGCTGGccacccaacagc                                                                                                                  >  2:358034/1‑240 (MQ=255)
                                                                                                                                   gcaTAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCg                                                                                                               <  1:315916/251‑1 (MQ=255)
                                                                                                                                    caTAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAAAGGCGACGATTTCACGCCCTTCCTGGTACGCCCGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTTTTGTTGGTCCGCTTCGCCAGCCCGCGTTTCCTATCTGCGCCGTTGGCAATCAGaaatccgg                                                                                                              >  1:107067/1‑245 (MQ=255)
                                                                                                                                       aaCGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTTACGCCTTCTACTTCGCTTAAACGTTGCGGGTGAACGGCGACGATTTTACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGGGATGGAGTATTTGCCCTGAGCATTATTGTTGGTCAGCTTTGCCAGCCAGCGTTTTATATCTGCGCCCGTGGCAATCCGAATGCCGgtgg                                                                                                           >  1:193427/1‑250 (MQ=255)
                                                                                                                                                   cagcagTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGGTGCGGGTGAAAGGCGCCGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTTGGTGGTGTAGGATTCCCCCTTAGCCTTATTGTTTGTCAGCTTTGCCAGCCCGCGGTTCATATCTGCTCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCtta                                                                                                >  2:178953/1‑248 (MQ=255)
                                                                                                                                                        gTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTTCACCTCGCTTAAAAGTTGCGGATGAAAGGCGACGATTTCCATCCCTTCCTGATACGCCAGCGCAATAATGTCGGGGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAACGGTTCATATCTGCGCCGTTGGCAATACGAATACCGGTGTTTATCTCCTGAATCTg                                                                                           >  1:184981/1‑250 (MQ=255)
                                                                                                                                                           tttCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTGCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACg                                                                                        <  2:280062/250‑1 (MQ=255)
                                                                                                                                                                      tCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTTACCCCTTCTCCTTCGCCTAAACGTTGCGAGTGAACGGCGACGATTTCACCCCCCTCCTGATACCCCCCCGCCATAATGTCCGGGGTGGAGGATTTTCCCCGAGCATTATTTTTTGTACGCCTTCCCCGCCAGCGGTTTCTTTCTCCGCCGTTGGCCAATAGAATGCCGGTGGTGATCTTCTTAATCTCACgccgcccctcggg                                                                            >  2:235129/1‑240 (MQ=255)
                                                                                                                                                                      tCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCgg                                                                             <  1:149225/250‑1 (MQ=255)
                                                                                                                                                                          aCTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTTTTTACCCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGAAGATTTCACGCCCTTCCTGATACGCCCGCGCCATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTTTTTTTGGTCAGCTTCGCCCGCCCGCGTTTCCTATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGCAGCTGCTCGTCGGTGGc                                                                         >  1:19996/1‑250 (MQ=255)
                                                                                                                                                                                aaaCACGCGCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTtg                                                                  <  2:431019/251‑1 (MQ=255)
                                                                                                                                                                                       cTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTTTACTTCGCTTAAACGTTTCGGGTGAACGGCGACGATTTTACGCCCTTCCCGATAAGCCCGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCCGCCAGCGTTTCATATCTGCGCCGTTTGCAATCAGAAAACCGGGGTTGGTCTCCTGAATCTGACGCTGCTCGTCGGGGGGCTCTTTGTGCTcc                                                            >  2:253859/1‑249 (MQ=255)
                                                                                                                                                                                          cAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTcaacaa                                                        <  1:5397/251‑1 (MQ=255)
                                                                                                                                                                                            gACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATg                                                      <  1:66338/251‑1 (MQ=255)
                                                                                                                                                                                                    agTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAAc                                              <  1:155436/251‑1 (MQ=255)
                                                                                                                                                                                                        gCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACttt                                           <  2:133955/250‑1 (MQ=255)
                                                                                                                                                                                                        gCAGGCCGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACttt                                           <  2:90888/250‑1 (MQ=255)
                                                                                                                                                                                                             cGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGtt                                     <  1:386337/251‑1 (MQ=255)
                                                                                                                                                                                                               gTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGtttt                                   >  2:282220/1‑251 (MQ=255)
                                                                                                                                                                                                                ttATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCATGATACGCCAGCGCAATAATGTCGGAGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGGGCTCCACAATGCCGGTACCTTTGCCGttttt                                  >  2:252195/1‑250 (MQ=255)
                                                                                                                                                                                                                 tATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTc                                  <  1:32526/250‑1 (MQ=255)
                                                                                                                                                                                                                  aTTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCAc                                <  1:1178/251‑1 (MQ=255)
                                                                                                                                                                                                                         ccTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCACGGGTGAt                         >  1:174266/1‑251 (MQ=255)
                                                                                                                                                                                                                           ttCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCACGGGTGATAc                       >  2:351855/1‑251 (MQ=255)
                                                                                                                                                                                                                            tCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCACGGGTGATAc                       <  1:448016/250‑1 (MQ=255)
                                                                                                                                                                                                                               ccTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCACGGGTGATACGTcc                   <  1:8312/250‑1 (MQ=255)
                                                                                                                                                                                                                               aCTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCAGCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCACGGGTGATACGTc                    <  1:452145/250‑1 (MQ=255)
                                                                                                                                                                                                                                   cGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCACGGGTGATACCTCCATa                >  2:159414/1‑250 (MQ=255)
                                                                                                                                                                                                                                        tAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCCGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCACTCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCACGGGTGATACGTCCATAACCgg           >  2:384904/1‑251 (MQ=255)
                                                                                                                                                                                                                                        gaaaCGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATGCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCGTACGCTGCGCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTGTTCACGGGTGAGACGTCCATAACCg            <  2:20249/249‑1 (MQ=255)
                                                                                                                                                                                                                                          aaCGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGACAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAAAAATGCCGGTAAATTTGCCGTTTTCACGGGTGATACGTCCATAACCGGTc         >  1:418284/1‑251 (MQ=255)
                                                                                                                                                                                                                                               ttGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCACGGGTGATACGTCCATAACCGGTCGGat     >  1:33916/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  cGGATGAACGGCGCCGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGGTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTTCTCGTCGGTGGCATCTTTGTGCGCAACAATGCCGTTATCTTGGCCGTTTTCACGGGTGATACGTCCATAACCGGTCGGatcat  <  1:272371/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                           
CACGCAACCGGTGCCAATTTTCACGCGATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCACGGGTGATACGTCCATAACCGGTCGGATCAT  >  NC_000913/3914310‑3914800

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