Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,846,7050GC22.6% 113.5 / 10.5 62V236V (GTC→GTG) hycDformate hydrogenlyase subunit HycD
Reads supporting (aligned to +/- strand):  ref base G (15/33);  new base C (13/1);  total (28/34)
Fisher's exact test for biased strand distribution p-value = 4.52e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CAGACATGCTCAAACTCTCTTTAATCACGCCGCCAGCAAGGAGACGAACGCTAAAAATGCAAAGCCAAACCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGTTTCCATTTGT‑CCCCACGGAATAAACACCCCGACGAACATCTGCAACACCA‑CCAGCTGTTTCAGGCTGATACCCCA‑TTTCATG‑ACGCCAAAGCCGCTGCCGCTGTATT‑CAGAGAGCGGGCCTTCCTGCAGCTCCTGCTCGGCTTCCGCC‑AGGTCGAACGGCAGTTTGCCC‑ATTTCGATAAAGGTGGCGAACGCACAGGCACAAAGCGCCAGTACCAGCGGGATGCTCTGGCTCAGCGGCCAGTGATAAACGGTGTCGGTGATGTTGCTGATGT  >  NC_000913/2846465‑2846946
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       
cAGACATGCTCAAACTCTCTTTAATCACGCCGCCAGCAAGGAGACGAACGCTAAAAATGCAAAGCCAAACCCAGCCCAGGTCATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AATGCTTTCAATTGT‑CCCCCAGGAATAAAAACCCCCACCAACATcc                                                                                                                                                                                                                                               >  1:29016/1‑250 (MQ=255)
      tGCTCAAACTCTCTTTAATCACGCCGCCAGCAAGGAGACGAACGCTAAAAATGCAAAGCCAAACCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGCACAGCGCGATAACCAGGACGCCGACCACCAGTTT‑TACGCTGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGAG‑AAGGTTTCCCTCTGT‑CCCCCCAGACACAACACCCCGACGCCCATCTGcccac                                                                                                                                                                                                                                         >  2:90912/1‑247 (MQ=255)
          cAAACTCTCTTGAATCACGCCGCCAGCAAGGAGACGAACGCTAAAAATGCAAAGCCAAACCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGGCCACCGGCGGTG‑AAGGTTTCCATTTGT‑CCCCACGGAATAAACACCCCGACGAACATCTGCAAcacc                                                                                                                                                                                                                                       <  2:216635/249‑1 (MQ=255)
           aaaCACTCTTAAATCACGCCGCCAGCAAGGAGACGAACGCTAAAAATGCAAAGCCAAACCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGTTTCCATTTGT‑CCCCACGGAATAAACACCCCGACGAACATCTGCAACacca‑cc                                                                                                                                                                                                                                   <  2:348950/251‑1 (MQ=255)
               tctctTTAATCACGACGCCAGCAAGGAGACGCACGCTAAAAATGCAAAGCCAAACCCAGCCCAGGTACTGCGCGGAGTAATATCAAGAAGCAGACGGGCCATGCTGTTTTCGCACAGCGCGATAACCAGGACGCAGACCACCAGTTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGGTTCCCTTTGT‑CCCCACAGAATAAAAAACCCCCCGAAAATCTGCAACACCC‑CCAGCt                                                                                                                                                                                                                               >  1:430855/1‑251 (MQ=255)
                     taaacACGCAGCCAGCAAGGATACTAACGCTAAAAATTCAAAGCCAAACCCAGCCCAGGTAATGCGCGGAGTAAGATCAAGACGCAGACGGGCCATGCTGTTTGCGAACAGCGCGATAACCAGGACGGCGACCACCAGTTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGTTTCCATTTGT‑CCCCACGGAATAAACACCCCGACGAACATCTGCAACACCA‑CCAGCTGTTTc                                                                                                                                                                                                                          <  2:761/246‑1 (MQ=255)
                          aCGCCGCCAGCAAGTAGACGAACGCGAAAAATGCAAAGCCAAACCCAGCCCAGGTATTGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTGCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGGCCACCGTCGGTG‑AAGGTTTCCATTTGT‑CCCCACGGAATAAACACCCCGACTAACATCTGCAACCCCA‑CCAGCTGTTTCAGGCTg                                                                                                                                                                                                                    <  2:320046/251‑1 (MQ=255)
                                              aaCGCTAAAAATGCAAAGCCAAACCCAGCCCAGGTAATGCGAGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGCG‑AAGGGTTTCATTTGT‑CCCCCCGGGATAAACACCCCGACGAACATCTGCCACACCA‑CCACCTGTATTAGGCCGATTCCCCC‑TTTTATG‑aaccc                                                                                                                                                                                              >  1:310996/1‑247 (MQ=255)
                                                     aaaaTGCAAAGCCAAACCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAAGGTTCCCTTTGTTCCCCACG‑AATAAAAACCCCCACGGACAACTGCAACACCA‑CCAACCGTTTTAGGCTGATACCCCC‑TTTCATG‑AAGCCAAAaccg                                                                                                                                                                                       >  1:298387/1‑251 (MQ=255)
                                                      aaaTGCTAAGCCAAACCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGTTTCCATTTTT‑CCCCCCGGAATAAAAACCCCCACGAACATCTTCCACACCC‑CCCGCTGTTTCAGGCTGTTACCCCC‑TTTCATG‑ACGCCCAAaccgc                                                                                                                                                                                      >  1:137723/1‑251 (MQ=255)
                                                       aaTGCAAAGCCAAACCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGTTTCCATTTGT‑CCCCACGGAATAAACACCCCGACGAACATCTGCAACACCA‑CCAGCTGTTTCAGGCTGATACCCCA‑TTTCATG‑ACGCCAAAgccgct                                                                                                                                                                                     <  2:452254/251‑1 (MQ=255)
                                                             aaGCCAAACCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGTTTCCATTTGT‑CCCCACGGAATAAACACCCCGACGCACATCTGCAACACCC‑CCCGCTGGTTCAGGCTGATACCCCA‑TTTTATG‑ACGCCCAAaccgctgccccc                                                                                                                                                                               >  1:83014/1‑248 (MQ=255)
                                                               gCCAAACCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCCTGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGGTTCCATTTGT‑CCCCACGGAATAAACACCCCCACGGACATCTGCAACACCAACCA‑CCGGTTCAGGCTGATAACCCA‑TTTTATT‑ACGCCCAAACCGCCGGCGCTGt                                                                                                                                                                             >  1:324169/1‑251 (MQ=255)
                                                                  aaaCCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGTTTCCATTTGT‑CCCCACGGAATAAACACCCCGACGAACATCTGCAACACCA‑CCAGCTGTTTCAGGCTGATACACCA‑TTTCATG‑ACGCCAAAGCCGCTGCCGCTGTAtt                                                                                                                                                                          <  2:137723/251‑1 (MQ=255)
                                                                    aCCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCGGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGGTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGTTTCCATTTGT‑CCCCACGGAATAAACACCCCGACGAACATCTGCAACACCA‑CCAGCTGTTTCAGGCTGATACCCCA‑TTTCATG‑ACGCCAAAGCCGCTGCCGCTGTATT‑Ca                                                                                                                                                                       <  2:431281/251‑1 (MQ=255)
                                                                        agcccagGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGTTTCCATTTGT‑CCCCACGGAATAAACACCCCGACGAACATCTGCAACACCA‑CCAGCTGTTTCAGGCTGATACCCCA‑TTTCATG‑ACGCCAAAGCCGCTGCCGCTGTATT‑Cagaga                                                                                                                                                                   <  1:90912/251‑1 (MQ=255)
                                                                         gcccagGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGAG‑AAGGTTTCCATTTGT‑CCCCACCGAATAAAAAACCCCACGCACATCTTCAACACCA‑CCACCCGTTTCAAGCTGATACCCCC‑TTTTATGCAC‑CCAAAACCGCCGCCGCCGTATT‑ccgcg                                                                                                                                                                    >  1:52486/1‑245 (MQ=255)
                                                                              ggTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGTTTCCATTTGT‑CCCCACGGAATAAACACCCCGACGAACATCTGCAACACCA‑CCAGCTGTTTCAGGCTGAt                                                                                                                                                                                                                  <  1:382831/201‑1 (MQ=255)
                                                                               gTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGTTTCCATTTGT‑CCCCACGGAATAAACACCCCGACGAACATCTGCAACACCA‑CCAGCTGTTTCAGGCTGATACCCCA‑TTTCATG‑ACGCCAAAGCCGCTGCCGCTGTATT‑CAGAGAGCGGGcc                                                                                                                                                            <  1:399820/251‑1 (MQ=255)
                                                                                     cgcgGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGTTTCCATTTGT‑CCCCACGGAATAAACACCCCGACGAACATCTGCAACACCA‑CCAGCTGTTTCAGGCTGATACCCCA‑TTTCATG‑ACGCCAAAGCCGCTGCCGCTGGATT‑CAGATAGCGGGCCTTCCTg                                                                                                                                                      <  2:298387/251‑1 (MQ=255)
                                                                                       cgGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGTTTCCATTTTT‑CCCCCCGGAATAAACACCCCCACGAACATCTTCAAAACCC‑CCCGCTGTTTCAGGCTTATATCCCC‑TTTTATG‑ACGCCCACAACGCTGCCGCTGTGTT‑CAGAGAGCGGGCCTTCCtccc                                                                                                                                                    >  1:17659/1‑248 (MQ=255)
                                                                                           gTAATATCAAGACGCCGACGGGCCATGCTGTTTGCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGTTTCCATTTGT‑CCCCACGGAATAAACACCCCGACGAACATCTGCAACACCG‑CCAGCTGTTTCAGGCTGATACCCCA‑TTCCATG‑ACGCCAAAGCCGCTGCCGCTGTATT‑CAGAGAGCGGGCCTTCCTGCAGCTc                                                                                                                                                <  2:133282/251‑1 (MQ=255)
                                                                                            tAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGAG‑GAGGGTTCCATTTTT‑CCCCACCGAGTAAAAACCCCCGCGACAAACCACAACACCA‑CCAGCTGTTTCCGTCTGATCCCCCA‑ATTAATG‑ACGCCAACAACGCAGCCGCTGGATT‑Cagaagacgggcgatcccgcacgtcc                                                                                                                                               >  2:319078/1‑229 (MQ=255)
                                                                                                atCAAGACGCAGACGGGCTATGCTGTTTACGAACAGCGCGATAACCAGGACGCCGAACACCAGTTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGTTTCCATTTGT‑CCCCACGGAATAAACACCCCGACGAACATCTGCAACACCA‑CCAGCTGTTTCAGGCGGATACCCCA‑TTTCATG‑ACGCCAAAGCCGCTGCCGCTGTAGT‑AAGAGAGCGGGACTTCCTGCAGCTCCTGc                                                                                                                                            <  2:398586/250‑1 (MQ=255)
                                                                                                      acgcagtagtgccaagcgttagtCGAACAGCACGATAACCCGGACGCCGACC‑CCTCTTTGTACGAGGGCAAGCACCAGCGCCAGCAGCTGTCCGCCGGCGATT‑AAAGTTTCCATTTGT‑CCCCACGGAATAAACACCCCGACGAACAGCTGTAACACCT‑CCAGCTTCTACAGGCTGAGTCCCCA‑TGACATG‑TCCCCAAATCCTCTGACGCGGTATT‑CAGTGCGCGGGCCTTCGTGCAGCTCCTGCGCGGCtt                                                                                                                                     <  2:211175/229‑1 (MQ=255)
                                                                                                           gACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGGG‑AAGGGTTCCATTTGT‑CCCCCCGGGATAAAAACCCCCAAGAAAAACCGCAACACCC‑CCCGCCGTTTCCGGCTGATACCCCA‑ATTCACG‑GCGCCACAACCGCTGCCGCTGTATT‑CCGAGAGCGGGCCTTCCTGCAGCCCCCGCCCCGCTTCCGcc                                                                                                                                >  1:254984/1‑251 (MQ=255)
                                                                                                                   aTGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGTTTCCATTTGT‑CCCCACGGAATAAACACCCCGACGAACATCTGCAACACCA‑CCAGCTGTTTCAGGCTGATACCCCA‑TTTCa                                                                                                                                                                                                      <  2:46786/175‑1 (MQ=255)
                                                                                                                   aTGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGTTTCCATTTGT‑CCCCACGGAATAAACACCCCGACGAACATCTGCAACACCA‑CCAGCTGTTTCAGGCTGATACCCCA‑TTTCATG‑ACGCCAAAGCCGCTGCCGCTGTATT‑CAGAGAGCGGGCCTTCCTGCAGCTCCTGCt                                                                                                                                           <  2:364965/232‑1 (MQ=255)
                                                                                                                   aTGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGGTTCCATTTTT‑CCCCCCCGAATAAAAACCCCGACGAACATCTGCAACACCA‑CCAGCTGTTTTAGGCTGATACCCCC‑ttttc                                                                                                                                                                                                      >  1:46786/1‑173 (MQ=255)
                                                                                                                   aTGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGGTTCCATTTGT‑CCCCACGGAATAAAAACCCCGACGAACATCTGCAACACCC‑CCAGCTGTTTCAGGCTGATACCCCA‑TTTCATG‑ACGCCCAAACCGCTGCCCCTGTATT‑CAGAGAGCGGGCCTTCCTGCAGCCCCTGCtagatcgggagagcacacgc                                                                                                                        >  1:364965/1‑232 (MQ=255)
                                                                                                                        gTTTTCGAACAGCGCGATAACCAGGACGCCGACCCCCAGTTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGTTTCCATTTGT‑CCCCACGGAATAAACACCCCGACGAACATCTGCAACACCA‑CCAGCTGTTTCAGGCTGATACCCCA‑TTTCATG‑ACGCCAAAGCCGCTGCCGCTGTATT‑CAGAGAGCGGGCCTTCCTGCAGCTCCTGCTCGGCTTCCGCC‑AGGTCGAACGGc                                                                                                                   <  2:100245/250‑1 (MQ=255)
                                                                                                                                     ccgATAACCAGGACGCCGACCACCAGTTT‑TACGATGGCAATCTCCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGTTTCCATTTGT‑CCCCTCGGAATAAACCCCCCGACGAACATCTGCAACACCA‑CCAGCTGTTTCAGGCTGATACCCCA‑TGTCATG‑ACGCCAAAGCCGCTGCCGCTGTATT‑CAGAGAGCGGGCCTTCCTGCAGCTCCTGCTCGGCTTCCGCC‑AGGTCGAACGGCAGTTTGCCC‑ATTTc                                                                                                    <  2:367585/250‑1 (MQ=255)
                                                                                                                                          aaCCAGGACGCCGACCCCCAGTTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGTTCCCATTTGT‑CCCCACGGAATAAACACCCCGACGAACATCTGCAACACCA‑CCAGCTGTTTCAGGCTGATACCCCA‑TTTCATG‑ACGCCAAAGCCGCTGCCGCTGTATT‑CAGAGAGCGGGCCTTCCTGCAGCTCCTGCTCGGCTTCCGCC‑AGGTCGAACGGCAGTTTGCCC‑ATTTCGATaa                                                                                               <  1:458146/251‑1 (MQ=255)
                                                                                                                                           aCCAGGACGCCGACCACCAGTTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGTTTCCATTTGT‑CCCCACGGAATAAACACCCCGACGAACATCTGCAACACCA‑CCAGCTGTTTCAGGCTGATACCCCA‑TTTCATG‑ACGCCAAAGCCGCTGCCGCTGTATT‑CAGAGAGCGGGCCTTCCTGCAGCTCCTGCTCGGCTTCCGCC‑AGGTCGAACGGCAGTTTGCCC‑ATTTCGATaaa                                                                                              <  1:17047/251‑1 (MQ=255)
                                                                                                                                                 aCGCCGACCACCAGTTT‑TACGATGGCAATCACCAGCGCCAGCACCAGTCCACAGGCGGTG‑AAGGTTTCCATTTGT‑CCCCACGGGATAAACACACCCCCGAACATATGCAAAAACA‑CCAGCCGTTTAAGGCTGATACCCCC‑TTTCATT‑ACGCCCAACACGCCGCCCCTGTATT‑CAGCGAGGGGGCCTTCCTGCAGCTCCTGCCCGGCCTCCCCC‑AGGGCGGACGGCAGTTTGccccattctgtaaagggggg                                                                                         >  1:159336/1‑234 (MQ=255)
                                                                                                                                                    ccgtcccccaGTTG‑TACGAGAGCAATCACCATCGCCAGCAGCAGTCCGCCGGCGGTT‑AAGGTTTCCATTTGT‑CCCCACGGAATAAACCCCCCCCCGACCATCTGCAACACCA‑CCAGCTGTTTCAGGCTGATACCCCA‑CTTCATT‑ACGCCAAAGCCGCCGCCGCTGTATT‑CAGAGAGCGGGCCTTCCTGCTGCGCCTGCTCGGCTTCCGCC‑AGGTCGAACGGCGGTTTGCCC‑TTTTCGATAATGGTTGCGaa                                                                                     <  2:226101/244‑1 (MQ=255)
                                                                                                                                                    ccGACCACCAGTTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGTTTCCATTTGT‑CCCCACGCAAAAAAAAACCCGACGGACATCTGCAAAACCA‑CCAGACGTTTCAGGCTGATACCCCA‑TTTCATG‑ACGCCAAAGCCGCTGCCGCTGTATT‑CACAGAACGGGGCCTCCTGCAGCTCCTGCTCCGCGTCCGCC‑CGGGCGAACGGCCGGTTGCCC‑CTTTCGATAAAGGAGGCgg                                                                                      >  1:473583/1‑249 (MQ=255)
                                                                                                                                                         caccaTTTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGTTTCCATTTGT‑CCCCACGGAATAAACACCCCGACGAACAGCTGCAACACCA‑CCAGCTGTTTCAGGCTGATACCCCA‑TTTCATG‑ACGCCAAAGCCGCTGCCGCTGTATT‑CAGAGAGCGGGCCTTCCTGCAGCGCCTGCTCGGCTTCCGCC‑AGGTCGAACGGCAGTTTGCCC‑ATTTCGATAAAGGTGGCGAACGcac                                                                                <  2:231603/251‑1 (MQ=255)
                                                                                                                                                         caccaGTTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGGAGGGTTTCC‑TTTTT‑CCCCCCCGAATAAACACCCCCACGAACACCTGCAACACCC‑CCACCCGTTTCAGGCTGATACCCCCCTTTC‑TG‑ACGCCCAAACCCCCGCCCCTGTATT‑CCGAGAGCGGGCCTTCCCTCAGCTCCCGCCCCGCCTCCGCC‑CGGGCGGACGGCAGGTTGCCCCATTT‑GCTAAAGGTgggggacccac                                                                                >  2:103383/1‑242 (MQ=255)
                                                                                                                                                           ccaGTTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGTTTCCATTTGT‑CCCCACGGAATAAACACCCCGACGAACATCTGCAACACCA‑CCAGCTGTTTCAGGCTGATACCCCA‑TTTCATG‑ACGCCAAAGCCGCTGCCGCTGTATT‑CAGAGAGCGGGCCTTCCTGCAGCTCCTGCTCGGCTTCCGCC‑AGGTCGAACGGCAGTTTGCCC‑ATTTCGATAAAGGTGGCGAACGCACAg                                                                              <  1:25195/251‑1 (MQ=255)
                                                                                                                                                             aGTTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGTTTCCATTTGT‑CCCCACGGAATAAACACCCCGACGAACATCTGCAACACCA‑CCAGCTGTTTCAGGCTGATACCCCA‑TTTCATG‑ACGCCAAAGCCGCTGCCGCTGTATT‑CAGAGAGCGGGCCTTCCTGCAGCTCCTGCTCGGCTTCCGCC‑AGGTCGAACGGCAGTTTGCCC‑ATTTCGATAAAGGTGGCGAACGCACAGGc                                                                            <  1:353115/251‑1 (MQ=255)
                                                                                                                                                                tt‑tACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGTTTCCATTTGT‑CCCCACGGAATAAACACCCCGACGAACATCTGCAACACCA‑CCAGCTGTTTCAGGCTGATACCCCA‑TTTCATG‑ACGCCAAAGCCGCTGCCGCTGTATT‑CAGAGAGCGGGCCTTCCTGCAGCTCCTGCTCGGCTTCCGCC‑AGGTCGAACGGCAGTTTGCCC‑ATTTCGATAAAGGTGGCGAACGCACAGGcac                                                                          <  1:103383/250‑1 (MQ=255)
                                                                                                                                                                    aCGATGGCACTCACCAGCGCCAGCAGTAGACCGCCGGCGGTG‑AAGGTTTCCATTTGT‑CCCCACGGAATAAACACCCCGACGAGCATCTGCAACACCA‑CCAGCTGTTTCAGGTTGATACCCCA‑TTTCATG‑ACGCCAAAGCCGCTGCCGCTGTAGT‑CAGGTGGCGGGTCTTTCTGCAGCTCCTGGACTGGTTCCGGG‑AGGTTGAACGTCAGGTTGCCC‑AGTTCGTTATAGGGGGCGAACGCGCTGGCTCAAAg                                                                      <  2:232359/251‑1 (MQ=255)
                                                                                                                                                                      gATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGTTTCCATTTGT‑CCCCACGGAATAAACACCCCGACGAACATCTGCAACACCA‑CCAGCTGGTTCAGGCTGATACCCCA‑TTTCATG‑ACGCCAAAGCCGCTGCCGCTGTATT‑CAGAGAGCTGGCCTTCCTGCAGGTCCTGCGCGGCTTCCGCC‑AGGTCTAACGGCAGTTTGCCC‑ATTTCGAGAAAGGTGGCGAACGCACAGGCACAAAg                                                                      <  2:299124/249‑1 (MQ=255)
                                                                                                                                                                         ggCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGTTTCCATTTGT‑CCCCACGGAATAAACACCCCGACGAACATCTGCAACACCA‑CCAGCTGTTTCAGGCTGATACCCCA‑TTTCATG‑ACGCCAAAGCCGCTGCCGCTGTATT‑CAGAGAGCGGGCCTTCCTGCAGCTCCTGCTCGGCTTCCGCC‑AGGTCGAACGGCAGTTTGCCC‑ATTTCGATAAAGGTGGCGAACGCACAGGCACAAAGCGCCa                                                                 <  1:429021/251‑1 (MQ=255)
                                                                                                                                                                              ccACCATCGCCCGCAGCAGTCCACCGGCTGGG‑AAGGTTTCCATTTGT‑CCCCACGGAAGAAACACCCCGACGAACATCTGCAACACCA‑CCAGCTGTTTCAGGCTGATACCCCA‑TTGCATG‑ACGCCAAAGCCGCTGCCGCTGTATT‑CAGAGAGCGGGCCTTCCTGCAGCTCCTGCTCGGCTTCCGCC‑AGGTCGAACGGCAGTTTGCCC‑ATTTCGATAAAGGTGGCGAACGCACAGGCACAAAGCGCCAGttcc                                                            <  2:181885/250‑4 (MQ=255)
                                                                                                                                                                                aCCAGCGCCAGCAGCAGTCCACCGGCGGTG‑GAGGTTTCCCTTTGT‑CCCCCCGGGATAAAAACCCCCACGAAACTCTTCAACCCCC‑CCAGCTGTTTCAAGCTGATACCCCC‑TTTCATG‑ACGCCCAAGCCGCCGCCGCTGTATT‑CACAGAGCGGGCCCTCCTGCCGCTCCTTCTCGGCTTCCGCC‑CGGTTGGACGGCCGTTTGCCC‑ATTTCGCTAAAAGGGGCGGACGAACAAGCACAAAAGGCCAGTaaca                                                           >  1:102035/1‑247 (MQ=255)
                                                                                                                                                                                 ccAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGTTTCCATTTGT‑CCCCACGGAATAAACACCCCGACGAACATCTGCAACACCA‑CCAGCTGTTGCAGGCTGATACCCCA‑TTTCATG‑ACGCCAAAGCCGCTGCCGCTGTATT‑CAGAGAGCGGGCCTTCCTGCAGCTCCTGCTCGGCTTCCGCC‑AGGTCGAACGGCAGTTTGCCC‑ATTTCGATAAAGGTGGCGAACGCACAGGCACAAAGCGCCAGTACCAGc                                                         <  1:372455/251‑1 (MQ=255)
                                                                                                                                                                                 ccAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGGTTCCACTTTT‑TCCCACGGAAAAAAAACCCCCACGAAAATCTTCAAAACCA‑CCCGCTGTTTCAGGGTGATACCCCA‑TTTTTTT‑AAGCCAAAGCCGCTGGCGCCGGATT‑CAGAGAGCGGGCCTTCCCGCCGCTCCCCCTCGGCTTCCGCCCAGCGCGA‑CGGCAGCTTGCCC‑CTTTTGATAAAGGCGGCGGACGCCCAGGCACAAAAGGCCAGTaacacc                                                         >  1:16469/1‑246 (MQ=255)
                                                                                                                                                                                  cAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGTTTTCCTTTGT‑CCCCCCCGGATAAACACCCCGCCGAACAACTGCCACACCA‑CCACCCGTTTCCGGCCGATACCCCC‑TTTCATG‑ACGCCCAAACCCCCCCCGCCGTATTTCA‑AGAGCGGGCCTTCCTGCCGCCCCCTCCCGGCCTCCTCC‑CGGTCCAACGGCCGTTTGCCC‑CTTTCGCTAAAAGTGGGGAACGCCCAGGCCCAAAGCGCCCGTAACAGCg                                                        >  2:80871/1‑251 (MQ=255)
                                                                                                                                                                                                   cACCGGCGGTG‑AAGGTTGCCATTTGT‑CCCCACGGAATAAACCCCCCGACTAACATCTGCAACACCA‑CCAGTTGTTTCAGGCTGATCCCCCA‑TTTCATG‑ACGCCAAAGCCGCTGCCGCTGTATT‑CAGAGAGCTGGCCTTCCTGCAGCTCCTGCTCGGCTTCCGCC‑AGGTCGAACGGCAGTTTGCCC‑ATTTCGATAAAGGTGTCGAACGCACAGGCACAAAGCGCCAGTTCCAGCGGGATGCGCGGGCTCAGc                                       <  2:430855/251‑1 (MQ=255)
                                                                                                                                                                                                          ggTG‑AAGGTTTCCATTTGT‑CCCCACGGAATAAACACCCCGACGAACATCTGCAACACCA‑CCAGCTGTTTCAGGCTGATACCCCA‑TTTCATG‑ACGCCAAAGCCGCTGCCGCTGTATT‑CAGAGAGCGGGCCTTCCTGCAGCTCCTGCTCGGCTTCCGCC‑AGGTCGAACGGCAGTTTGCCC‑ATTTCGATAAAGGTGGCGAACGCACAGGCACAAAGCGCCAGTACCAGCGGGATGCTCTGGCTCAGCGGCCAGt                                >  1:205867/1‑251 (MQ=255)
                                                                                                                                                                                                             g‑AAGGTTTCAATTTGA‑CCCAACGGAATAAACACCCCGACGAACATCTGCAACACCA‑CCAGCTGTTTCAGGCTGATACCCCA‑TTTCATG‑ACGCCAAAGCCGCTGACGCTGTATT‑CAGAGAGCGGGCCTTCCTGCAGCTCCTGCTCGGCTTCCGCC‑AGGTCGAACGGCAGTTTGCCC‑ATTTCGATAAAGGTGGCGAACGCACCAGCCCAAAGCGACAGGACCAGCGGGACGCTCTGGCTCAGCGGCCAGTGAt                             >  2:281827/1‑251 (MQ=255)
                                                                                                                                                                                                                 ggTTTCCATTTGT‑CCCCTCGGAATAAACACCCCGACGAACATCTGCAACACCA‑CCAGCTGTTTCAGGCTGATACCCCA‑TTTCATG‑ACGCCAAAGCCGCTGCCGCTGTATT‑CAGAGAGCGGGCCTTCCTGCAGCTCCTGCTCGGCTTCCGCC‑AGGTCGAACGGCAGTTTGCCC‑ATTTCGATAAAGGTGGCGAACTCACAGGCACAAAGCTCCAGTACCAGCGGGATGCTCTGGCTCAGCGGCCAGTGATaaa                          <  1:208790/251‑1 (MQ=255)
                                                                                                                                                                                                                 ggTTTCCATTTGT‑CCCCACGGAATAAACACCCCGACGAACATCTGCAACACCA‑CCAGCTGTTTCAGGCTGATACCCCA‑TTTCATG‑ACGCCAAAGCCGCTGCCGCTGTATT‑CAGAGAGCGGGCCTTCCTGCAGCTCCTGCTCGGCTTCCGCC‑Agg                                                                                                                            <  1:96851/154‑1 (MQ=255)
                                                                                                                                                                                                                 ggTTTCCATTTGT‑CCCCACGGAATAAACACCCCGACGAACATCTGCAACACCA‑CCAGCTGTTTCAGGCTGATACCCCA‑TTTCATG‑ACGCCAAAGCCGCTGCCGCTGTATT‑CAGAGAGCGGGCCTTCCTGCAGCTCCTGCTCGGCTTCCGCC‑Agg                                                                                                                            >  2:96851/1‑154 (MQ=255)
                                                                                                                                                                                                                       cATTTGT‑CCCCACGGAATAAACACCCCGACGAACATCTGCAACACCA‑CCAGCTGTTTCAGGCTGATACCCCA‑TTTCATG‑ACGCCAAAGCCGCTGCCGCTGTATT‑CAGAGAGCGGGCCTTCCTGCAGCTCCTGCTCGGCTTCCGCC‑AGGTCGAACGGCAGTTTGCCC‑ATTTCGATAAAGGTGGCGAACGCACAGGCACAAAGCGCCAGTACCAGCGGGATGCTCTGGCTCAGCGGCCAGTGATAAACGgtgt                    >  1:438457/1‑251 (MQ=255)
                                                                                                                                                                                                                       cATTTGT‑CCACACGGCATAAACAACACGACGAACATCCGCAACCCCA‑CCAGCTGTTTCAGCCTGATACCCCA‑TTTCATG‑AAGCCACAGCCGCTGCCGCTGTATT‑CCGACAGCGCGACTTCCCGCAGCTCCTGCTCGCCTTCCGCC‑AGGTAGACCGGCCGTGTGCCA‑CTTTCGCTAAAGGTGGCGAACGCACAGGCACAAAGCGCCAGTACCCGCGGGATGCTCTGGCTCAGCGGCCAgcgcacaacag                       >  2:21163/1‑238 (MQ=255)
                                                                                                                                                                                                                          ttGT‑CCCCACGGAATAAACACCCCGACGAACATCTGCAACACCA‑CCAGCTGTTTCAGGCTGATACCCCA‑TTTCATG‑ACGCCAAAGCCGCTGCCGCTGTATT‑CAGAGAGCGGGCCTTCCTGCAGCTCCTGCTCGGCTTCCGCC‑AGGTCGAACGGCAGTTTGCCC‑ATTTCGATAAAGGTGGCGAACGCACAGGCACAAAGCGCCAGTACCAGCGGGATGCTCTGGCTCAGCGGCCAGTGATAAACGGTGTCgg                 >  2:74122/1‑251 (MQ=255)
                                                                                                                                                                                                                               ccccACGGAATAAACACCCCGACGAACATCTGCAACACCA‑CCAGCTGTTTCAGGCTGATACCCCA‑TTTCATG‑ACGCCAAAGCCGCTGCCGCTGTATT‑CAGAGAGCGGGCCTTCCTGCAGCTCCTGCTCGGCTTCCGCC‑AGGTCGAACGGCAGTTTGCCC‑ATTTCGATAAAGGTGGCGAACGcaca                                                                               <  1:480970/185‑1 (MQ=255)
                                                                                                                                                                                                                               ccccACGGAATAAACACCCCGACGAACATCTGCAACACCA‑CCAGCTGTTTCAGGCTGATACCCCA‑TTTCATG‑ACGCCAAAGCCGCTGCCGCTGTATT‑CAGAGAGCGGGCCTTCCTGCAGCTCCTGCTCGGCTTCCGAC‑AGGTCGAACGGCAGTTTGCCC‑ATTTCGATAACGGTGGCGAACGcaca                                                                               >  2:480970/1‑185 (MQ=255)
                                                                                                                                                                                                                                          aaaCACCCCGACGAACATCTGCAACACCA‑CCAGCTGTTTCAGGATGATACCCCA‑TTTAATG‑ACGCCAAAGCCGCTGCCGCTGTATT‑CAGAGCGCGGGCCTTCCTGCAGCTCCTGCTCGGCTTCCGCC‑AGGTCGAACGGCAGTTTGCCC‑ATTTCGATACAGGTGGCGAACGCACAGGCACACAGCGCCAGTACCAGCGGGATGCTCTGGCTCAGCGGCCAGTGATACACGGTGTCGGTGCTGTTGCTGATGt  >  2:71798/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       
CAGACATGCTCAAACTCTCTTTAATCACGCCGCCAGCAAGGAGACGAACGCTAAAAATGCAAAGCCAAACCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTT‑TACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTG‑AAGGTTTCCATTTGT‑CCCCACGGAATAAACACCCCGACGAACATCTGCAACACCA‑CCAGCTGTTTCAGGCTGATACCCCA‑TTTCATG‑ACGCCAAAGCCGCTGCCGCTGTATT‑CAGAGAGCGGGCCTTCCTGCAGCTCCTGCTCGGCTTCCGCC‑AGGTCGAACGGCAGTTTGCCC‑ATTTCGATAAAGGTGGCGAACGCACAGGCACAAAGCGCCAGTACCAGCGGGATGCTCTGGCTCAGCGGCCAGTGATAAACGGTGTCGGTGATGTTGCTGATGT  >  NC_000913/2846465‑2846946

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