Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,258,9960TG30.8% 131.0 / 20.0 78V69G (GTG→GGG) ychHDUF2583 family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (11/43);  new base G (24/0);  total (35/43)
Fisher's exact test for biased strand distribution p-value = 5.28e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCCAGAAAGTTCCGGAGGTGACATATGAAACGCAAAAACGCTTCGTTACTCGGTAACGTGCTCATGGGGTTGGGTCTGGTGGTAATGGTGGTCGGCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAATCTCGTACTTATAAATCCGCCATCGCCGCGCGACGA‑TTCGGGAAGAACGCCAGGCGTCCCGGGATCGGTTGAATGCCAGCGCGAGCCATAGTGCGCAGTGGC  >  NC_000913/1258767‑1259217
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                               
cccAGAAAGTTCCGGAGGTGACATATGAAACGCAAAAACGCTTCGTTACTCGGTAACGTGCTCATGGGGTTGGGTCTGGTGGTAATGGTGGTCGGCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGGGGCCGGGTGTGGGGGGGTTGAAAAGGTGTC‑GGAcc                                                                                                                                                                                                                        >  2:200523/1‑239 (MQ=255)
cccAGAAAGTTCCGGAGGAGACATATGAAACGCAAAAACGCTTCGTTACTCGGTAACGTGCTCATGGGGTTGGGTCTGGCGGTAATGGTGGTCGGCGTGGGGGATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTCCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGGGCCCCGGC‑TGGAGGGCATTGACAGTTGTG‑Cg                                                                                                                                                                                                                           >  1:937584/1‑235 (MQ=255)
 ccAGAAAGTTCCGGAGGTGACATATGAAACGCAAAAACGCTTCGTTACTCGGTAACGTGCTCATGGGGTTGGGTCTGGTGGTAATGGTGGTCGGCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑CGACCGt                                                                                                                                                                                                                      <  1:837523/239‑1 (MQ=255)
    gAAAGTTCCGGAGGTGACATATGAAACGCAAAAACGCTTCGTTACTCGGTAACGTGCTCATGGGGTTGGGTCTGGTGGTAATGGTGGTCGGCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGGGCCCGTGT‑TGGCGGGGGTGTACAGGTGTG‑GGGCCGttt                                                                                                                                                                                                                    >  2:1143041/1‑237 (MQ=255)
           ccGGAGGTGACATATGAAACGCAAAAACGCTTCGTTACTCGGTAACGTGCTCATGGGGTTGGGTCTGGTGGTAATGGTGGTCGGCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTggg                                                                                                                                                                                                            <  1:200523/239‑1 (MQ=255)
                 gTGACATATGAAACGCAAAAACGCTTCGTTACTCGGTAACGTGCTCATGGGGTTGGGTCTGGTGGTAATGGTGGTCGGCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGCGT‑TTGCGGGGATGGACAAGGGGG‑GTGCCGTTACTgggggggttgc                                                                                                                                                                                                       >  2:1285897/1‑233 (MQ=255)
                                         ttCGTTACTCGGTAACGTGCTCAGGGGGTTGGGTCTGGTGGTAATGGTGGTCGGGGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCGGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑GGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGt                                                                                                                                                                               <  1:247802/238‑1 (MQ=255)
                                                     tAACGTGCTCATGGGGTTGGGTCTGGTGGTAATGGTGGTCGGCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGTCATGAACAGGTGTG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCg                                                                                                                                                                    <  1:1401735/237‑1 (MQ=255)
                                                        cGTGCTCATGGGGTTGGGTCTGGTGGTAATGGTGGTCGGCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTGGGTTCGCCACTATGACa                                                                                                                                                                                             >  1:331728/1‑209 (MQ=255)
                                                        cGTGCTCATGGGGTTGGGTCTGGTGGTAATGGTGGTCGGCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTGGGTTCGCCACTATGACa                                                                                                                                                                                             <  2:331728/209‑1 (MQ=255)
                                                               aTGGGGTTGGGTCTGGTGGTAATGGTGGTCGGCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGGGTGAACAGGTGGG‑CGGCCGTTACTGGGGGGGTCGCCACCATGGCAAACGCTGCCGCCGGCGCGGTAAACGCCGTCATAGCt                                                                                                                                                         >  2:1401735/1‑238 (MQ=255)
                                                               aTGGGGTTGGGTCTGGTGGTAATGGTGGTCGGCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTa                                                                                                                                                        <  2:905964/239‑1 (MQ=255)
                                                               aTGGGGTTGGGTCTGGTGGTAATGGTGGTCGGCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTa                                                                                                                                                        <  1:1091921/239‑1 (MQ=255)
                                                               aTGGGGTTGGGTCTGGTGGTAATGGTGGTCGGCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTa                                                                                                                                                        <  1:1047396/239‑1 (MQ=255)
                                                                   ggTTGGGTCTGGTGGTAATGGTGGTCGGCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCTCGTGT‑TGGCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTa                                                                                                                                                        <  1:1058974/235‑1 (MQ=255)
                                                                    gTTGGGTCTGGTGGTAATGGTGGTCGGCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGGCTTAACAGGGGTG‑CGACCGTTACTGGTGGGGTCGCCCCTATGACCAACGGTGCCGGCGTAGCGGTAATTGCCGTCCTAGCTaaaaaa                                                                                                                                                   >  1:752064/1‑239 (MQ=255)
                                                                           tggggggtatgggtggtCGGCGTGGGTTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTTTTTCGCCCATGGTGCAGTGCTAAGTATTTTCGTCGGGGCCATTCTCTGGCTGGCGGGTGCCCGGGT‑TGGCGGGCATGAACAGGTGGG‑CGACCGTTACGGGTGGGTTCGCCACTATGAAAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCa                                                                                                                                            <  2:146074/228‑1 (MQ=255)
                                                                                         ggtCGGCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACcag                                                                                                                               <  1:1285897/238‑1 (MQ=255)
                                                                                                        ttCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGa                                                                                                                 <  1:715942/237‑1 (MQ=255)
                                                                                                         tCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTg                                                                                                              <  1:694077/239‑1 (MQ=255)
                                                                                                           aaTCCTCACCCAGTTACCCAAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGt                                                                                                            <  2:1407980/239‑1 (MQ=255)
                                                                                                           aaTCCTCAACCAGTTACCACATTTTAATATGCCCCAGTTTTTCGCACTTGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGt                                                                                                            <  2:752064/239‑1 (MQ=255)
                                                                                                             tCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGGGTG‑CGACCGTTACTGGGGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAAAAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCa                                                                                                          >  1:567204/1‑239 (MQ=255)
                                                                                                                  aaCCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGGATGAAACGGTGTG‑CGACCGTTACTGGGGGGTTCGCCCCTATGGCAAACGTTGCCGCCGGAGCGATAATCGCCGGCATAGCTAACAAATAATGCAGGTTTCTGAACAGTCAGTTTGCGCCGACCGGGCaa                                                                                                         >  2:715942/1‑235 (MQ=255)
                                                                                                                                 tttAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑CGa                                                                                                                                                                                                                          <  1:650671/107‑1 (MQ=255)
                                                                                                                                 tttAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAaata                                                                                                                                                 <  1:392734/180‑1 (MQ=255)
                                                                                                                                     atttGCCCCAGTATTTCGCACATGGTGCGGTGCTAAGTTTTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAATCTCGTACTTATAAATCCGCCAt                                                                                  <  1:1069184/236‑1 (MQ=255)
                                                                                                                                           cccAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGGGCCCGTGT‑TGGCGGGCATGAAAAGGGGTG‑CGACCGTTACTGGGGGGGTCGCCACTATGACAAACGGTGCCGCCCTAGCGATAATCGCCGGCATAGCTTAAAAAAAATGCAGGTTg                                                                                                                                       >  2:392734/1‑180 (MQ=255)
                                                                                                                                           cccAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGGGCCCCTGT‑TGGCGGGCATGAAAAGGGGTG‑CGACCCGTACTGGTGGGGTCGCCCCCATGACAAAAGTTGCCGCCGGCGCGCTAATCGCCGTCATAGCTAACAAAAAATTCCGGTTTCTTACCAGTCAGGTTGCGCTGACTGGGCAATCCCGTACTTTTAAAATCGcccaccccc                                                                             >  2:1380229/1‑231 (MQ=255)
                                                                                                                                                tATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTgc                                                                                                                                                                                     <  1:1129817/129‑1 (MQ=255)
                                                                                                                                                  tttCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCt                                                                                                                                     <  1:639364/175‑1 (MQ=255)
                                                                                                                                                  tttCGCACATGGGGCAGTGCTAAGTATTTTCGTCGGGGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGACCGGGTGTG‑CGACCGTTATTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGt                                                                                                            <  1:1021353/200‑1 (MQ=255)
                                                                                                                                                       cacaTGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTGGGGTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAATCTCGGACTTATAAATCCGCCATCGCCGCGCGACGA‑GTCgg                                                               >  1:494652/1‑239 (MQ=255)
                                                                                                                                                          aTGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGGGT‑TGGCGGGCATGAAAAGGGGGG‑CGACCGGTACTGGGGGGGTTGCCACTATGACAAAAGGTTGCGCCGGAGCg                                                                                                                                                                           >  2:1129817/1‑129 (MQ=255)
                                                                                                                                                            ggTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGGGGGCATGAAAAGGGGTG‑CGGCCGTTACTGGTGGGGTTGCCCCCTTGAAAAACGTTTCCGCCGGAGCGGAAATATCCGGCCCTGCTAACAAAAAATGCAGGTTTCTGACCAGTCCGTTTTTGCCGACTTGGCCATCTCGTa                                                                                                  >  2:1021353/1‑200 (MQ=255)
                                                                                                                                                            ggTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTGGGTTTGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACcagtcagt                                                                                                                          >  1:1019033/1‑176 (MQ=255)
                                                                                                                                                            ggTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACcagtcagt                                                                                                                          <  2:1019033/176‑1 (MQ=255)
                                                                                                                                                            ggTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGGGCCCGTGT‑TGGCGGGCATGAACAGGGGTG‑CGACCGGTTTTGGGGGGGTCGCCACCATGGAAAACGTTGCCGCCGGAGCGATAAATGCCGTCAGAACTAACAAAAAATGAAGTTTTCTGACcagtcag                                                                                                                           >  2:639364/1‑175 (MQ=255)
                                                                                                                                                              tGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCg                                                                                                                                                                    <  1:389755/132‑1 (MQ=255)
                                                                                                                                                              tGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCg                                                                                                                                                                    <  1:922504/132‑1 (MQ=255)
                                                                                                                                                              tGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTgcg                                                                                                                     <  1:427423/179‑1 (MQ=255)
                                                                                                                                                                 aGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTa                                                                                                                                                        <  1:701408/141‑1 (MQ=255)
                                                                                                                                                                   tGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑c                                                                                                                                                                                                                            <  1:904891/71‑1 (MQ=255)
                                                                                                                                                                   tGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAAc                                                                                                                                                      <  1:1294736/141‑1 (MQ=255)
                                                                                                                                                                    gCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGCGT‑TGGCGGGGCTGGACCGGGGGG‑CGGCCGGTACTGGGGGGGTTGCCCCTCCGGCAAACGGTGCCGCCGGGGCGGTAATCGCCGTCCCAGCTAACAAAAAATGCGGTTTGCTGACCAGTAAGGTTGCGCTGACCGGTCAATTCCGGACTTTTAAATCCGCCCTCGCCCCGCGCCGACTTCgggcgaacccccgc                                                   >  2:409525/1‑228 (MQ=255)
                                                                                                                                                                     cTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGt                                                                                                            <  1:756880/181‑1 (MQ=255)
                                                                                                                                                                     cTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAATCTCGTACTTATAAATCCGCCATCGCCGCGCGACGA‑TTCGGGAAGAACGCCAgg                                                  <  1:783569/238‑1 (MQ=255)
                                                                                                                                                                     cTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑CGACCGTTACTGGGGGGTTCGCCACTATGACAAACCTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTTCGCTGACTGGTCAATCTCGTACTTATAAATCCGCCATCGCCGCGCGACGA‑TTCGGGAAGAACGCCAAGc                                                 >  1:529639/1‑239 (MQ=255)
                                                                                                                                                                       aaGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGGCTGAACAGGGGGG‑CGGCCGTTACTGGTGGGTTCGCCCCCATGACAAACCTTGCCGCCGGAGGGATAAACGCCGGCCTAGCTAACAAAAAATGCAGTTTTCTGGCCCATCAGTTTTCGCTGACTGGTCAATCTCGTAATTTTAAATCCCCCATCCCCGCGTGCGCG‑TTTGGGGAGAAAGCCAGGCg                                                >  2:262825/1‑238 (MQ=255)
                                                                                                                                                                        aGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGGCCGTGT‑TGGCGGGCCTGAACCGGGGGG‑GGGCCGGTACTGGTGGGGTCGCCACTATGACAAACGTTGCCGCCGGAGCGATAATCGCCGGCATAGCTAAAAAAAAAAGCAGTTTTCTGACCAATCCGGTTTGGCTGACTGGTc                                                                                                           >  2:427423/1‑179 (MQ=255)
                                                                                                                                                                        aGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGGGGGGCTGAACAGGGGTG‑CGACCGTTACTGGGGGGGTCGCCACCATGGCCAACGTTGCCGCCGTAGCGGTAATCGCCGGCCTAGGTAAAAAATAAAGAAGTTTTCTGACCCGTCAGTTTGCGCTGACTGGGCAATCTCGTACTTATAAATCCCCCATCGCCGCCCGGCGC‑TTCCGGGAGAAAGCCAGGCgg                                               >  2:70210/1‑237 (MQ=255)
                                                                                                                                                                        aGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAATCTCGTACTTATAAATCCGCCATCGCCGCGCGACGA‑TTCGGGAAGAACGCCAGGCGTc                                              <  1:409525/239‑1 (MQ=255)
                                                                                                                                                                        aGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGGGTG‑CGGCCGGTACTGGGGGGGTCGCCACTATGACAAACGTTGCCGCCGGAGCGATAATCGCCGTCATAGCTAACAAAAAATGCAGGTTGCTGACCAGGCAGTTTTCGCTGACTGGTCAAACTCGTACTTATAAAACCCCCATCCCCCCCCGACGA‑TTCGGGGAGAAAGcccagcgg                                               >  2:782583/1‑232 (MQ=255)
                                                                                                                                                                        aGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAAC‑GGGGTG‑CGGCCGTTACTGGGGGGTTCGCCCCTATGACAAACGTTGCCGCCGTAGCGATAAACGCCGTCCTAGCTAACAAAAAATGCAGGTTTCTGACCAGGCAGGTTTCGCTGGCTGGTCAAACTCGTACTTTTAAAACCCCCCACTCCGCCCGGCCC‑ATCCGGAAGAAAGCCacgtggc                                              >  2:159691/1‑232 (MQ=255)
                                                                                                                                                                        aGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGGGT‑TGGCGGGGATGAACAGGGGTG‑CGGCCGGTACTGGGGGGGTCGCCCCTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGGCATAGCTAACAAATAATTCAGTTTTCTGGCCAGACCGGTTGCGCTGACTGGGCAATCTCGTACTTTTAAATCCGCCCTCCCCCCGCGGCCG‑TTTGGGAAGAAAGCCCGGGGt                                               >  2:438236/1‑238 (MQ=255)
                                                                                                                                                                        aGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGGGCCCGTGT‑TGGCGGGCATGAACAGGGGTG‑CGACCGTTACTGGGGGGGTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAATCTCGTACTTATAAATCCGCCCTCGCCGCGCGACGA‑TTCGGGAAGAACGCCAGGCGt                                               >  1:1038782/1‑238 (MQ=255)
                                                                                                                                                                        aGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGGGCCCCGGT‑TGGGGGGCATGAACAGGGGTG‑CGGCCGTTACTGGTGGGGTCGCCACCATGACAAACGTTGCCGCCGTAGCGATAATCGCCGGCATAGCTCACAAAAAATGCCGGTTGCTGACCAGTCAGTTTTCGCCGACTGGTCAAACTCGTACCTATAAATCCGCCATCGCCGCGCGACGA‑TTCGGGGAGAACGCCCGGGGt                                               >  2:783569/1‑238 (MQ=255)
                                                                                                                                                                          tATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑CGAc                                                                                                                                                                                                                         <  2:1393811/67‑1 (MQ=255)
                                                                                                                                                                          tATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCa                                                                                                                                            <  1:607317/144‑1 (MQ=255)
                                                                                                                                                                           aTTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAATCTCGTACTTATAAATCCGCCATCGCCGCGCGACGA‑TTCGGGAAGAACGCCAGGCGTccc                                            <  1:782583/238‑1 (MQ=255)
                                                                                                                                                                           aTTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAAC‑GGTGTG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAATCTCGTACTTATAAATCCGCCATCGCCGCGCGACGA‑TTCGGGAAGAACGCCAGGCGTccc                                            <  1:159691/237‑1 (MQ=255)
                                                                                                                                                                           aTTTTCGTCGGTGCCATTCTCTGGCTGGCGGGGGCCCGTGG‑TGGCGGGGATGAAAAGGGGTG‑GGACCGGTACTGGGGGGGTCGCCACTATGACAAAAGGTGCCGCCGGAGCGATAATCGCCGTCATAGCTAACAAAAAATg                                                                                                                                              >  2:701408/1‑141 (MQ=255)
                                                                                                                                                                             tttCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGGATGGAAAGGGGTG‑GGGCCGGTACt                                                                                                                                                                                                                  >  2:904891/1‑71 (MQ=255)
                                                                                                                                                                             tttCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCCTGT‑TGGCGGGGCTGAAAAGGGGTG‑GGGCCGGTTCTGGGGGGGTCGCCCCTATGACCAACGTTTCCCCCGTAGCGGTAATCCCCCGCCTAGCCAACaaaaaaatca                                                                                                                                            >  2:1294736/1‑137 (MQ=255)
                                                                                                                                                                               tcgtcgGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGGATGAAAAGGGGTG‑CGACCGGTTCTGGGGGGGTCGCCCCCTTGACAAACGTTGCCGCCGGAGCGATAATCGCCGTCCTAGCTAAAAAATAAAGCAGTTTTCTGACCAGGCAGTTTGCGCTCGCTGGTCAAACTCGTa                                                                                                  >  2:756880/1‑181 (MQ=255)
                                                                                                                                                                                   cgGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAATCTCGTACTTATAAATCCGCCATCGCCGCGCGACGA‑TTCGGGAAGAACGCCAGGCGTCCCGGGATCgg                                    <  1:438236/238‑1 (MQ=255)
                                                                                                                                                                                    ggTGCCATTCTCTGGCTGGCGGGGGCCCGTGT‑TGGCGGGGATGAACAGGGGTG‑CGACCGGTACtggt                                                                                                                                                                                                               >  1:1393811/1‑67 (MQ=255)
                                                                                                                                                                                    ggTGCCATTCTCTGGCTGGCGGGGGCCCGTGT‑TGGCGGGCATGGAAAGGGGTGCCGCCCGTT‑CTGGGGGGGTCGCCACTCTGACAAACGTTTCCCCCGTAGGGATAATCGCCGTCCTAGCTAACAAATAATGCAGGTTTCTGAc                                                                                                                                  >  2:607317/1‑144 (MQ=255)
                                                                                                                                                                                     gtgccgtcttctGGCTGGAGGGTGCCGGTGT‑TGGCGGGCATGAACAGGTGTG‑CGACCGTGACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCAGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGGCAATCTCGTACTTATAAATCCGCCATCGCCGCGCGACGA‑TTCGGGAAGAACGCCAGGCGTCCCGGGATCGGtt                                  <  2:1042701/229‑1 (MQ=255)
                                                                                                                                                                                     gtcccATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAATCTCGTACTTATAAATCCGCCATCGCCGCGCGACGA‑TTCGGGAAGAACGCCAGGCGTCCCGGGATCGGtt                                  <  1:1143041/235‑1 (MQ=255)
                                                                                                                                                                                              tcgggctgggGGGTGCCCGTGT‑TGGCGGGAAGGAAAGGGTGTG‑GGACGGTTACTGGTGGGTTCGCCAATATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAATCTCGTACTTATAAATCCGCCATCGCCGCGCGACGA‑TTCGGGAAGAACGCCAGGTGTCCCGGGATCGGTTGAATGCCAgc                        <  2:567204/236‑1 (MQ=255)
                                                                                                                                                                                                      gcGGGTGCCCCTGT‑TGGCGGGCATGAAAAGGGGTG‑GGACCGTTACTGGTGGGTTCGCCCCTATGGCAAAAGTTGGCGCCGTAGCGGTAATCGCCGGCATTGCTAACAAATAAAGCAGTTTTCTGACCCGTCAGTTTGCGCTGACTGGTCAAACTCGTACTTTTAAATCCGCCCTCCCCCCGCGGCGG‑TTCGGGGAGAAAGCCCGGCGTCCCGGGCTCGGTTGAATGCCacacggagcc                 >  2:411221/1‑229 (MQ=255)
                                                                                                                                                                                                               cGTGT‑TGGCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAATCTCGTACTTATAAATCCGCCATCGCCGCGCGACGA‑TTCGGGAAGAACGCCAGGCGTCCCGGGATCGGTTGAATGCCAGCGCGAGCCATAGTGCGCa       <  1:174621/239‑1 (MQ=255)
                                                                                                                                                                                                                     tGGCGGGCATGAACGGGTGGG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAATCTCGTACTTATAAATCCGCCATCGCCGCGCGACGA‑TTCGGGAAGAACGCCAGGCGTCCCGGGATCGGTTGAATGCCAGCGCGAGCCATAGTGCGCAGTGGc  <  2:1038782/239‑1 (MQ=255)
                                                                                                                                                                                                                      ggCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTGGGTTCGCCAc                                                                                                                                                                                                    >  2:470219/1‑45 (MQ=255)
                                                                                                                                                                                                                      ggCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTGGGTTCGCCAc                                                                                                                                                                                                    <  1:470219/45‑1 (MQ=255)
                                                                                                                                                                                                                      ggCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAATCTCGTACTTATAAATCCGCCATCGCCGCGCGACGA‑TTCGGGAAGAACGCCAgg                                                  >  1:300808/1‑190 (MQ=255)
                                                                                                                                                                                                                      ggCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAATCTCGTACTTATAAATCCGCCATCGCCGCGCGACGA‑TTCGGGAAGAACGCCAgg                                                  <  2:300808/190‑1 (MQ=255)
                                                                                                                                                                                                                                     gtgtG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGTAGTTTGCTGACCAGTCAGTTTGCGCTGa                                                                                                                 <  2:274096/113‑1 (MQ=255)
                                                                                                                                                                                                                                     gtgtG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGa                                                                                                                 >  1:274096/1‑113 (MQ=255)
                                                                                                                                                                                                                                     gtgtG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGa                                                                                                                 >  1:1268744/1‑113 (MQ=255)
                                                                                                                                                                                                                                     gtgtG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGa                                                                                                                 <  2:1268744/113‑1 (MQ=255)
                                                                                                                                                                                                                                     gtgtG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAAtctc                                                                                                     >  1:78870/1‑125 (MQ=255)
                                                                                                                                                                                                                                     gtgtG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAAtctc                                                                                                     <  2:78870/125‑1 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                               
CCCAGAAAGTTCCGGAGGTGACATATGAAACGCAAAAACGCTTCGTTACTCGGTAACGTGCTCATGGGGTTGGGTCTGGTGGTAATGGTGGTCGGCGTGGGGTATTCAATCCTCAACCAGTTACCACAGTTTAATATGCCCCAGTATTTCGCACATGGTGCAGTGCTAAGTATTTTCGTCGGTGCCATTCTCTGGCTGGCGGGTGCCCGTGT‑TGGCGGGCATGAACAGGTGTG‑CGACCGTTACTGGTGGGTTCGCCACTATGACAAACGTTGCCGCCGTAGCGATAATCGCCGTCATAGCTAACAAATAATGCAGTTTGCTGACCAGTCAGTTTGCGCTGACTGGTCAATCTCGTACTTATAAATCCGCCATCGCCGCGCGACGA‑TTCGGGAAGAACGCCAGGCGTCCCGGGATCGGTTGAATGCCAGCGCGAGCCATAGTGCGCAGTGGC  >  NC_000913/1258767‑1259217

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