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

CAATCATCGCCAATCACGCTGTTTTTAATCACGCAACCGGTGCCAATTTTCACGCGATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCAC  >  NC_000913/3914281‑3914770
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                            
caaacatcGCCAAGCACGCTGGTTTTAAGCACGCAACCGGTGCCAATTTTCACGCGATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTc                                                                                                                                                                                                                                                     <  2:522436/247‑1 (MQ=255)
   tcatcGCCAATCACGCTGTTTTTAATCACGCAACCGGTGCCAATTTTCACGCGATGACCGCGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTTTTTACGCCctctcc                                                                                                                                                                                                                                                  >  1:354079/1‑249 (MQ=255)
         ccAATCACGCTGTTTTTAATCACGCAACCGGTGCCAATTTTCACGCGATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTTACGCCCTCTTCTTTGCt                                                                                                                                                                                                                                           >  1:310337/1‑251 (MQ=255)
           aaTCCCGCTTTTTTTAATCACGCAACCGGTGCCAATTTTCCCGCGATGACCTAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGTGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCTCTTa                                                                                                                                                                                                                                         <  2:480828/251‑1 (MQ=255)
              cACGCTGTTTTTAATCACGCAACCGGTGCCAATTTTCACGCGATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTAATTTGCTTaaa                                                                                                                                                                                                                                       >  1:385702/1‑250 (MQ=255)
               aCGCTGTTTTTAATCACGCAACCGGTGCCAATTTTCACGCGATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTTACCGCCTCTCCTTCGCtttaa                                                                                                                                                                                                                                       >  1:556806/1‑246 (MQ=255)
               aCGCTGTTTTTAATCACGCAACCGGTGCCAATTTTCACGCGATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACg                                                                                                                                                                                                                                     <  1:464414/251‑1 (MQ=255)
               aCGCTGTTTTTAATCACGCAACCGGTGCCAATTTTCACGCGATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACg                                                                                                                                                                                                                                     <  1:304303/251‑1 (MQ=255)
                  cTGTTTTTAATCACGCAACCGGTGCCAATTTTCACGCGATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTTTTCACGCCCTCTACCTCgttttaacagtt                                                                                                                                                                                                                                  >  2:39909/1‑240 (MQ=255)
                          aaTCACGCAACCGGTGCCAATTTTCACGCGATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGaa                                                                                                                                                                                                                          <  1:483433/251‑1 (MQ=255)
                                 cAACCGGTGCCAATTTTCACGCGATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTTACGCCTTCTACTTTCCTTAAACGTTGCGGGTGGACGGcgac                                                                                                                                                                                                                   >  1:407355/1‑251 (MQ=255)
                                                  cccgcgATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCt                                                                                                                                                                                                 <  2:310337/249‑1 (MQ=255)
                                                  cACGCGATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTTTTTTCGCCCTCTACTTTGCTTAAACATTGCGGGGGAAAGGCGGCGGTTTCACACCCTTCCt                                                                                                                                                                                                 >  2:368063/1‑251 (MQ=255)
                                                      cgATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCCCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCTCTTCCTGATa                                                                                                                                                                                             <  1:177749/251‑1 (MQ=255)
                                                              gagaGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTTTTTACGCCTTCTACTTCCCTTAAACGGTGCGGATGAACGGCGACGATTTTACGCCCTTCCTGATACGCCcgcg                                                                                                                                                                                    >  1:449952/1‑251 (MQ=255)
                                                              gagaGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAgcg                                                                                                                                                                                    <  1:351459/251‑1 (MQ=255)
                                                                     aCGTTGCCCTCGATTATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTATCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCaataat                                                                                                                                                                             <  2:339792/251‑1 (MQ=255)
                                                                         tGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACTCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACTGCGACGATGTCACGCCCTTCCTGATACGCCAGCGCAATAATGTc                                                                                                                                                                          <  2:577574/250‑1 (MQ=255)
                                                                          gCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGAGCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATTCGCCAGCGCAATAATGTCg                                                                                                                                                                        <  2:308691/250‑1 (MQ=255)
                                                                            ccTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTTACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTTACGCCCTTTCTGCTTCGCCAGCGCAATAATGTCGGTg                                                                                                                                                                     >  1:58893/1‑251 (MQ=255)
                                                                            ccTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGATTTTTACTACTTCTTCTTCGCTTAAAACTTGCGGATGAAAGGCCACCATTTTCCGCCCTTCCTgtttatccagacgaatcatgtcggtg                                                                                                                                                                       >  1:528809/1‑226 (MQ=255)
                                                                             cTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTTACGCCCTCTACTTTGCTTAAACATTTCGGAGGAACGGCGGCGATTTTCCACCCTTCCTGATaacccacccccaatatgttggggg                                                                                                                                                                      >  2:334070/1‑228 (MQ=255)
                                                                                gatgatAACGTTAGTACCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATgt                                                                                                                                                                 <  1:275036/251‑1 (MQ=255)
                                                                                 atgatAACGTTATTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATgta                                                                                                                                                                <  1:105393/251‑1 (MQ=255)
                                                                                    atAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACCCCTTCTACTTTGCTTAAACGTTTCGGGTGAACCGCGACGAATTCACCCCCTTCCTGCTAACCCAGCGCACTAATGTAGGTGgtgtt                                                                                                                                                                >  1:569796/1‑247 (MQ=255)
                                                                                      aaCGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCGTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTAt                                                                                                                                                            <  1:457324/250‑1 (MQ=255)
                                                                                      aaCGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAAAGTTGCGGGTGAAAGGGGACGATTTCACGCCCTTCCTGCAACGCCAGCGCAATAATGTCGGGGAGGTAGAAtt                                                                                                                                                           >  1:40492/1‑251 (MQ=255)
                                                                                        cGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTTACGCCATCTAATTCTCTTAAACGTTGTGGATTAAAGGCGGAGCATTCACGCCCTTCCTGTTACGCCAGCGCTATAATGTCGGTGAGGTAGGAtttc                                                                                                                                                         >  1:54596/1‑249 (MQ=255)
                                                                                        cGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGGTATTTCCTCCTTCTACTTCCCTTCAACGTTTGGGGTGAAAAGCGACGATTTTACGCCCTTCCTGGTAAGCCCGAGGAATAATGTTGGGGATGTGATATTCg                                                                                                                                                         >  1:37628/1‑251 (MQ=255)
                                                                                        cGTTAGTATCAAATTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGAGCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGTCGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCg                                                                                                                                                         <  1:290210/251‑1 (MQ=255)
                                                                                                 cAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTTACGCCTTCTTCTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTTACGCCCTTCCTGATACGCCAGCGCCATAATTTCGGTGGTGTAGTATTTGCCCTTAGcc                                                                                                                                                >  1:52117/1‑250 (MQ=255)
                                                                                                     ttCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTGTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCatta                                                                                                                                             <  2:243550/250‑1 (MQ=255)
                                                                                                      tCAACATCGAGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCattat                                                                                                                                            <  2:291428/250‑1 (MQ=255)
                                                                                                         aCATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTAttgtt                                                                                                                                        <  1:37255/251‑1 (MQ=255)
                                                                                                           aTCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAAAGTTGCGGGAGAACGGGGACGATTTTACCCCCTTCCTGATACGCCAGCGCCATAATGTCCGGGGTTGAGTATTCGCCCTgggcctttttttttg                                                                                                                                      >  1:490363/1‑237 (MQ=255)
                                                                                                              gctcccGTGAGTTAGCGTACCCCGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACAGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTAATTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGAGACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCa                                                                                                                                   <  2:446798/248‑1 (MQ=255)
                                                                                                                   cgtgatgtcgCGTACCACTCAGATGCAAACGCGCTGTATCGCGCAACATAACGCCTGCTAATAGCAGTTTTTCAGCCTGTTCGTACGTATAAACACGTTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTCATCGTTGCGGATGATCGGCGAAGATTTCACGCCCTTCCTGATACTCCAGCTCAATAATGGCGGTGAGGTAGTATTCGCCCTTAGCATTATTGTTGGTCTGCTTc                                                                                                                              <  2:580274/242‑1 (MQ=255)
                                                                                                                      gAGTTAGCGTACCTCGCAGATCAAAACGCGCTGGAGCTGGCAGCAGAACGCCTGCTAACAGCAGTTTGTCATCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGATGCAGGCGGTTAATAACGCCTTCTACTTCGCTTAAACGTTGCGGATGAATGGCGACGATTTCACGCCCTTGCTGATACGCCTGCGCAATAATGTCGGTGATGTAGCATTCGCCCTATGCATTAGTGTCGGTCAGCATCgcc                                                                                                                           <  2:35760/251‑1 (MQ=255)
                                                                                                                                  cccGCAGATCAAAACGCGCTGGATCTCGCAGCATAACGCCTTCTAACAGCAGGTGTGCAGCCTGTTCGGACTGATAAACACTCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCGTAAACGTTGCGGATGAACGGCGACTATTTCACGCCCTTCCTGATACGCCAGCGCAATAATTTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTc                                                                                                               <  2:549777/249‑1 (MQ=255)
                                                                                                                                    cGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCat                                                                                                             <  2:174590/251‑1 (MQ=255)
                                                                                                                                      cAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCatat                                                                                                           <  1:334070/251‑1 (MQ=255)
                                                                                                                                      cAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCatat                                                                                                           <  2:505039/251‑1 (MQ=255)
                                                                                                                                                   gcTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTTATAAACACGCTCCAGACGTGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAAGAATGTCGGTGATGTAGTATTCGCCCTGATCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGc                                                                                              <  1:28601/251‑1 (MQ=255)
                                                                                                                                                   gcTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGc                                                                                              <  1:419366/251‑1 (MQ=255)
                                                                                                                                                        aTCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCCCTTAAACGGTGCGGGTGAACGGCGACGATTTTACGCCCTTCCTGATAAGCCAGCGCAATAATGTCGGGGATGGAGTATTTGCCCCGAGCATTTTTTTTTGTCAGCTTCGCCCGCCAGCGGTTTATATCTTCGcccgttgccatcc                                                                                         >  1:258268/1‑240 (MQ=255)
                                                                                                                                                          tgcctatcaTAACGCCTGCTACCAGCAGTTTTGCAGCCTGTTCGGACTTATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTACGTTTAAACGTTGCGGATGAACGGCGACGATTTCGCGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCGGAGCATTATTGTTGGTCAGCTTCGCCCTCCAGCGTTTCATATCTGCGCCGTTGGCAATCAg                                                                                        <  2:366165/243‑1 (MQ=255)
                                                                                                                                                             gcagcaTAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTTAAAGTTGCGGATGGACGGCGACGATTTCACGCCCTTTCTGATACGCCAGCGCAATAATGTCGGGGATGTAGTATTCGCCCCGAGCCTTATTGTTGGTGAGCTTTGCCCGCCAGCGTTTCATATTTGCGCCGGTTGCAATCAGAATg                                                                                    >  1:481641/1‑251 (MQ=255)
                                                                                                                                                               agaaTAACGCCTGCTAACAGCAGTTTGTCAGCCTGTTCGGACTGATAAACCCGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACGTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCATCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGGCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGcc                                                                                  <  2:354079/248‑1 (MQ=255)
                                                                                                                                                                   tAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCaa                                                                                                                                                                                 <  1:144326/153‑1 (MQ=255)
                                                                                                                                                                    aaCGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGgtg                                                                              <  2:407355/250‑1 (MQ=255)
                                                                                                                                                                     aCGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGtt                                                                            <  1:134260/251‑1 (MQ=255)
                                                                                                                                                                      cGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTg                                                                           <  2:449952/251‑1 (MQ=255)
                                                                                                                                                                       gCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTTACGCCCTCTTCTTCCCTTAAAAGTTGAGGGTGAACGGCGACGATTTCACGCCCCTCCTGGTACGCCCGCGCCATAATGTAGGTGTGGTAGGTTTCGCCCCGAGCATTATTGTTGGTTCGCTTCGCCCGCCCGCGGTTTATTTCTGCGCCGGTGGCAATCAGAAtacaggggttgt                                                                          >  1:249280/1‑240 (MQ=255)
                                                                                                                                                                         cTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCCTCTACTTCGCTTGAACAGTGCGGATGAAAGGCGACGATTTCACGCCCCTACTGGTCCCGCCGCGCAATAATGACGGTGGTGTAAAATTCGCCCCGGGCATTATTGTTTGTCAACTTCGCCACACCGCGATTCATATTTACGCCGCTGGCAATCAAAATGGCGGGGTTGAtc                                                                        >  1:551462/1‑251 (MQ=255)
                                                                                                                                                                         cTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCCTCTACTTCCCTTAAACGTTGCGGATGAAAGGCGAAGATTTTAAGCCCTTTCTGATACACCCACGCAATAAAGTCGGTAATGTAGGATTCGCCCTGGGCCTTATTGTTTGGCAGCTATGCCAGCCAGCGGTTCATATATTCCCCGGTGGGAATTAGAATCCCGGTGTTgttc                                                                        >  1:353064/1‑248 (MQ=255)
                                                                                                                                                                                 agcagTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTGGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAAt                                                                <  2:58893/251‑1 (MQ=255)
                                                                                                                                                                                       ttttCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGAc                                                           <  1:198075/250‑1 (MQ=255)
                                                                                                                                                                                            aGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTTACCCCCTCTACCTTGCCTTAAACTTTCAGAAGGACGGGGACGATTTAACCCCCCTTCCGGTACGCCAGCGCAATAAATTCGGTGATGGAGGATTAGCCCAGCGCAATATTTTTTGTCAACTTCTCCCGCCCGCGTTTCCAATCTGCGCCGTTTGCAAACAGAAAGGCGGGGTTGGTCTCCttcattatagccctcc                                                     >  1:542287/1‑236 (MQ=255)
                                                                                                                                                                                             gCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTTACGCCTTCTACTTCGCTTTAACGTTGCGGGTGAAAGGCGACGGTTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGGGGTGTAGTATTTGCCCCGAGCCTTATTTTTTGTCAGCTTCGCCCGCCAGCGGTTCTTATCTTCGCCGTTTGCGAACAAAATGCCCGGGTTGATCTCCTGAATTTGACGCTGCtc                                                    >  1:553410/1‑251 (MQ=255)
                                                                                                                                                                                             gCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTTCTTCGCTTTAAAGTTGCGGATGGACGGCGACGATTTCACGCCCTTCCTGATACCCCCGCGCCATAATGTCGGGGATGTAGTATTCGCCCTTAGCATTATTTTTTGTCAACTTCGCCCGCCAGCGGTTCTTATCTGCGCCGTTGGCAATTAGAATACCGGGGTTGGTCTCCTTGATCTGGCGCTGCtc                                                    >  1:38341/1‑251 (MQ=255)
                                                                                                                                                                                                   tCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCgg                                               <  2:52117/250‑1 (MQ=255)
                                                                                                                                                                                                                     cccAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTcaac                            <  2:129923/250‑1 (MQ=255)
                                                                                                                                                                                                                      ccAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTcaac                            <  1:368063/250‑1 (MQ=255)
                                                                                                                                                                                                                         gACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATg                        <  1:39909/251‑1 (MQ=255)
                                                                                                                                                                                                                                    tGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACtt              <  1:121672/250‑1 (MQ=255)
                                                                                                                                                                                                                                     gCAGGCGGTTATTCACGCCTTCTACTTCGCTTAACCGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCGTTGTGCTCATCAATGCCGGTAACTTTg            <  2:376287/251‑1 (MQ=255)
                                                                                                                                                                                                                                       aGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGcc          >  1:82734/1‑251 (MQ=255)
                                                                                                                                                                                                                                         gCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGt        <  1:548600/251‑1 (MQ=255)
                                                                                                                                                                                                                                           ggTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGttt      >  2:176540/1‑251 (MQ=255)
                                                                                                                                                                                                                                            tgtATTCACGCCTTCTACTTCGCTTAAACGTTGCGGAGGAACGGCGACGATTTCACTCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGttt      <  2:474603/248‑1 (MQ=255)
                                                                                                                                                                                                                                               aTTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCAc  >  1:87940/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                            
CAATCATCGCCAATCACGCTGTTTTTAATCACGCAACCGGTGCCAATTTTCACGCGATGACCGAGAGTCACGTTGCCCTCGATGATAACGTTAGTATCAATTTCAACATCGCGCCCGTGAGTTAGCGTACCACGCAGATCAAAACGCGCTGGATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCAC  >  NC_000913/3914281‑3914770

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: