Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,604,2370CT23.7% 123.7 / 11.2 59S64L (TCG→TTG) lsrDAutoinducer‑2 ABC transporter membrane subunit LsrD
Reads supporting (aligned to +/- strand):  ref base C (6/39);  new base T (13/1);  total (19/40)
Fisher's exact test for biased strand distribution p-value = 8.25e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACGCCACCGCCATCCGTTAAACCCGCTTCGTCAGGTAAAAAACGGGAGGCCGCATAATGCGTATTCGCTACGGTTGGGAACTGGCTCTTGCCGCACTGCTCGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGATCG‑ATATTTC‑GTTT‑GGTTCGA‑CCATCGGCCTCTGCGCCATTGCATTGGGCGTACTG‑TTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTAACCCGCTGGTGATTACGCTTGGCACGCTGTATCTGTTTGCCGGAAGCGCTCTGCTGCTTTCCGGTATGGCCGGAGCGACGGGGTACGAAGGTATTGGTGGA  >  NC_000913/1603991‑1604487
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                              
aCGCCACCGCCATCCGTTAAACCCGTTTCGTCAGGTAAAAAACGGGAGGCCGCATAATGCGTATTCGCTACGGTTGGGAACTGGCTCTTGCCGCACTGCTCGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGATCG‑ATATTTC‑Gttt                                                                                                                                                                                                                                                           <  2:237175/251‑1 (MQ=255)
  gCCACCGCCATCCGTTAAACCCGCTTCGTCAGGTAAAAAACGGGAGGCCGCATAATGCGTATTCGCTACGGTTGGGAACTGGCTCTTGCCGCACTGCTCGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGGTCG‑GTTTTTT‑GTTT‑gg                                                                                                                                                                                                                                                        >  1:364198/1‑251 (MQ=255)
       cGCCATCCGTTAAACCCGCTTCGTCAGGTAAAAAACGGGAGGCCGCATAATGCGTATTCGCTACGGTTGGGAACTGGCTCTTGCCGCACTGCTCGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGATCG‑ATATTTC‑GTTT‑GGTTCGa                                                                                                                                                                                                                                                   <  2:153978/251‑1 (MQ=255)
       cGCCATCCGTTAAACCCGCTTCGTCAGGTAAAAAACGGGAGGCCGCATAATGCGTATTCGCTACGGTTGGGAACTGGCTCTTGCCGCACTGCTCGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGATCG‑ATATTTC‑GTTT‑GGTTCGa                                                                                                                                                                                                                                                   <  1:230681/251‑1 (MQ=255)
             ccGTTAAACCCGCTTCGTCAGGTAAAAAACGGGAGGCCGCATAATGCGTATTCGCTACGGTTGGGAACTGGCTCTTGCCGCACTGCTCGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGATCG‑ATATTTC‑GTTT‑GGTTCGA‑CCATc                                                                                                                                                                                                                                             <  1:151624/250‑1 (MQ=255)
               gTTAAACCCGCTTCGTCAGGTAAAAAACGGGAGGCCGCATAATGCGTATTCGCTACGGTTGGGAACTGGCTCTTGCCGCACTGCTCGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGATCG‑ATATTTC‑TTTT‑gtgtctcacatcagc                                                                                                                                                                                                                                           >  1:146742/1‑237 (MQ=255)
                 tAAACCCGCTTCGTCAGGTAAAAAACGGGAGGCCGCATAATGCGTATTCGCTACGGTTGGGAACTGGCTCTTGCCGCACTGCTCGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGATCG‑ATATTTC‑GTTT‑GGTTCGA‑CCATCGGcc                                                                                                                                                                                                                                         <  1:370301/250‑1 (MQ=255)
                  aaaCCCGCTTCGTCAGGTAAAAAACGGGAGGCCGCATAATGCGTATTCGCTACGGTTGGGAACTGGCTCTTGCCGCACTGCTCGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGGTCG‑ATATTTC‑GTTT‑GGTTTGA‑CCATCGGCctc                                                                                                                                                                                                                                       >  1:145556/1‑251 (MQ=255)
                   aaCCCGCTTCGTCAGGTAAAAAACGGGAGGCCGCATAATGCGTATTCGCTACGGTTGGGAACTGGCTCTTGCCGCACTGCTCGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGATCG‑ATATTTC‑GTTT‑GGTTCGA‑CCAGCGGCctct                                                                                                                                                                                                                                      <  1:237004/251‑1 (MQ=255)
                            cGTCAGGTAAAAAACGGGAGGCCGCATAATGCGTATTCGCTTCGGTTGGGAACTGGCTCTTGCCGCACTGCTCGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGATCG‑ATATTTC‑GTTT‑GGTTCGA‑CCATCGGCCTCTGCGCcattg                                                                                                                                                                                                                             <  1:404690/251‑1 (MQ=255)
                             gTCAGGTAAAAAACGGGAGGCCGCATAATGCGTATTCGCTACGGTTGGGAACTGGCTCTTGCCGCACTGCTCGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGGACG‑GTATTTT‑GGTT‑TGTTCGG‑CCATCCGCCTCTCCGccctttc                                                                                                                                                                                                                            >  2:53173/1‑246 (MQ=255)
                             gTCAGGTAAAAAACGGGAGGCCGCATAATGCGTATTCGCTACGGTTGGGAACTGGCTCTTGCCGCACTGCTCGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGATCG‑ATATTTC‑GTTT‑GGTTCGA‑CCATCGGCCTCTGCGCcattgc                                                                                                                                                                                                                            <  1:334772/251‑1 (MQ=255)
                                 ggTAAAAAACGGGAGGCCGCATAATGCGTATGCGCTACGGTTGGGAACTGGCTCTTGCCGCACTGCTCGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGATCG‑ATATTTC‑GTTT‑GGTTCGA‑CCATCGGCCTCTGCGCcattgcattg                                                                                                                                                                                                                        <  2:211961/251‑1 (MQ=255)
                                  gTAAAAAACGGGAGGCCGCATAATGCGTATTCGCTACGGTTGGGAACTGGCTCTTGCCGCACTGCTCGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGGTG‑GTGATTGTCAGTGGCGGGATCG‑ATATTTT‑GTTT‑GGGTCGA‑CCATCGGCCTCTGCGCcattgcatttg                                                                                                                                                                                                                       >  1:431830/1‑249 (MQ=255)
                                  gTAAAAAACGGGAGGCCGCATAATGCGTATTCGCTACGGTTGGGAACTGGCTCTTGCCGCACTGCTCGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGATCG‑ATATTTC‑GTTT‑GGTTCGA‑CCATCGGCCTCTGCGCCATTGCATTgg                                                                                                                                                                                                                       <  1:502205/251‑1 (MQ=255)
                                    aaaaaaCGGGAGGCCGCATAATGCGTATTCGCTACGGTTGGGAACTGGCTCTTGCCGCACTGCTCGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGATCG‑ATATTTC‑GTTT‑GGTTCGA‑CCATCGGCCTCTGCGCCATTGCATTGGGc                                                                                                                                                                                                                     <  1:49715/251‑1 (MQ=255)
                                                   gCATAATGCGTATTCGCTACGGTTGGGAACTGGCTCTTGCCGCACTGCTCGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCTGGACCG‑ATATTTC‑GTTT‑GGTTCGA‑CCATCGGCCTCTGCGCCATTGCATTGGGCGTACTG‑TTTCAAAgt                                                                                                                                                                                                     <  1:432466/251‑1 (MQ=255)
                                                    cATAATGCGTATTCGGTACGTTGGGGATATGGCTCTTGCCGCTCTGCTCTTTATTGAGATTGTCGCATTTGGTGCACTTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGATCG‑ATATTTC‑GTTT‑GGTTCGA‑CCATCGGCCTCTGCGCCATTGCATTGGGCGTACTG‑TTGCAAAg                                                                                                                                                                                                      <  2:63503/249‑1 (MQ=255)
                                                       aaTGCGTATTCGCTACGGTTGGGAACTGGCTCTTGCCGCACTGCTCGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGGTCC‑GTCTTTT‑GTTT‑TGTTCGT‑CCATCCGCCCCTCCGCCCTTGCCTTTGGCGATCTG‑TTTCTAAggggtt                                                                                                                                                                                                 >  1:487344/1‑250 (MQ=255)
                                                        aggcttATTCGCTACGGTTGGGAACGGGCTCTTGCCGCACTGCTCGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCTAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGAGTGTCAGTGGCGGTATCG‑ATATTTC‑GTTT‑GGTTCGA‑CCATCGGCCTCTGCGCCATTGCATTGGGCGTACTG‑TTTCAAAGTGgtgt                                                                                                                                                                                                <  1:430745/246‑1 (MQ=255)
                                                         tGCGTATTCGCTACGGTTGGGAACTGGCTCTTGCCGCACTGCTCGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGATAG‑ATTTTTT‑GTTT‑GGTTTGG‑ACACCGGCCTCTGCGGCATTGCCTTgggggtgcctttttaaagggggggg                                                                                                                                                                                                >  1:303896/1‑229 (MQ=255)
                                                                 cGCTTCGGTTGGGAACTTGCGTTTGCCGCACTGCTCGTTATTGAGATTGTCGCATTTTGTGCAATTAACCCGCGAATGTTAGATCGCAAGATGTTGCTGTTAAGCACCTGTGACTTTATCTGCATTGGTATTGTAGCCCGGCCGCGAACG‑TGTGTGATTGTCAGCGGCGCGATCG‑AAGTTGT‑GTTT‑GGTGAAA‑CCATCGGCCTCTTCGCCATTGCAGTGGGCGTATTG‑TTTAAAAGTGGTTGGCCGAttcc                                                                                                                                                                                       <  2:134586/251‑4 (MQ=255)
                                                                      cGGTTGGGAACTGGCTCTTGCCGCACTGCTCGTTAGTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGTCTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGATCG‑ATATTTC‑GTTT‑GGTTCGA‑CCATCGGCCTCTGCGCCATTGCATTGGGCGTACTG‑TTTCAAAGTGGTGTGCCGATGCCGCTgg                                                                                                                                                                                  <  2:119078/251‑1 (MQ=255)
                                                                           gggAACTGGCTCTTGCCGCACTGCTCGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGG‑TAGCATATTTT‑GGTT‑GGTTCGAGCCTTCGGC‑TTTGCGCCATTGCACTGGGTGTTCTG‑TTTCAAAGTGGTGTGTCGGTGcccccggggga                                                                                                                                                                              >  1:133312/1‑241 (MQ=255)
                                                                             gAACTGGCTCTTGCCGCACTGCTCGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGGTCG‑ATATTTT‑GTTT‑TGTTTGG‑CCATCCGCCCCTGCGCCCTTGCATTTGGCGGACTG‑TTTCTAAGTGGGGTGGCGATGCCGCCGGCggtaat                                                                                                                                                                           >  1:422332/1‑246 (MQ=255)
                                                                                  ggCTCTTGCCGCACTGCTCGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGATCG‑ATATTTC‑GTTT‑GGTTCGA‑CCATCGGCCTCTGCGCCATTGCATTGGGCGTACTG‑TTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCt                                                                                                                                                                        <  1:454403/249‑1 (MQ=255)
                                                                                                  cTCGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGATCG‑ATATTTT‑GTTT‑TGTTTGG‑CCATCGGCCTCTGCGCCCTTTCATTTGGGGTGCTG‑TTTCAAAGTGGTGTGCCGATGCCGCCGGCGGTACTCCTGCCCtttatgctccgggc                                                                                                                                                      >  1:343375/1‑239 (MQ=255)
                                                                                                    cGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGATCG‑AAATTTC‑GTTT‑TGGTCGA‑CCACCGGCCTCTGCGCCCTTGCATTGGGGGGACTG‑TTTTAAAGTGGGGTGCCGAGGACGCTGGGGAGACTCCTGACCTTACTGCTCGgcgcct                                                                                                                                                    >  1:236905/1‑249 (MQ=255)
                                                                                                           gagaTTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATCGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGAGCG‑ATATTTC‑GTTT‑GGTTCGA‑CCATCGGCCTCTGCGCCATTGCATTGGGCGTACTG‑TTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCgg                                                                                                                                             <  1:482945/251‑1 (MQ=255)
                                                                                                                 gTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGATCG‑ATATTTC‑GTTT‑GGTTCGA‑CCATCGGCCTCTGCGCCATTGCATTGGGCGTACTG‑TTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGAt                                                                                                                                       <  2:145556/251‑1 (MQ=255)
                                                                                                                   cGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGATCG‑ATATTTC‑GTTT‑GGTTCGA‑CCATCGGCCTCTGCGCCATTGCATTGGGCGTACTG‑TTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATc                                                                                                                                      <  1:322345/250‑1 (MQ=255)
                                                                                                                     cATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGATCG‑ATATTTC‑GTTT‑GGTTCGA‑CCATCGGCCTCTGCGCCATTGCATTGGGCGTACTG‑TTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAAc                                                                                                                                   <  1:212578/251‑1 (MQ=255)
                                                                                                                         tGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGATCG‑ATATTTC‑GTTT‑GGTTCGA‑CCATCGGCCTCTGCGCCATTGCATTGGGCGTACTG‑TTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGcc                                                                                                                                <  1:234022/250‑1 (MQ=255)
                                                                                                                          ggTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGGTTG‑ATAGTTT‑TTTT‑GGTTTGA‑CCATCGGCCTCTGCGCCCCTGCATTGGGGGTGGTG‑TTTCAAAGAAGTGGGCCGATGCCGCCGGCGACACTCCCGACCCTACTGCTCGGCGCACTGTGCGGGCTGGTCAACCCCgg                                                                                                                              >  2:73290/1‑251 (MQ=255)
                                                                                                                           gTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGGTCG‑ATATTTT‑GTTT‑TGTTTGT‑CCAACGGCCTCTTCGCCCTTTCATTTGGCGGCCCG‑TTTCCAAGGGGGGTGCCGGTGCCGCTGGCGGTACCTCTTACCTTACTGCTCCGCGGCTTTTGTGGGCTGGTCAACGcccgg                                                                                                                             >  1:228613/1‑248 (MQ=255)
                                                                                                                              cAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGGTCG‑ATATTTTTGTTT‑G‑TTCGG‑CCATCCGCCCCTGCGCCATTTCATTTGGCGGACTG‑TTTTAAAGTGGGGTGGCGGTGCCGCTGGCGATACTCCTGTCCCTTCTGCTCGGCGCCTTTTGCGGGGTGATCAACGCCgggtt                                                                                                                           >  1:336202/1‑247 (MQ=255)
                                                                                                                               aaTTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGATCG‑GTATTTT‑GTTT‑TGTTTGA‑CCATCCGCCTCTGCGCCCTTTGATTGTGGGTACTG‑TTTTAAAGTGGGGTGTCCATGCCGCCGGGGGTAATCCCTACCCTACTGCTCGGCGGCTTGTTCGGGCGGATCAACGCCGGgtta                                                                                                                          >  1:201114/1‑250 (MQ=255)
                                                                                                                                   aaCCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGATCG‑ATATTTC‑GTTT‑GGTTCGA‑CCATCGGCCTCTGCGCCATTGCATTGGGCGTACTG‑TTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGattaatta                                                                                                                      <  1:53173/250‑1 (MQ=255)
                                                                                                                                    aCCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGGCCG‑ATATTTT‑GTTTTGGTT‑GA‑CCCTCGGGCTCTGGCCCCCTGCTTTGGGGGTACTGGTTTCAAAG‑GGGGTGCCGACGCCGCTGGCGGTACTCCTTACCCTACTGCTCGGGGCACTGTGCGGGGTGGTCAACGACCGATTAATTATc                                                                                                                    >  1:166987/1‑251 (MQ=255)
                                                                                                                                     cccGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGATCG‑ATATTTC‑GTTT‑GGTTCGA‑CCATCGGCCTCTGCGCCATTGCATTGGGCGTACTG‑TTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCt                                                                                                                   <  1:496385/251‑1 (MQ=255)
                                                                                                                                                        aaTATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGTGATCG‑ATATTTC‑GTTT‑GGTTCGA‑CCATCGGCCTCTGCGCCATTGCATTGGGCGTACTG‑TTTCAAAGTGGTGTGCCGATGCCGCTTGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTAACCCGCt                                                                                                <  1:338537/251‑1 (MQ=255)
                                                                                                                                                         atatGTTGCTGTGCAGCACCATTGACTTTAGCTGCATTGTCATTGGCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGATCG‑ATATTTC‑GTTT‑GGTTCGA‑CCATCGGCCTCTGCGCCATTGCATTGGGCGTACTG‑TTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGTCGCATTTTGCGGGCTGATCAACGCCGGATTAATTATCTATCCCAAAGTTAACCCGCTg                                                                                               <  2:10879/251‑1 (MQ=255)
                                                                                                                                                                      cAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATGGTCAGTGGCGGGATCG‑ATATTTC‑GTTT‑GGTTCGA‑CCATCGGCCTCTGCGCCATTGCATTGGGCGTACTG‑TTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGTATTAATTATCTATACCAAAGTTAACCCGCTGGTGATTACGCtt                                                                                   <  1:344112/250‑1 (MQ=255)
                                                                                                                                                                                 cctgtttCTGCAATGGCATTGTCGCCCTACCTCTAACGATG‑GTGATTGTGAGGGGCGGGATAG‑AGATTTC‑GTTT‑GGTTCGA‑CCATCGGCCCGTGCGCCATTGCATTGGGCGTACTG‑TTTCAATGTGGTGTGACGATGCAGCGGGCGATACTCCTGACCTTACGGCTCGGCGCATTGTGCGGGATGATCAACGCCGGATTAATTATTTAGACCAAAGTTAACCCGCTGGTGATTACGCTTTGCACGCTggat                                                                       <  2:71092/245‑4 (MQ=255)
                                                                                                                                                                                 aCTTTAGCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGATCG‑ATATTTC‑GTTT‑GGTTCGA‑CCATCGGCCTCTGCGCCATTGCATTGGGCGTACTG‑TGTCAAAGTGGTGGGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTAACCCGCTGGTGATTACGCTTGGCACGCTGTAt                                                                       <  2:59773/251‑1 (MQ=255)
                                                                                                                                                                                 aCTTTACCTGCATTGGCATTGTCTCCCTTCCGCTAACGGTG‑GTGATTGTGAGTGGCGGGATCG‑ATATTTC‑GTTT‑GGTTCGG‑CCATCGGCCTCTGCGCCATTGCATTGGGCGTACTG‑TTTCAAAGTGGTGTGCCGTGGCCGCTGGCGATACTCCTGTCCATACTGCTCGGCGCATTGTGCGGGCTGATCAGCGCCGGATTAATTATCTATACCAAAGTTATCCCGCTGGTGATTACGCTTGGCACTCTGTAt                                                                       <  2:133312/251‑1 (MQ=255)
                                                                                                                                                                                            attggcattgTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGATCG‑ATATTTC‑GTTT‑GGTTCGA‑CCATCGGCCTCTGCGCCATTGCATTGGGCGTACTG‑TTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTACCCCGCTGGTGATTACGCTTGGCACGCTGTATCTGTTTGCCgg                                                            <  1:151260/251‑1 (MQ=255)
                                                                                                                                                                                                     gTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGAGCG‑ATATTTC‑GTTT‑GGTTCGA‑CCATCGGCCTCTGCGCCATTGCATTGGGCGTACTG‑TTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTAACCCGCTGGTGATTACGCTTGGCACGCTGTATCTGTTTGCCGGAAGCGctct                                                   <  1:207809/251‑1 (MQ=255)
                                                                                                                                                                                                             aCCGCTAACGATG‑GTGATTGTCAGTGGCGGGGTCG‑ATATTTC‑GTTT‑GGTTCGA‑CCATCGGCCTCTGCGCCATTGCAGTGGGCGTACTG‑TTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAATTTAACCCGCTGGTGATTACGCTTGGCACGCTGTATCTGTTTGCCGGAAGCGCTCTGCTGCtt                                            <  2:343375/250‑1 (MQ=255)
                                                                                                                                                                                                                   aaCGGTG‑GTGATTGTCAGTGGCGGGATCG‑ATATTAC‑GTTT‑GGTGCGG‑CCATCGGCCTCTGCGCCATTGCATTGGGCGTCCTG‑TTTCAAAGTGGTGTGACGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTATCCCGCTGGTGATTACGCTTGGCACGCTGTATCTGTTTGCCGGAAGCGCTCTGCGGCTTTCCGGt                                      <  2:32362/250‑1 (MQ=255)
                                                                                                                                                                                                                         g‑gTGATTGTCAGTGGCGGGATCG‑ATATTTC‑GTTT‑GGTTCGA‑CCATCGGCCTCTGCGCCATTGCATTGGGCGTACTG‑TTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTAACCCGCTGGTGATTACGCTTGGCACGCTGTATCTGTTTGCCGGAAGCGCTCTGCTGCTTTCCGGTATGGcc                                <  1:216737/250‑1 (MQ=255)
                                                                                                                                                                                                                          tgtggtagTCAGTGTTGGGATCT‑ATATTTC‑GTTT‑GTTTTGC‑CCATATGTCTCCTCTCCATTGCATTGTGCGTGATG‑TTTCAAAGTTGTGTGCCGTTGCCGCTGGCGACACTCCTGACGTTACTGCCCGTCGCATTGTGCGTGCTTTCCATCGCCGGAGTAAGTATCTATACCAAAGTTAACCCGCTGGTGATTACGCGTGGCACGCGGTATCTGTTTGCCGGAAGCGCTCTGCTGCGTTCCGGTATGGcc                                <  2:348964/243‑1 (MQ=255)
                                                                                                                                                                                                                                tgtcaggtgCGGGATCG‑ATAGTTC‑GTTT‑GGTTCGA‑CCGTCGGCCGCTGCGCCATTGCATTGGGCGTACTG‑TTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGGCCTTACTTCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTAACCCGCTGGTGATTACGCTTGGCACGCTGTATCTGTTTGCCGGAAGCGCTCTGCTGCTTTCCGGTATGGCCGGAGCg                          <  2:228613/242‑1 (MQ=255)
                                                                                                                                                                                                                                 gTCAGTGGCGGGATCG‑ATATTTC‑GTTT‑GGTTCGA‑CCATCGGCCTCTGCGCCATTGCATTGGGCGTACTG‑TTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTAATCAACGCCGGATTAATTATCTATACCAAAGTTAACCCGCTGGTGATTACGCTTGGCACGCTGTATCTGTTTGCCGGAAGCGCTCTGCTGCTTTCCGGTATGGCCGGAGCGa                         <  1:282160/250‑1 (MQ=255)
                                                                                                                                                                                                                                       ggCGGGATCG‑ATATTTC‑GTTT‑GGTTCGA‑CCATCGGCCTCTGCGCCATTGCATTGGGCGTACTG‑TTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTAACCCGCTGGTGATTACGCTTGGCACGCTGTATCTGTTTGCCGGAAGCGCTCTGCTGCTTTCCGGTATGGCCGGAGCGACGGGGTa                  <  1:325981/251‑1 (MQ=255)
                                                                                                                                                                                                                                        gcgtgaccG‑ATATTTC‑GTTT‑GGTTCGA‑CCATCGGCCTCTGCGCCATTGCAGTGGGCGTACTG‑TTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTAACCCGCTGGTGATTACGCTTGGCACGCTGTATCTGTTTGCCGGAAGCGCGCTGCTGCGTTCCGGTATGGCCGGAGCGACGTGGTAc                 <  2:420462/244‑1 (MQ=255)
                                                                                                                                                                                                                                                   tatTTC‑GTTT‑GGTTCGA‑CCATCGGCCTCTGCGCCATTGCATTGGGCGTACTG‑TTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTAACCCGCTGGTGATTACGCTTGGCACGCTGTATCTGTTTGCCGGAAGCGCTCTGCTGCTTTCCGGTATGGCCGGAGCGACGGGGTACGAAGGTATtg       >  1:37275/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    atTTC‑GTTT‑GGTTCGA‑CCATCGGCCTCTGCGCCATTGCATTGGGCGTACTG‑TTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTAACCCGCTGGTGATTACGCTTGGCACGCTGTATCTGTTTGCCGGAAGCGCTCTGCTGCTTTCCGGTATGGCCGGAGCGACGGGGTACGAAGGTATtgg      >  2:123581/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      ttC‑GTTT‑GGTTCGA‑CCATCGGCCTCTGCGCCATTGCATTGGGCGTACTG‑TTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTAACCCGCTGGTGATTACGCTTGGCACGCTGTATCTGTTTGCCGGAAGCGCTCTGCTGCTTTCCGGTATGGCCGGAGCGACGGGGTACGAAGGTATtggtg    >  1:145047/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        c‑GTTT‑GGTTCGA‑CCATCGGCCTCTGCGCCATTGCATTGGGCGTACTG‑TTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTAACCCGCTGGTGATTACGCTTGGCACGCTGTATCTGTTTGCCGGAAGCGCTCTGCTGCTTTCCGGTATGGCCGGAGCGACGGGGTACGAAGGTATTGGTGGa  <  1:180991/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                              
ACGCCACCGCCATCCGTTAAACCCGCTTCGTCAGGTAAAAAACGGGAGGCCGCATAATGCGTATTCGCTACGGTTGGGAACTGGCTCTTGCCGCACTGCTCGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATG‑GTGATTGTCAGTGGCGGGATCG‑ATATTTC‑GTTT‑GGTTCGA‑CCATCGGCCTCTGCGCCATTGCATTGGGCGTACTG‑TTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTAACCCGCTGGTGATTACGCTTGGCACGCTGTATCTGTTTGCCGGAAGCGCTCTGCTGCTTTCCGGTATGGCCGGAGCGACGGGGTACGAAGGTATTGGTGGA  >  NC_000913/1603991‑1604487

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: ‑