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.

CATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTG‑TCG‑GATGCGCCGTC‑AGAGTGGCGTATCC‑GATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGC‑TCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATG‑ACGTTTACTTCGCTTTACCCTGGTTTG‑CAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAA  >  NC_000913/786542‑787035
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                       
cATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTG‑TCG‑GATGCGCCGTC‑AGAGTGGCGTATCC‑GATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGC‑TCCATTGTTAGcaaca                                                                                                                                                                                                                                                        <  2:504539/251‑1 (MQ=255)
 aTCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTG‑TCG‑GATGCGCCGTC‑AGAGTGGCGTATCC‑GATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGC‑TCCATTGTTAGcaaca                                                                                                                                                                                                                                                        >  2:505939/1‑250 (MQ=255)
               aaTCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTG‑TCG‑GATGCGCCGTC‑AGAGTGGCGTATCC‑GATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGC‑TCCATTGTTAGCAACAAAAAAGCCGACTCAc                                                                                                                                                                                                                                         >  1:493073/1‑251 (MQ=255)
                aTCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTG‑TCG‑GATGCGCCGTC‑AGAGTGGCGTATCC‑GATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGC‑TCCATTGTTAGCAACAAAAAAGCCGACTCACt                                                                                                                                                                                                                                        >  2:224704/1‑251 (MQ=255)
                        cTCGAGCGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGATGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTG‑TCG‑GATGCGCCGTC‑AGAGTGGCGTATCC‑GATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGC‑TCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCg                                                                                                                                                                                                                                >  2:98448/1‑251 (MQ=255)
                               gCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTG‑TCG‑GATGCGCCGTC‑AGAGTGGCGTATCC‑GATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGC‑TCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTtct                                                                                                                                                                                                                         <  2:477000/251‑1 (MQ=255)
                               gCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTG‑TCG‑GATGCGCCGTC‑AGAGTGGCGTATCC‑GATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGC‑TCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTtct                                                                                                                                                                                                                         <  1:525411/251‑1 (MQ=255)
                                 ggggAGCCGCTGGCCTACGGTAAGAGCATCTCCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGGCGCGGGTAAGGTTTAATTG‑TCG‑GATGCGCCGTC‑AGAGTGGCGTATCC‑GATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAAGGTGC‑TCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGGCGGCTTTCTCa                                                                                                                                                                                                                       <  2:388256/251‑1 (MQ=255)
                                 ggggAGCCGCTGGCCTACGGTAAGAGCATCCCCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTG‑TCG‑GATGCGCCGTC‑AGAGTGGCGTATCC‑GATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGC‑TCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTtctc                                                                                                                                                                                                                        >  1:454207/1‑250 (MQ=255)
                                       ccGCTGGCCTACGGTAAGAGTATCTCCGATGCCTGCATCGGCAGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTG‑TCG‑GATGCGCCGTC‑AGAGTGGCGTATCC‑GATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGC‑TCCATTGTTAGCAACAAAAAAGCCGACTCACCTGCGGTCGGCTTTCTCATTTTaa                                                                                                                                                                                                                 <  2:94275/251‑1 (MQ=255)
                                            ggCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTG‑TCG‑GATGCGCCGTC‑AGAGTGGCGTATCC‑GATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGC‑TCCCTTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACgc                                                                                                                                                                                                             >  2:7949/1‑249 (MQ=255)
                                             gCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTG‑TCG‑GATGCGCCGTC‑AGAGTGGCGTATCC‑GATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGC‑TCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAAt                                                                                                                                                                                                           >  2:276842/1‑251 (MQ=255)
                                                    gTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGATCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTG‑TCG‑GATGCGCCGTC‑AGAGTGGCGTATCC‑GATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGC‑TCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATG‑ACGttt                                                                                                                                                                                                   >  1:540597/1‑251 (MQ=255)
                                                     tAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTG‑TCG‑GATGCGCCGTC‑AGAGTGGCGTATCC‑GATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGC‑TCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATG‑ACGttt                                                                                                                                                                                                   >  1:120363/1‑250 (MQ=255)
                                                       agagCATCACCGATGCCTGCATCGGCTGGGCAGGTACCGATGCTCCGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTG‑TCG‑GATGCGCCGTA‑AGAGTGGCGTATCC‑GATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGC‑TCCATTGTTCGCAACAACAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTTAACGCATA‑ACGTTTACt                                                                                                                                                                                                >  2:412318/1‑251 (MQ=255)
                                                        gagCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTG‑TCG‑GATGCGCCGTC‑AGAGTGGCGTATCC‑GATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGC‑TCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATG‑ACGTTTACtt                                                                                                                                                                                               >  1:162119/1‑251 (MQ=255)
                                                           cATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTG‑TCG‑GATGCGCCGTC‑AGAGTGGCGTATCC‑GATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGC‑TCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAATCGAATG‑ACGTTTACTTCGc                                                                                                                                                                                            <  1:93874/251‑1 (MQ=255)
                                                            aTCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTG‑TCG‑GATGCGCCGTC‑AGAGTGGCGTATCC‑GATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGC‑TCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATG‑ACGTTTACTTCGCt                                                                                                                                                                                           <  1:2141/251‑1 (MQ=255)
                                                                ccGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTG‑TCG‑GATGCGCCGTC‑AGAGTGGCGTATCC‑GATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGC‑TCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATG‑ACGTTTACTTCGCTTTAc                                                                                                                                                                                       <  1:229632/251‑1 (MQ=255)
                                                                         gCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTG‑TCG‑GATGCGCCGTC‑AGAGTGGCGTATCC‑GATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGC‑TCCATTGTTAGCAACAAAAAATCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATG‑ACGTTTACTTCGCTTTACCCTGGttt                                                                                                                                                                               <  1:515543/250‑1 (MQ=255)
                                                                                   ggAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTG‑TCG‑GATGCGCCGTC‑AGAGTGGCGTATCC‑GATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGC‑TCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATG‑ACGTTTACTTCGCTTTACCCTGGTTTG‑CAA‑Ccgccgc                                                                                                                                                                  >  2:30131/1‑251 (MQ=255)
                                                                                            cGATGCTCTGTTACGTCACCTGGCGAATGCAGTACAAGCGCGTCGCGGGTAAGGTTTAATTG‑TCG‑GATGCGCCGTC‑AGAGTGGCGTATCC‑GATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCCGGCAATGTGC‑TCCATTGTTAGCACCAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATG‑ACGTTTACTTCGCTTTACCCTGGTTTG‑CAA‑CCGCCGCTGCTTTCGc                                                                                                                                                         >  2:111847/1‑251 (MQ=255)
                                                                                                     gTTACGTCATCTGGAGAATGCGGAAAAAGCGCGTCGGGGGTAAGGTTTAATTG‑TGG‑GATGGGCCGTC‑AGAGGGGCGTATCC‑GATGAATCCCCACAGGCCTGATAAGCCGCGCAGCGTCGAATCAGGCAATGTGC‑TCCATTGTTAGCAACAAAAAAGCCGCCTCCCTTGCAGTCGGTTTTCTCATTTAAAACGAATG‑ACGTTTCCTTCGCTTTCCCCTGGTTTG‑CAA‑CCGCCGCTGCTTTCGCTGCGAtctc                                                                                                                                                <  1:398340/251‑1 (MQ=255)
                                                                                                                gggCGAAGGCAGAAAAAGCGCGGCG‑GGGAAAAGTTTTATTGGCCG‑GATGCGCGGCC‑AGAGTGGCG‑ATCCCGATAAATCACCACAGCCCTGATAAGTCGCGCGGCGTCGCAACAGGAAATGTCC‑TCCATTGTTAGCAAAAAAAAGGCCGCCGCCCTTGCAGCCGTCTTCCTCTTTTAAAACGAATG‑CCGTTTACTCCGCTTT‑CCCCGGTTTGGCAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAc                                                                                                                                     <  2:120363/250‑1 (MQ=255)
                                                                                                                          aGTAAAAGCGCGTCGCGGGTAAGGTTTAATTG‑TCG‑GATGCGCCGTC‑AGAGTGGCGTATCC‑GATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGC‑TCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATG‑ACGTTTACTTCGCTTTACCCTGGTTTG‑CAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGataat                                                                                                                            >  2:97830/1‑250 (MQ=255)
                                                                                                                             aaaaGCGCGTCGCGGGTAAGGTTTAATTG‑TCG‑GATGCGCCGTC‑AGAGTGGCGTATCC‑GATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGC‑TCCATTGTTAGCAACACAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATG‑ACGTTTACTTCGCTTTACCCTGGTTTG‑CAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCg                                                                                                                        >  1:245273/1‑251 (MQ=255)
                                                                                                                                   gcgTCGGGGGTAAGGTTTAATTG‑‑CGTGATGC‑CCCTCTAGGGTGGCTTTTCC‑AATAAATCCCCACAGGCCTGAAAATGCGCGCAGCGTCGCATCGGGCTTGGTGC‑TCCATTTTTAGAAAAAAAAAAGCCGCCTCACTTGCAGTCGTTTTTCTCATTTAAACCGATGG‑ACGTTTCCTTC‑CCTTGCCCTGGGTTG‑CAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAg                                                                                                                  <  2:585582/251‑1 (MQ=255)
                                                                                                                                                    tAATTG‑TCG‑GATGCGCCGTC‑AGAGTGGCGTATCC‑GATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGC‑TCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATG‑ACGTTTACTTCGCTTTACCCTGGTTTG‑CAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgt                                                                                                 >  1:330489/1‑251 (MQ=255)
                                                                                                                                                                        cc‑AGAGTGCCTTATCC‑GATTAACCCCCACAGGCCTGTAAAGCCGGCCGGCTTCGCATCGGGCATTTTGC‑GCAATTGTAACCAACAAAAAGCCCGCCTCTCTTGCAGTCGGCTTTCTCATTTTAAACTAATG‑ACGTTTCCTTCTCTTTCCCCTGGTTTC‑TAC‑CCGTCGCTGCTGTCGCTGCGATCGCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCg                                                                               <  1:271782/250‑1 (MQ=255)
                                                                                                                                                                             agTGGCGTATCC‑GATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGC‑TCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGCATG‑ACGTTTACTTCGCTTTACCCTGGTTTG‑CAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCAc                                                                           >  1:211107/1‑251 (MQ=255)
                                                                                                                                                                                       cc‑GATGAATCCCCACAGGCTGGATAAGTCGCGCCGCGCCGCATCAGGCAATGTCC‑TCAATTTTTAGAAAAAAAAAAGCCGCCTCACTGGCAGTCGTCTTTCTCATTTAAAACGATTG‑ACGTTTACTTCGCTTTACCCTGGTTTG‑CAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCggg                                                                 <  2:437574/251‑1 (MQ=255)
                                                                                                                                                                                                acaccaaatgCCTGATAAGCCGCGCAGCTTCGCATCAGGCAATTTTC‑TCCATTGTTAGCAAAAAAAAAGCAGACCCATTTGCAGGCGGCTTTCTCATTTTAAACGAATG‑ACGTTTACTCCGTTTTCCCCTGTTTTG‑CAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAg                                                         <  1:505939/242‑1 (MQ=255)
                                                                                                                                                                                                       aGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGC‑TCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATG‑ACGTTTACTTCGCTTTACCCTGGTTTG‑CAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCa                                                     >  1:132788/1‑248 (MQ=255)
                                                                                                                                                                                                         gCCTGAAAAGTCGCGCAGCTTCGCATCAGGCATTGGGC‑TCTATTGTTAAAAAAAAAAACGCCGACTCACTTGAAGTCGGCTTTCTCATTTTAAA‑GAAAGTACGTTTACTTCGTTTTCCCCGGGTTTG‑CAT‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAAt                                                <  1:77958/251‑1 (MQ=255)
                                                                                                                                                                                                               tactccccgcAGCATCTCATCAGGAAATTTGC‑CCCTTTTTTAGCAACAAAAAAGCCCCCTCACTTGCTGTCGGCTTTCTCATTTTAAACAATTG‑ACGTTTCCTTCGCTTTCCCCTGGTTTG‑CAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCGCGCCAGTCGGGATATTAAGCTCAAGAAGCtctt                                           <  1:169428/243‑1 (MQ=255)
                                                                                                                                                                                                                aaGTCGCGAAGCGACGCTCCAGGCAA‑GTGCCACCATTTTAAGCAAAAAAAAAGCCGCCCCCCTGCCAGCCGGCTTCCTCATTTTAAACGATGG‑CCTTTTGCTTCTCTTTACCCTGGTTTG‑CAT‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCGCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCtcttct                                         <  2:528373/251‑1 (MQ=255)
                                                                                                                                                                                                                  gccgcccggCGTCGCATCAGGCAATGTGC‑TCCTTTTTAAGCAAAAAAAAAGCCCCCTCCCTGGCGGTCGTCTTTCTCATTTTAAACGAATG‑CCGTTTCCTCCGCTTTACCCTGTTTTG‑CAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCttcttc                                       <  2:132788/243‑1 (MQ=255)
                                                                                                                                                                                                                  gTCGCGCAGCGCCGATTCAGGAAATGTGC‑TCCATTGTTAGAAAAAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTAAAACGAATG‑ACGTTTACTTCGCTTTCCCCTGGTTTG‑CAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCttcttc                                       <  1:276842/251‑1 (MQ=255)
                                                                                                                                                                                                                      cgcAGCGTCGCATCAGGCAATGTGC‑TCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATG‑ACGTTTACTTCGCTTTACCCTGGTTTG‑CAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTc                                   >  2:218906/1‑251 (MQ=255)
                                                                                                                                                                                                                      cgcAGCGTCGCATCAGGCAATGTGC‑TCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATG‑ACGTTTACTTCGCTTTACCCTGGTTTG‑CAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTc                                   >  1:405669/1‑251 (MQ=255)
                                                                                                                                                                                                                                   cAGGAAATGTCC‑TCCATTTTTAGCAAAAAAAAACCCGTCCCACTTGCAGTCGGCTTTCTCATTTAAACCGAATG‑CCGTTTACTTCGCTTTCCCCTGTTTTG‑CAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGat                      <  1:358791/251‑1 (MQ=255)
                                                                                                                                                                                                                                     ggCAATGTGC‑GCCATTGTTAGCAAAAAAAAAGCCGACTCTCTTGCAGTCGGTTTTCTCATTTTAAACGAATG‑ACTTTTACTTCGCTTTCCCCTGGTTTG‑CAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGatat                    <  1:224704/251‑1 (MQ=255)
                                                                                                                                                                                                                                     aggcaaggggccccATTGTAAGCAAAAAAAAAGCGGACTCACTTGCAGTCGGTTTTCTAATTTAAAACGAAGG‑ACTTTTACTTCGCTTTCCCCTGGTTGG‑CAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGata                     <  1:97830/239‑1 (MQ=255)
                                                                                                                                                                                                                                        aagggcc‑TCCATTGTTAGCAAAAAAAAAGCCGACTCACTTGCAGTCGGCTTTCGCATTTAAAACGAATG‑ACGTTTCCTTCGCTTCCCCCGGGTTTG‑CAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTc                 <  1:98448/245‑1 (MQ=255)
                                                                                                                                                                                                                                                      gTTAGCAAAAAAAAAGCCGACCCCCTTGCAGTCGTCTTTCTCATTTAAAACGAATG‑CCTTTTTCTTCGCTTTACCCTGTTTTG‑CAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTa    <  1:44328/251‑1 (MQ=255)
                                                                                                                                                                                                                                                         aGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATG‑ACGTTTACTTCGCTTTACCCTGGTTTG‑CAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTaa   >  2:342316/1‑249 (MQ=255)
                                                                                                                                                                                                                                                         aGCAAAAAAAAAGCCGACTCACTTGCAGTCGGTTTTCTCATTTTAAACGAATG‑ACGTTTACTTCGCTTTACCCTGGTTTG‑CAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATATACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTaaa  <  1:7949/250‑1 (MQ=255)
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                       
CATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTG‑TCG‑GATGCGCCGTC‑AGAGTGGCGTATCC‑GATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGC‑TCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATG‑ACGTTTACTTCGCTTTACCCTGGTTTG‑CAA‑CCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATGTTATCAAGATATTTCACCAGCGCACGTAAA  >  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: ‑