Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913786,7910CA31.1% 68.6 / 13.9 45intergenic (+106/+52)aroG/gpmA3‑deoxy‑7‑phosphoheptulonate synthase, Phe‑sensitive/2,3‑bisphosphoglycerate‑dependent phosphoglycerate mutase
Reads supporting (aligned to +/- strand):  ref base C (20/11);  new base A (0/14);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 3.86e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAA  >  NC_000913/786542‑787035
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                       
cATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGcaaca                                                                                                                                                                                                                                                        <  2:504539/251‑1 (MQ=255)
 aTCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGcaaca                                                                                                                                                                                                                                                        >  2:505939/1‑250 (MQ=255)
               aaTCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCAc                                                                                                                                                                                                                                         >  1:493073/1‑251 (MQ=255)
                aTCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACt                                                                                                                                                                                                                                        >  2:224704/1‑251 (MQ=255)
                        cTCGAGCGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGATGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCg                                                                                                                                                                                                                                >  2:98448/1‑251 (MQ=255)
                               gCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTtct                                                                                                                                                                                                                         <  2:477000/251‑1 (MQ=255)
                               gCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTtct                                                                                                                                                                                                                         <  1:525411/251‑1 (MQ=255)
                                 ggggAGCCGCTGGCCTACGGTAAGAGCATCTCCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGGCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAAGGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGGCGGCTTTCTCa                                                                                                                                                                                                                       <  2:388256/251‑1 (MQ=255)
                                 ggggAGCCGCTGGCCTACGGTAAGAGCATCCCCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTtctc                                                                                                                                                                                                                        >  1:454207/1‑250 (MQ=255)
                                       ccGCTGGCCTACGGTAAGAGTATCTCCGATGCCTGCATCGGCAGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACCTGCGGTCGGCTTTCTCATTTTaa                                                                                                                                                                                                                 <  2:94275/251‑1 (MQ=255)
                                            ggCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCCTTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACgc                                                                                                                                                                                                             >  2:7949/1‑249 (MQ=255)
                                             gCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAAt                                                                                                                                                                                                           >  2:276842/1‑251 (MQ=255)
                                                    gTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGATCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGttt                                                                                                                                                                                                   >  1:540597/1‑251 (MQ=255)
                                                     tAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGttt                                                                                                                                                                                                   >  1:120363/1‑250 (MQ=255)
                                                       agagCATCACCGATGCCTGCATCGGCTGGGCAGGTACCGATGCTCCGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTAAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTCGCAACAACAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTTAACGCATAACGTTTACt                                                                                                                                                                                                >  2:412318/1‑251 (MQ=255)
                                                        gagCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACtt                                                                                                                                                                                               >  1:162119/1‑251 (MQ=255)
                                                           cATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAATCGAATGACGTTTACTTCGc                                                                                                                                                                                            <  1:93874/251‑1 (MQ=255)
                                                            aTCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCt                                                                                                                                                                                           <  1:2141/251‑1 (MQ=255)
                                                                ccGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTAc                                                                                                                                                                                       <  1:229632/251‑1 (MQ=255)
                                                                         gCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAATCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGttt                                                                                                                                                                               <  1:515543/250‑1 (MQ=255)
                                                                                   ggAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACcgccgc                                                                                                                                                                  >  2:30131/1‑251 (MQ=255)
                                                                                            cGATGCTCTGTTACGTCACCTGGCGAATGCAGTACAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCCGGCAATGTGCTCCATTGTTAGCACCAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGc                                                                                                                                                         >  2:111847/1‑251 (MQ=255)
                                                                                                     gTTACGTCATCTGGAGAATGCGGAAAAAGCGCGTCGGGGGTAAGGTTTAATTGTGGGATGGGCCGTCAGAGGGGCGTATCCGATGAATCCCCACAGGCCTGATAAGCCGCGCAGCGTCGAATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGCCTCCCTTGCAGTCGGTTTTCTCATTTAAAACGAATGACGTTTCCTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGAtctc                                                                                                                                                <  1:398340/251‑1 (MQ=255)
                                                                                                                gggCGAAGGCAGAAAAAGCGCGGCGGGGAAAAGTTTTATTGGCCGGATGCGCGGCCAGAGTGGCGATCCCGATAAATCACCACAGCCCTGATAAGTCGCGCGGCGTCGCAACAGGAAATGTCCTCCATTGTTAGCAAAAAAAAGGCCGCCGCCCTTGCAGCCGTCTTCCTCTTTTAAAACGAATGCCGTTTACTCCGCTTTCCCCGGTTTGGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAc                                                                                                                                     <  2:120363/250‑1 (MQ=255)
                                                                                                                          aGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGataat                                                                                                                            >  2:97830/1‑250 (MQ=255)
                                                                                                                             aaaaGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACACAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCg                                                                                                                        >  1:245273/1‑251 (MQ=255)
                                                                                                                                   gcgTCGGGGGTAAGGTTTAATTGCGTGATGCCCCTCTAGGGTGGCTTTTCCAATAAATCCCCACAGGCCTGAAAATGCGCGCAGCGTCGCATCGGGCTTGGTGCTCCATTTTTAGAAAAAAAAAAGCCGCCTCACTTGCAGTCGTTTTTCTCATTTAAACCGATGGACGTTTCCTTCCCTTGCCCTGGGTTGCAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAg                                                                                                                  <  2:585582/251‑1 (MQ=255)
                                                                                                                                                    tAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgt                                                                                                 >  1:330489/1‑251 (MQ=255)
                                                                                                                                                                        ccAGAGTGCCTTATCCGATTAACCCCCACAGGCCTGTAAAGCCGGCCGGCTTCGCATCGGGCATTTTGCGCAATTGTAACCAACAAAAAGCCCGCCTCTCTTGCAGTCGGCTTTCTCATTTTAAACTAATGACGTTTCCTTCTCTTTCCCCTGGTTTCTACCCGTCGCTGCTGTCGCTGCGATCGCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCg                                                                               <  1:271782/250‑1 (MQ=255)
                                                                                                                                                                             agTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGCATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCAc                                                                           >  1:211107/1‑251 (MQ=255)
                                                                                                                                                                                       ccGATGAATCCCCACAGGCTGGATAAGTCGCGCCGCGCCGCATCAGGCAATGTCCTCAATTTTTAGAAAAAAAAAAGCCGCCTCACTGGCAGTCGTCTTTCTCATTTAAAACGATTGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCggg                                                                 <  2:437574/251‑1 (MQ=255)
                                                                                                                                                                                                acaccaaatgCCTGATAAGCCGCGCAGCTTCGCATCAGGCAATTTTCTCCATTGTTAGCAAAAAAAAAGCAGACCCATTTGCAGGCGGCTTTCTCATTTTAAACGAATGACGTTTACTCCGTTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAg                                                         <  1:505939/242‑1 (MQ=255)
                                                                                                                                                                                                       aGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCa                                                     >  1:132788/1‑248 (MQ=255)
                                                                                                                                                                                                         gCCTGAAAAGTCGCGCAGCTTCGCATCAGGCATTGGGCTCTATTGTTAAAAAAAAAAACGCCGACTCACTTGAAGTCGGCTTTCTCATTTTAAAGAAAGTACGTTTACTTCGTTTTCCCCGGGTTTGCATCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAAt                                                <  1:77958/251‑1 (MQ=255)
                                                                                                                                                                                                               tactccccgcAGCATCTCATCAGGAAATTTGCCCCTTTTTTAGCAACAAAAAAGCCCCCTCACTTGCTGTCGGCTTTCTCATTTTAAACAATTGACGTTTCCTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCGCGCCAGTCGGGATATTAAGCTCAAGAAGCtctt                                           <  1:169428/243‑1 (MQ=255)
                                                                                                                                                                                                                aaGTCGCGAAGCGACGCTCCAGGCAAGTGCCACCATTTTAAGCAAAAAAAAAGCCGCCCCCCTGCCAGCCGGCTTCCTCATTTTAAACGATGGCCTTTTGCTTCTCTTTACCCTGGTTTGCATCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCGCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCtcttct                                         <  2:528373/251‑1 (MQ=255)
                                                                                                                                                                                                                  gccgcccggCGTCGCATCAGGCAATGTGCTCCTTTTTAAGCAAAAAAAAAGCCCCCTCCCTGGCGGTCGTCTTTCTCATTTTAAACGAATGCCGTTTCCTCCGCTTTACCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCttcttc                                       <  2:132788/243‑1 (MQ=255)
                                                                                                                                                                                                                  gTCGCGCAGCGCCGATTCAGGAAATGTGCTCCATTGTTAGAAAAAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTAAAACGAATGACGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCttcttc                                       <  1:276842/251‑1 (MQ=255)
                                                                                                                                                                                                                      cgcAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTc                                   >  2:218906/1‑251 (MQ=255)
                                                                                                                                                                                                                      cgcAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTc                                   >  1:405669/1‑251 (MQ=255)
                                                                                                                                                                                                                                   cAGGAAATGTCCTCCATTTTTAGCAAAAAAAAACCCGTCCCACTTGCAGTCGGCTTTCTCATTTAAACCGAATGCCGTTTACTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGat                      <  1:358791/251‑1 (MQ=255)
                                                                                                                                                                                                                                     ggCAATGTGCGCCATTGTTAGCAAAAAAAAAGCCGACTCTCTTGCAGTCGGTTTTCTCATTTTAAACGAATGACTTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGatat                    <  1:224704/251‑1 (MQ=255)
                                                                                                                                                                                                                                     aggcaaggggccccATTGTAAGCAAAAAAAAAGCGGACTCACTTGCAGTCGGTTTTCTAATTTAAAACGAAGGACTTTTACTTCGCTTTCCCCTGGTTGGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGata                     <  1:97830/239‑1 (MQ=255)
                                                                                                                                                                                                                                        aagggccTCCATTGTTAGCAAAAAAAAAGCCGACTCACTTGCAGTCGGCTTTCGCATTTAAAACGAATGACGTTTCCTTCGCTTCCCCCGGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTc                 <  1:98448/245‑1 (MQ=255)
                                                                                                                                                                                                                                                      gTTAGCAAAAAAAAAGCCGACCCCCTTGCAGTCGTCTTTCTCATTTAAAACGAATGCCTTTTTCTTCGCTTTACCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTa    <  1:44328/251‑1 (MQ=255)
                                                                                                                                                                                                                                                         aGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTaa   >  2:342316/1‑249 (MQ=255)
                                                                                                                                                                                                                                                         aGCAAAAAAAAAGCCGACTCACTTGCAGTCGGTTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATATACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTaaa  <  1:7949/250‑1 (MQ=255)
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                       
CATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAA  >  NC_000913/786542‑787035

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: