Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009641,778,3950TG36.8% 36.5 / 11.3 38V20G (GTT→GGT) mutLDNA mismatch repair endonuclease
Reads supporting (aligned to +/- strand):  ref base T (3/21);  new base G (14/0);  total (17/21)
Fisher's exact test for biased strand distribution p-value = 7.03e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCCGGCACAGCTGTCCTTTTTTGACGAAGCGGAAAAGCCGGCTGAAACACCAAAGCTTTCAAAAAAAGAAAAGCAAGTGATCGATGCTTTCAAATCACTTAATATACTGGATATGACACCGCTTGAAGCGATGAATGAAATGTACAAGCTGCAAAAGAAATTACATTAAAACGGGGTGATGAATGTGGCAAAAGTCATCCAACTGTCAGATGAGCTTTCAAATAAAATAGCGGCGGGCGAGGTTGTGG‑AACGGCCCGCCTCAGTCGTC‑AAAGAATTGGTGGAAAATGCG‑ATCGACGCTGACAGCACAGTCATTGAAATCGATATTGAGGAAGCAGGCCTTGCATCCATTCGGGTCTTGG‑ATAACGGCGAAGGAATGGAGAATGAAGATTGCAAGCGAGCTTTCC‑GGCGCCACGCAACGAGTAAAATAAAAGATGAAAATGATTTATTCAGGGTGAGAACGCTCGGTTTCAGGGGAGAAGCCC  >  NC_000964/1778153‑1778634
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
gCCGGCACAGCGGTCCTTTTTTGACAAAGCGGAAAAGCCGGCTGAAACACCAAAGCTTTCAAAAAAAGAAAAGCAAGTGATCGATGCTTTCAAATCACTTAATATACTGGATATGACACCGCTTGAAGCGATGAATTAAATGTACAAGCTGTAAAATAAATTACATTAAAACGTGGGGATGAATGTGGCAAAAGTCGACCAGCGGTCAGATGAGCTTTCAAATAAAATATCGGCGGGGGAGGTTGTGG‑AAc                                                                                                                                                                                                                                             <  2:95939/251‑1 (MQ=255)
                 tttttGACGAAGCGGAAAAGCCGGCTGAAACACCAAAGCTTTCAAAAAAAGAAAAGCAAGTGATCGATGCTTTCAAATCACTTAATATACTGGATATGACACCGCTTGAAGCGATGAATGAAATGTACAAGCTGCAAAAGAAATTACATTAAAACGGGGTGATGAATGTGGCAAAAGTCATCCACCTGTCAGATGAGCTTTCAAATAAAATAGCGGCGGGCGAGGTTGTGG‑AACGGCCCGCCTCAgtcgt                                                                                                                                                                                                                             <  1:162207/250‑1 (MQ=255)
                  ttttGACGAAGCGGAAAAGCCGGCTGAAACACCAAAGCTTTCAAAAAAAGAAAAGCAAGTGATCGAGGCTTTCAAATCACTGAATATACTGGATATGACACTGCTTGAAGCGATGAATGAAATGTACAAGCTGCAAAAGAAATTACATTAAAACGGGGTGATGAATGTGGCAAAAGTCATCCAACTGTCAGATGTGCTTTCAAATAAAATAGCGGCGGGCGAGGTTGTGG‑AACGGCCCGCCTCAgtcgtc                                                                                                                                                                                                                            <  1:271241/250‑1 (MQ=255)
                       aCGAAGCGGAAAAGCCGGCTGAAACACCAAAGCTTTCAAAAAAAGAAAAGCAAGTGATCGATGCTTTCAAATCACTTAATATACTGGATATGACACCGCTTGAAGCGATGAATGAAATGTACAAGCTGCAAAAGAAATGACATTAAAACGGGGTGATGAATGTGGCAAAAGTCATCCAACTGTAAGATGAGCTTTCAAATAAAATAGCGGCGGGGGGGGTCGTGG‑AGCGGCCCTCCTCAgtcggataaga                                                                                                                                                                                                                       >  2:25514/1‑243 (MQ=255)
                               gAAAAGCCGGCTGAAACACCAAAGCTTTCAAAAAAAGAAAAGCAAGTGATCGATGCTTTCAAATCACTTAATATACTGGATATGACACCGCTTGAAGCGATGAATGAAATGTACAAGCTGCAAAAGAAATTACATTAAAACGGGGTGATGAATGTGGCAAAAGTCATCCAACTGTCAGATGAGCTTTCAAATAAAATAGCGGCGGGCGGGGGGGGGG‑GAAGGCCCGCCTCcgtcgtgcaaaagatggtggg                                                                                                                                                                                                              >  2:63175/1‑236 (MQ=255)
                                           gAAACACCAAAGCTTTCAAAAAAAGAAAAGCAAGTGATCGATGCTTTCAAATCACTTAATATACTGGATATGACACCGCTTGAAGCGATGAATGAAATGTACAAGCTGCAAAAGAAATTACATTAAAACGGGGTGATGAATGTGGCAAAAGTCATCCAACTGTCAGATGAGCTTTCAAATAAAATAGCGGCGGGCGGGGGTGGGG‑GAAGGGCCGCCCCTGTCGGC‑AAAGAAATGGGGGAAAAAGCG‑Atagc                                                                                                                                                                                                >  1:6898/1‑248 (MQ=255)
                                              acacCAAAGCTTTCAAAAAAAGAAAAGCAAGTGATCGATGCTTTCAAATCACTTAATATACTGGATATGACACCGCTTGAAGCGATGAATGAAATGTACAAGCTGCAAAAGAAATTACATTAAAACGGGGTGATGAATGTGGCAAAAGTCATCCAACTGTCAGATGAGCTTTCAAATAAAATAGCGGCGGGGGAGGGTGTGG‑GACGGCCCCCCCCCGTCGTC‑CAAGAATTGGGGGGAAAAACG‑ATCGACGc                                                                                                                                                                                             >  1:261817/1‑251 (MQ=255)
                                                           cAAAAAAAGAAAAGCAAGTGATCGATGCTTTCAAATCACTTAATATACTGGATATGACACCGCTTGAAGCGATGAATGAAATGTACAAGCTGCAAAAGAAATTACATTAAAACGGGGTGATGAATGTGGCAAAAGTCATCCAACTGTCAGATGAGCTTTCAAATAAAATAGCGGCGGGCGAGGGTGTGGGAAAGGCC‑GCCTCACGCGTC‑CAAGAATTGGTGGGAAAAGCG‑GTCGACGCCGCCAACACAGTc                                                                                                                                                                                >  1:247894/1‑251 (MQ=255)
                                                               aaaataaaagAAAGTGATCGATGCTTTCAAATCCCTTAATATACGGAATATGACCCCGCTGGAAGCGAGAAATGAAATGTAAAAGCTGCAAAAGAAATTACATTAAAACGGGGTGATAAATTTGGAAAAAGGCATCCC‑CTGTCAGATGAGTTTTCAAATAAAATAGCGGCGGGCGAGGTTGTGG‑AACGGCCCGCCTCAGTCGTC‑AAAGAATTGGTGGAAAATGCG‑ATCGACGCTGACAGCACAGTCATTGa                                                                                                                                                                           <  2:184442/251‑1 (MQ=255)
                                                                            gTGATCGATGCTTTCAAATCACTTAATATACTGGATATGACACCGCTTGAAGCGATGAATGAAATGTACAAGCTGCAAAAGAAATTACATTAAAACGGGGTGATGAATGTGGCAAAAGTCATCCAACTGTCAGATGAGCTTTCAAATAAAATAGCGGCGGGCGAGGTTGTGG‑AACGGCCCGCCTCAGTCGTC‑AAAGAATTGGTGGAAAATGCG‑ATCGACGCTGACAGCACAGTCATTGAAATCGATATTGa                                                                                                                                                               <  1:235759/251‑1 (MQ=255)
                                                                                          cAAATCACTTAATATACTGGATATGACACCGCTTGAAGCGATGAATGAAATGTACAAGCTGCAAAAGAAATTACATTAAAACGGGGTGATGAATGTGGCAAAAGTCATCCAACTGTCAGATGAGCTTTCAAATAAAATAGCGGCGGGCGAGGTTGTGG‑AACGGCCCGCCTCAGTCGTC‑AAAGAATTGGTGGAAAATGCG‑ATCGACGCTGACAGCACAGTCATTGAAATCGATATTGAGGAAGCAGGCCtt                                                                                                                                                  <  1:29750/250‑1 (MQ=255)
                                                                                                 cTTAATATACTGGATATGACACCGCTTGAAGCGATGAATGAAATGTACAAGCTGCAAAAGAAATTACATTAAAACGGGGTGATGAATGTGGCAAAAGTCATCCAACTGTCAGATGAGCTTTCAAATAAAATAGCGGCGGGGGAGGGTGTGG‑GACGGCCCGCCTCAGGCGGC‑CAAAAATTGGGGGGAAAAGAG‑GACGACGCTGACAGCACAAGCATTGAAATATAAAATGAGGGAGCAGGGCTTGCAACCAt                                                                                                                                          >  1:11464/1‑251 (MQ=255)
                                                                                                    aaTATACTGGATATGACACCGCTTGAAGCGATGAATGAAATGTACAAGCTGCAAAAGAAATTACATTAAAACGGGGTGATGAATGTGGCAAAAGTCATCCAACTGTCAGATGAGCTTTCAAATAAAATAGCGGCGGGCGGGGGTGGGG‑GACGGCCCGCCCCAATAGTC‑AAAAAATTGGGGGAAAAAGCG‑GTCGACGCTGACAGCACAGTCATTGAAATCGAAATTGAGGAAGGAGGCCTTGCATCCATTCg                                                                                                                                       >  1:119125/1‑251 (MQ=255)
                                                                                                     atatACTGGATATGACACCGCTTGAAGCGATGAATGAAATGTACAAGCTGCAAAAGAAATTACATTAAAACGGGGTGATGAATGTGGCAAAAGTCATCCAACTGTCAGATGAGCTTTCAAATAAAATAGCGGCGGGCGGGGGTGTGG‑AACGGCCCCCCCCCGTCGTC‑CAAGAATTGGGGGGAAATGGG‑GTCGACGCTGGCAGCACAGGCCTTTAAAACGATATTTAGGAAACAGGCCTTGCATCCATTCgg                                                                                                                                      >  1:24022/1‑251 (MQ=255)
                                                                                                        taCTGGATATGACACCGCTTGAAGCGATGAATGAAATGTACAAGCTGCAAAAGAAATTACATTAAAACGGGGTGATGAATGTGGCAAAAGTCATCCAACTGTCAGATGAGCTTTCAAATAAAATAGCGGCGGGCGGGGTTGGGG‑AAAGGCCCCCCTCAGTCGTC‑AAAGAATTGGTGGGAAAAGCG‑AACGACGCTGACAGCAAAGTCCTTGAAAACGGTATTGAGGAAACAGGCCTTGCCTCCATTCggggc                                                                                                                                   >  1:73178/1‑249 (MQ=255)
                                                                                                                       cGCTTGAAGAGATGAATGAAATGTACAAGCTGCAAAAGAAATTACATTAAAACGGGGTGATGAATGTGGCAAAAGTCATCCAACTGTCAGATGAGCTTTCAAATAAAATAGCGGCGGGCGAGGTTGTGG‑AACGGCCCGCCTCAGTCGTC‑AAAGAATTGGTGGAAAATGCG‑ATCGACGCTGACAGCACAGTCATTGAAATCGATATTGAGGAAGCAGGCCTTGCATCCATTCGGGTCTTGG‑ATAACGGCGaa                                                                                                                   <  1:137157/251‑1 (MQ=255)
                                                                                                                            gAAGCGATGAAGGAAAGGTACAAGCGGCAAAAGAAATTACATTAAAACGGGGTGATAAATGTGGAAAAAGTAATCCAACTGTCAGATGAGCTTTCAAATAAAATGGGGGCGGGCGAGGTTGTGG‑AACGGCCCGCCTCAGTCGTC‑AAAGAATTGGTGGAAAATGCG‑ATCGACGCTGACAGCACAGTCATTGAAATCGATATTGAGGAAGCAGGCCTTGCATCCATTCGGGTCTTGG‑ATAACGGCGAAGGAAt                                                                                                              <  2:193729/251‑1 (MQ=255)
                                                                                                                                     agtgaaATGTACAAGCTGCAAAAGAAATTACATTAAAACGGGGATATGAATGTGGCAAAAGTCATCCAACTGTCAGATGAGCTTTCAAATAAAATAGCGGCGGGCGAGGTGGTGG‑AACGGCCCGCCTCAGTCGTC‑AAAGAATTGTTGGAAAATGCG‑ATCGACGCTGACAGCACAGTCATTGAAATCTATATTGAGTAAGCAGGCCTTGCATCCATTCGGGTCTTGT‑ATAACGGCGAAGGAATGGAGAATGa                                                                                                     <  1:89600/249‑1 (MQ=255)
                                                                                                                                             gTACAAGCTGCAAAAGAAATTACATTAAAACGGGGTGATGAATGTGGCAAAAGTCATCCAACTGTCAGATGAGCTTTCAAATAAAATAGCGGCGGGCGAGGTTGTGG‑AACGGCCCGCCTCAGTCGTC‑AAAGAATTGGTGGAAAATGCG‑ATCGACGCTGACAGCACAGTCATTGAAATCGATATTGAGGAAGCAGGCCTTGCATCCATTCGGGTCTTGG‑ATAACGGCGAAGGAATGGAGAATGAAGATTGc                                                                                              <  1:164387/250‑1 (MQ=255)
                                                                                                                                               aCAAGCTGCAAAAGAAATTACATTAAAACGGGGTGATGAATGTGGCAAAAGTCATCCAACTGTCAGATGAGCTTTCAAATAAAATAGCGGCGGGCGAGGTTGTGG‑AACGGCCCGCCTCAGTCGTC‑AAAGAATTGGTGGAAAATGCG‑ATCGACGCTGACAGCACAGTCATTGAAATCGATATTGAGGAAGCAGGCCTTGCATCCATTCGGGTCTTGG‑ATAACGGCGAAGGAATGGAGAATGAAGATTGCAAg                                                                                           <  1:94345/251‑1 (MQ=255)
                                                                                                                                                    cTGCAAAAGAAATTACATTAAAACGGGGTGATGAATGTGGCAAAAGTCATCCAACTGTCAGATGAGCTTTCAAATAAAATAGCGGCGGGCGAGGGTGTGG‑AACGGCCCGCCTCCGTCGTC‑CAAGAATTGGTGGAAAATGCGGAT‑GACGCTGACAGCACAGTCCTTTAAATCGATAATGAGGGAACAGGCCTTGCATCCATTTGGGTCTTGG‑ATAACGGCGAAGGGATGGGGAAAGaaagatgcaaaggagc                                                                                      >  1:284418/1‑237 (MQ=255)
                                                                                                                                                         aaaGAAATTACATTAAAACGGGGTGATGAATGTGGCAAAAGTCATCCAACTGTCAGATGAGCTTTCAAATAAAATAGCGGCGGGCGAGGTTGTGG‑AACGGCCCCCCCCAGACGTC‑CAAGAAATTGCGGAAAATTCG‑ATCGACGCTAAAAGCACAGAAATTCAAATCGATATTGCGGAAGCCGGCCTTGCATCCATTCGGGGCTTGG‑CTAAACGCGGAGGAATAGAGAATGAAGATTGCAAGCGAGCTTTcc                                                                                 >  1:277199/1‑251 (MQ=255)
                                                                                                                                                          aaGAAATTACATTAAAACGGGGTGATGAATGTGGCAAAAGCCATCCAACTGTCAGATGAGCTTTCAAATAAAATAGCGGCGGGCGAGGTTGTGG‑AACGGCCCGCCTCAGTCGTC‑AAAGAATTGGTGGAAAATGCG‑ATCGACGCTGACAGCACAGTCATTGAAATCGATATTGAGGAAGCAGGCCTTGCATCCATTCGGGTCTTGG‑ATAACGGCGAAGGAATGGAGAATGAAGATTGCAAGCGTGCTTTCC‑g                                                                               <  1:199125/251‑1 (MQ=255)
                                                                                                                                                             aaaTTACATTAAAACGGGGTGATGAATGTGGCAAAAGTCATCCAACTGTCAGATGAGCTTTCAAATAAAATAGCGGCGGGCGAGGGTGTGG‑AACGGCCCGCCCCAGTCGTC‑CAAGAATTGGGGGAAAAAGGG‑ATCGACCCTTGCAGCACAGTCATTTAAATCGATATTGAGGAAGCAGGCCCTGCATCCATTCGGGGCTTGG‑ATAAAGGCGAAGGAATGGGGAATGAAGATTGCaaacaaggtttacggcg                                                                             >  1:244201/1‑236 (MQ=255)
                                                                                                                                                                ttACATTAAAACGGGGTGATGAATGTGGCAAAAGTCATCCAACTGTCAGATGAGCTTTCAAATAAAATAGCGGCGGGCGGGGGTGTGG‑GACGGCCCGCCCCAGTCGTC‑AAAGAATTGGTGGGAAAAGCG‑ATCGACGCTGACAGCACAGTCATTGAAATCGATATTGAGGGAGCAGGCCTTGCATCCAt                                                                                                                                          >  1:293889/1‑188 (MQ=255)
                                                                                                                                                                ttACATTAAAACGGGGTGATGAATGTGGCAAAAGTCATCCAACTGTCAGATGAGCTTTCAAATAAAATAGCGGCGGGCGAGGTTGTGG‑AACGGCCCGCCTCAGTCGTC‑AAAGAATTGGTGGAAAATGCG‑ATCGACGCTGACAGCACAGTCATTGAAATCGATATTGAGGAAGCAGGCCTTGCATCCAt                                                                                                                                          <  2:293889/188‑1 (MQ=255)
                                                                                                                                                                 tACATTAAAACGGGGTGATGAATGTGGCAAAAGTCATCCAACTGTCAGATGAGCTTTCAAATAAAATAGCGGCGGGCGGGGGTGGGG‑GACGGGCCCCCTCAGTCGTCCAAAAAATGGGGGGAAAA‑GCG‑GTCGACGCTGGCAGCACAGTCCTTGAAAACGATATTGGGGAAACAAGCCCTTCATCCAATCGGGGCTTGG‑ATAACGGCGGAGGAATGGGGGAAGAAAATTGCAAGCGAGCTTTCC‑GGCGCCAc                                                                        >  1:79549/1‑251 (MQ=255)
                                                                                                                                                                  aCATTAAAACGGGGTGATGAATGTGGCAAAAGTCATCCAACTGTCAGATGAGCTTTCAAATAAAATAGCGGCGGGGGAGGTTGTGG‑ACCGGCCCGCCTCAGTCGTC‑AAAGAATTGGTGGAAAATGCG‑ATCGACGCTGACAGCACAGTCATTGAAATCGATATTGAGGAAGCAGGCCTTGCATCCATTCGGGTCTTGG‑ATAACGGCGAAGGAATGGAGAATGAAGATTGCAAGCGAGCTTTCC‑GGCGCCAc                                                                        <  2:73178/250‑1 (MQ=255)
                                                                                                                                                                   aaatAAAACGGGGTGAGGAATGTGGAAAAAGTCATCCAACTGTCAGATGAGCTTTCAAATAAAATAGCGGCGGGCGAGGTTGTGG‑ACCGGCCCGCCTCAGTCGTC‑AAAGAATTGGTGGAAAATGCG‑ATCGACGCTGACAGCACAGTCATTGAAATCGATATTGAGGAAGCAGGCCTTGCATCCATTCGGGTCTTGG‑ATAACGGCGAAGGAATGGAGAATGAAGATTGCAAGCGAGCTTTCC‑GGCGCCACGc                                                                      <  2:31093/248‑1 (MQ=255)
                                                                                                                                                                      tAAAACGGGGTGATGAATGTGGCAAAAGTCATCCAACTGTCAGATGAGCTTTCAAAGAAAATAGCGGCGGGCGAGGTTGTGG‑AACGGCCCGCCTCAGTCGTC‑AAAGAATTGGTGGAAAATGCG‑ATCGACGCTGACAGCACAGTCATTGAAATCGATATTGAGGAAGCAGGCCTTGCATCCATTCGGGTCTTGG‑ATAACGGCGAAGGAATGGAGAATGAAGATTGCAAGCGAGCTTTCC‑GGCGCCACGCaa                                                                    <  2:88404/250‑1 (MQ=255)
                                                                                                                                                                        aatcggggttgtgaAGGTGGCAAAAAGCAACAAACTGTCCGCTGAGCTTTCAAATAAAATAGCGGCGGCCGAGGTTGTGG‑AGCGGCCCGCCTCAGTCGTC‑AAAGAATTGGTGGAAAATGCG‑ATCGACGCTGACAGCACAGTCATTGAAAGCGATATGGAGGAAGCAGGCCTTGCATCCATTCGGGTCTTGG‑ATAACGGCGAAGGAATGGAGAATGAAGATTGCAAGCGAGCTTTCC‑GGCGCCACGCAACGa                                                                 <  2:181039/240‑1 (MQ=255)
                                                                                                                                                                           cGGGGTGATGAATGTGGCAAAAGTCATCCAACTGTCAGATGAGCTTTCAAATAAAATAGCGGCGGGCGGGGGTGTGG‑GACGGCCCGCCCCAGTCGTC‑AAAGAATTGGTGGGAAAAGCG‑GGAGACGCTGGCAGCACAGTCATTGAAATCGATATTGAGGAAGCAGGCCTTGCATCCCTTCTGGTCTTGG‑GTAACGGCGGAGGAATGGAGAATGAAGATTTGAAGGGAGCTTTCC‑GGCGCCACGCAACGAgg                                                               >  1:160816/1‑249 (MQ=255)
                                                                                                                                                                                  atgaATGTGGCAAAAGTCATCCAACTGTCAGATGAGCTTTCAAATAAAACAGCGGCGGGCGAGGGTGTGG‑AACGGCCCGCCTCACTCGTC‑AAAGAATTGGTGGAAAAAGGG‑ATCGACGCGGACAGCACAGTCATTGAAATCGATATTGAGGAAGCAGGCCTTGCATCCATTCGGGGCTTGG‑AGAACGGCGGAGGAATGGAGAATGAAGATTGCAAGCGGGGTTTCC‑GGGGcccaaccacggggaaaataa                                                        >  2:120805/1‑232 (MQ=255)
                                                                                                                                                                                          ggCAAAAGACATCCAACTGCCAGATGAGCTTTAAAATAAAATAGGGGCGGGCGAGGTTGTGG‑AACGGCCCGCCTCAGTCGTC‑AAAGAATTGGTGGAAAATGCG‑ATCGACGCTGACAGCACAGTCATTGAAATCGATATTGAGGAAGCAGGCCTTGCATCCATTCGGGTCTTGG‑ATAACGGCGAAGGAATGGAGAATGAAGATTGCAAGCGAGCTTTCC‑GGCGCCACGCAACGAGTAAAATAAAAGATGaaa                                               <  2:305322/251‑1 (MQ=255)
                                                                                                                                                                                                            ggcAGATGAGCTTTCAAATAAAATAGCGGCGGGCGAGGTTGTGG‑AACGGCCCGCCTCAGTCGTC‑AAAGAATTGGTGGAAAATGCG‑ATCGACGCTGACAGCACAGTCATTGAAATCGATATTGAGGAAGCAGGCCTTGCATCCATTCGGGTCTTGG‑ATAACGGCGAAGGAATGGAGAATGAAGATTGCAAGCGAGCTTTCC‑GGCGCCACGCAACGAGTAAAATAAAAGATGAAAATGATTTATTCAGGGTg                              <  2:261817/248‑1 (MQ=255)
                                                                                                                                                                                                                   gAGCTTTCAAATAAAATAGCGGCGGGCGAGGGTGTGG‑GACCGGCCGCCTCCGACGTC‑AAAGAATTGGTGGGAAAAGCG‑GTAGCCGCGGACAGCACAATCCTTGAAAACGGTATTTAAGAAGAAGGCCTTGCATCCATTAGGGGCTTGGGTTAACG‑CGGAGGAGAGGGGGATGAAGAGTGTAAGCGAGGTTTCCCGGGGCC‑CGCCACGAGTAAAAAAAAAAAATAAAAAGATTTATTTAGGGTGAGAACGcc                      >  2:74578/1‑250 (MQ=255)
                                                                                                                                                                                                                        ttCAAATAAAATAGCGGCGGGCGAGGTTGTGG‑AACGGCCCGCCTCAGTCGTC‑AAAGAATTGGTGGAAAATGCG‑ATCGACGCTGACAGCACAGTCATTGAAATCGATATTGAGGAAGCAGGCCTTGCATCCATTCGGGTCTTGG‑ATAACGGCGAAGGAATGGAGAATGAAGATTGCAAGCGAGCTTTCC‑GGCGCCACGCAACGAGTAAAATAAAAGATGAAAATGATTTATTCAGGGTGAGAACGCTCGGt                  <  2:119085/250‑1 (MQ=255)
                                                                                                                                                                                                                                       ggcgtacgAGGTTGTGG‑AACGGCCCGCCTCAAGCGTC‑AAAGAAGTGGTGGAAAATGCG‑AGCGACGCTGACAGCACAGTCATTGAAATCGATATTGAGGAAGCAGGCCTTGCATCCATTCGGGGCTTGG‑ATAACGGCGAAGGAATGGAGAATTAAGATTGCAAGCTAGCTTTCC‑GGCGCCACGCAACGAGTAAAATAAAAGATGAAAATGATTTATTCAGGGTGAGAGCGCTCGGTTTCAGTGGAGAAGccc  <  2:168985/245‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
GCCGGCACAGCTGTCCTTTTTTGACGAAGCGGAAAAGCCGGCTGAAACACCAAAGCTTTCAAAAAAAGAAAAGCAAGTGATCGATGCTTTCAAATCACTTAATATACTGGATATGACACCGCTTGAAGCGATGAATGAAATGTACAAGCTGCAAAAGAAATTACATTAAAACGGGGTGATGAATGTGGCAAAAGTCATCCAACTGTCAGATGAGCTTTCAAATAAAATAGCGGCGGGCGAGGTTGTGG‑AACGGCCCGCCTCAGTCGTC‑AAAGAATTGGTGGAAAATGCG‑ATCGACGCTGACAGCACAGTCATTGAAATCGATATTGAGGAAGCAGGCCTTGCATCCATTCGGGTCTTGG‑ATAACGGCGAAGGAATGGAGAATGAAGATTGCAAGCGAGCTTTCC‑GGCGCCACGCAACGAGTAAAATAAAAGATGAAAATGATTTATTCAGGGTGAGAACGCTCGGTTTCAGGGGAGAAGCCC  >  NC_000964/1778153‑1778634

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