Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 762,113 T→C 23.7% S198P (TCC→CCC)  sucB → dihydrolipoyltranssuccinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913762,1130TC23.7% 112.8 / 7.0 65S198P (TCC→CCC) sucBdihydrolipoyltranssuccinase
Reads supporting (aligned to +/- strand):  ref base T (34/14);  new base C (3/14);  total (37/28)
Fisher's exact test for biased strand distribution p-value = 1.89e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCG‑AAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGC‑AACCGGCTCTGGCTGCACGTAGTGAAAAA‑CGTGTCCCGATGACT‑CGCCTGCGTAAGCGTGTGGCAGAGCG‑TCTGCTGGAAGCGAAAAA‑CTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGACGGCGATGACGTGGTTTACCACAACTATTTCGACGTCAGCATGGCGGT  >  NC_000913/761879‑762339
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                  
cGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCG‑AAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGC‑AACCGGCTCTGGCTGCACGTAGTGAAAAA‑CGTGTCCCGATGACT‑CGCCTGCGTAAGCGTGTGGCAGAGCG‑TCTGCTGGAAGCGAAAAA‑CTCCa                                                                                                                                                                                                                                 >  1:255286/1‑238 (MQ=255)
     gctggctgAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCG‑AAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGC‑AACCGGCTCTGGCTGCACGTAGTGAAAAA‑CGGGTCCCGATGACT‑CGCCTGCGTAAGCGTGTGGCAAAGGG‑GCTGCTGGAAGCGAAAAA‑CCCCACCGCCa                                                                                                                                                                                                                           >  1:791606/1‑239 (MQ=255)
                    tctcGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCG‑AAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGC‑AACCGGCTCTGGCTGCACGTAGTGAAAAA‑CGTGTCCCGATGACT‑CGCCTGCGTAAGCGTGTGGCAACGCG‑TCGGCTGGAAGCGAAAAA‑CCCCACCGCCATGCTGACCACGTTCa                                                                                                                                                                                                            >  1:785036/1‑239 (MQ=255)
                        gagccgggcgcCATGAAAGGCAGGGGGGGGGGTGGTCGGCGGACGGGTAAGGTGGGGGAAAAACGTGGGGGGAAAGCCCCGGGG‑AAAAGGTCTGCTCGGGTGCGGGCTGCCCCGGCGGCGC‑AACGGGCGCTGGCTGCAAGTAGTAAAAAA‑CGGGTCCGGTGGACG‑CCCCTGCGTAAGCGGGTGGCAGAGAG‑TCTGCGGGAAGGAAAAAC‑CTCCACCGCCATGCTGACCACGTTCAACGa                                                                                                                                                                                                        <  2:578277/232‑1 (MQ=255)
                         aCGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCG‑AAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGC‑AACCGGCTCTGGCTGCACGTAGTGAAAAA‑CGGGTCCCGATGACC‑CGCCTGCGTAACCGGGTGGCAGAGCG‑GCGGCTGGAAGCGAAAAA‑ATcccccccctaccgcaacccctccaaccaa                                                                                                                                                                                                       >  2:632676/1‑218 (MQ=255)
                              aGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCG‑AAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGC‑AACCGGCTCTGGCTGCACGTAGTGAAAAA‑CGTGTCCCGATGACT‑CGCCTGCGTAAGCGTGTGGCAGAGCG‑TCTGCTGGAAGCGAAAAA‑CTCCACCGCCATGCTGACCACGTTAAACGAAGTCaa                                                                                                                                                                                                  >  1:663723/1‑239 (MQ=255)
                                     ttAAAGGCCCGGGGGGGGGGGGGGGTCTGACTCGGGAAGGTGTGGAAAACCAGCTGGGGAAACCCCTGGGG‑AAAGGGGCTGCTCCGGCAGGGGCTGC‑CCGGGGGGGCAAACCGGCTCTGGCTGACGGGAGGGAAAAA‑GGTGTCCCGATGACT‑CGCCTGGTTAAGCGTGTGGCAGAGCG‑TCTGCTGGAAGCGAAAAA‑CTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAg                                                                                                                                                                                           <  2:279600/239‑1 (MQ=255)
                                         agcccccgggtggggtggTCCTCTGACGCGGGAAGGTGTGGAAAAAAATCGGGCGAAAGCCCCGGCG‑AAAGGGTCTGCTCGGGCAGGGGCTGCTCCGGCGGCGC‑ACCCGGCTCTGGCGGCACGTAGTGAAAAA‑CGTGTCCCGATGACT‑CGCCTGCGTAAGCGTGTGGCAGAGCG‑TCTGCGGGAAGCGAAAAC‑CCCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGa                                                                                                                                                                                       <  2:1069343/227‑1 (MQ=255)
                                             aCCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCG‑AAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGC‑AACCGGCCCTGGCTGCACGTAGTGAAAAA‑CGTGTCCCGATGACT‑CGCCTGCGGAAGCGTGGGGCAGAGCG‑TCTGCTGGAAGCGAAAAA‑CTCCACCGCCATGCGGACCACGTTCAACGAAGTCAAAAAAAAGCCGAtttg                                                                                                                                                                                   >  2:1030460/1‑237 (MQ=255)
                                               cGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCG‑AAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGC‑AACCGGCTCTGGCTGCACGTAGTGAAAAA‑CGTGTCCCGATGACT‑CGCCCGCGTACACGGGTGGCAGAGCG‑TCCGCTGGAAGCGAAAAA‑CCCCACCGCCATGCTGACCAACTCTAACAAAGTCAACAAAAAACCGATTAAgg                                                                                                                                                                                 >  2:186625/1‑239 (MQ=255)
                                                ggTGTGGGGGGTCGTTTGACCCGGGAAGATGGGGAAAAACATCGGGCGAAAGCCCGGGCG‑AAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGC‑AACCGGCTCTGGCTGCACGTAGTGAAAAA‑CGTGTCCCGATGACT‑CGCCTGCGTAAGCGTGTGGCAGAGCG‑TCTGCTGGAAGCGAAAAA‑CTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGa                                                                                                                                                                                <  1:557240/239‑1 (MQ=255)
                                                    tgggtggTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCG‑AAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGC‑AACCGGCTCTGGCTGCACGTAGTGAAAAA‑CGTGTCCCGATGACT‑CGCCTGCGTAAGCGTGGGGCAGAGCG‑TCTGCTGGAAGCGAAAAA‑CTCCACCGCCATGCTGACCACGTTCAACGAAGTCAAAAGGAAGCCGATTTAGGAACTg                                                                                                                                                                            >  2:672196/1‑239 (MQ=255)
                                                    tgggtggTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCG‑AAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGC‑AACCGGCTCTGGCTGCACGTAGTGAAAAA‑CGCGTCCCGATGACT‑CGCCTGCGTAAGCGTGTGGCAGAGCG‑TCTGCTGGAAGCGAAAAA‑CTCCACCGCCATGCTGCCCACGTTCAACGAAGTCAACATGAAGCCGATTATGGAttt                                                                                                                                                                             >  1:892543/1‑236 (MQ=255)
                                                           tcgtcTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCG‑AAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGC‑AACCGGCTCTGGCTGCACGTAGTGAAAAA‑CGTGTCCCGATGACT‑CGCCTGCGTAAGCGTGTGGCCGAGGG‑TTCGCTGGAAGCGCAAAA‑CTCCACCGCCATCCTGACCACGTTCAACGAAGTCAACATGAACCCCACTTTGGCTCTCCGaaaac                                                                                                                                                                     >  2:987396/1‑237 (MQ=255)
                                                           ccgtcTGACTCGTGAAGATGTGGAAAACCATCGGGCGAAAGCCCCGGCG‑AAAGAGTCTGCCCCGGCAGCGGCTGCTCCGGCGGCGC‑AACCGGCTCTGGCTGCACGTAGTGAAAAA‑CGTGTCCCGATGACT‑CGCCTGCGTAAGCGTGTGGCAGAGCG‑TCTGCTGGAAGCGAAAAA‑CCCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAg                                                                                                                                                                      <  1:591202/237‑1 (MQ=255)
                                                            cgtcTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCG‑AAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGC‑AACCGGCTCTGGCTGCACGTAGTGAAAAA‑CGTGTCCCGATGACT‑CGCCTGCGTAAGCGTGTGGCAGAGCG‑TCTGCTGGAAGCGAAAAA‑CTCCACCGCCATGCTGCCCACGTTCAACGCAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCa                                                                                                                                                                    >  1:919220/1‑239 (MQ=255)
                                                            cgtcTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCG‑AAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGC‑AACCGGCTCTGGCTGCACGTAGTGAAAAA‑CGTGTCCCGATGACT‑CGCCTGCGTAAGCGTGGGGCAGAGCG‑TCTGCTGGAAGCGAAAAA‑CTCCCCCGCCATGCTGCCCCCGTTCACAGAAGTCAACATGAAACCGATAATGGAATTCGGTAAGCa                                                                                                                                                                    >  2:840842/1‑239 (MQ=255)
                                                                       tGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCG‑AAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGC‑AACCGGCTCTGGCTGCACGTAGTGAAAAA‑CGTGTCCCGATGACT‑CGCCTGCGTAAGCGTGTGGCAGAGCG‑TCTTCTGGAAGCGAAAAA‑CTCCACCGCCATGCTGACCACGTTAAACGAAGTCACAATGAAGCCGATTAAGGATCTGCGTAAACAGAAAGGGGAAg                                                                                                                                                         >  2:65344/1‑239 (MQ=255)
                                                                              gtgGAAAAACATCTGGCGAAAGCCCCGGCG‑AAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGC‑AACCGGCTCTGGCTGCACGTAGTGAAAAA‑CGTGTCCCGATGACT‑CGCCTGCGTAAGCGTGTGGCAGAGCG‑TCTGCTGGAAGCGAAAAA‑CTCCACCGCCATGCTGAc                                                                                                                                                                                                                    >  2:336693/1‑173 (MQ=255)
                                                                              ggggAAAAACATCTGGCGAAAGCCCCGGCG‑AAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGC‑AACCGGCTCTGGCTGCACGTAGTGAAAAA‑CGTGTCCCGATGACT‑CGCCTGCGTAAGCGTGTGGCAGAGCG‑TCTGCTGGAAGCGAAAAA‑CTCCACCGCCATGCTGAc                                                                                                                                                                                                                    <  1:336693/171‑1 (MQ=255)
                                                                                      acctCTGGGGAAAACCCCGGGGAAAAGGGTCTCCCCGGGCAGGGGCTGCCCGGGGGGCGC‑AACCGGCTTTGGCTGCACGTAGTGAAAAA‑CGTGTCCCGATGACT‑CGCCTGCGTAAGCGTGTGGCAGGGCG‑TCTGCTGGAAGCGAAAAA‑CTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGc                                                                                                                                           <  2:255286/236‑1 (MQ=255)
                                                                                       cctggggagAAACCCCGGGGG‑AAGGGGCCTGCCCGGGCGGGGGTGGCGCCGGGGGCCA‑ACCCGGCTGTGGTGGCAGGGAGGGAAAAC‑CGTGTGCCGATGACC‑CCCCTGCGAAGGGGGGGGGCAGAGTG‑TCTGCGGGAAGGGAAAAA‑CCCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCa                                                                                                                                         <  2:148130/231‑1 (MQ=255)
                                                                                            ggCGAAAGCCCCGGCG‑AAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGC‑AACCGGCTCTGGCTGCACGTAGTGAAAAA‑CGGGTCCCGATGACT‑CGCCTGCGTAAGCGTGTGGCAGAGCG‑TCTGCTGGAAGCGAAAAA‑CTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACCTGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGCAAAAAGCCACGGc                                                                                                                                     >  1:643033/1‑238 (MQ=255)
                                                                                                  aGCCCCGGCG‑AAAGCGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGC‑AACCGGCTCTGGCTGCACGTAGTGAAAAA‑CGTGTCCCGATGACT‑CGCCTGCGTACGCGTGGGGCAGAGCG‑CCTGCTGGAAGGGAAAAA‑CTCCCCCCCCATGCTGGCCCCGTTCACCGAAGTCAACATAAAGCCGACTATGGACTTGCGTACACAGAACGGTAAAACGTTTGAAAAACCCCAAGGCATCCGTc                                                                                                                              >  2:62515/1‑239 (MQ=255)
                                                                                                       cGGCG‑AAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGC‑AACCGGCTCTGGCTGCACGTAGTGAAAAA‑CGTGTCCCGATGACT‑CGCCTGCGTAAGCGTGTGGCAGAGCG‑TCTGCTGGAAGCGAAAAA‑CTCCCCCGCCATGCTGACCCCGTTCAACGAAGCCAACATGAAGCCGATTAAGGATCTGCGTAAGCAATACGGGAAAGCGTTTGAAAAAACCCACGGACACCGCCTgggg                                                                                                                         >  2:482241/1‑238 (MQ=255)
                                                                                                             aagtggtttGCCCCGGGAGGGGTTGCTCGGGGGGCGC‑ACCCGGGTCGGGCGGCGCGTGGTGAAAAC‑CGGGTCCGGATGACT‑CGCCTGCGTAAGGGTGGGGCAGGTCG‑TCTGCTGGAAGCAAAAAA‑CTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTAt                                                                                                                    <  2:1044712/231‑1 (MQ=255)
                                                                                                              aagggttcgctcgggcgggggtggctCCGGGGGCAC‑AACGGGCCTTGGTGGCCCGTGGTGAAAAA‑CGTGTCCCGATGCCT‑CGCCTGGGTAGGGGTGGGGCAGGGCG‑TCTGTGGAAAGAAAAAAA‑CCCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATg                                                                                                                   <  2:83806/216‑1 (MQ=255)
                                                                                                               tggtgttGCGCCGGCGGGGGGTCCGCCGGGTGCGC‑ACCCGGGTCTGGCGACAGGGTGGGAAAAA‑CGTGTCCGGATAACT‑CGCCTGCGTAGGCTTGGGGCAGGGCT‑TCGGGTGGAAGCGAAAAA‑CCCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGt                                                                                                                  <  2:405339/233‑1 (MQ=255)
                                                                                                                     tGCTCCGGCAGCGGCTGCTCCGGCGGCGC‑AACCGGCTCTGGCTGCACGTAGTGAAAAA‑CGTGTCCCGATGACT‑CGCCTGCGTAAGCGTGTGGCAGAGCG‑TCTGCTGGAAGCGAAAAA‑CTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACCTGAAGCCGATTATGGTTCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCCCGGCAACCGCCTGGGCTTTATGTCCTTc                                                                                                             >  1:899380/1‑238 (MQ=255)
                                                                                                                            gCAGCGGCTGCTCCGGCGGCGC‑AACCGGCTCTGGCTGCACGTAGTGAAAAA‑CGTGTCCCGATGACT‑CGCCTGCGTAAGCGTGTGGCAGAGCG‑TCTGCTGGAAGCGAAAAA‑CTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGt                                                                                                                               >  1:13009/1‑213 (MQ=255)
                                                                                                                                 ggCTGCTCCGGGGGCGC‑ACCCGGCTCTGGCTGCACGTAGTGAAAAA‑CGTGTCCCGATGACT‑CGCTTGCGTAAGCGTGTGGCAGAGCT‑TCTGCTGGAAGCGAAAAA‑CTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCg                                                                                                 <  1:632676/238‑1 (MQ=255)
                                                                                                                                           gcggcgC‑AACCGGCTCTGGCTGCACGTAGTGAAAAA‑CGTGTCCCGATGACT‑CGCCTGCGTAAGCGTGTGGCAGAGCG‑TCTGCTGGAAGCGACAAA‑CTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAACGTTTGAAAAACGCCCCGCCATCCGTCCGGGGTTCAAGCCCTTCAACCGGGAAGAGGGGGTTGGAg                                                                                       >  2:575825/1‑238 (MQ=255)
                                                                                                                                                 c‑AACCGGCTCTGGCGGCCCGTAGTGAAAAA‑CGTTTCCCGATGTCC‑CGCCTGCGTAAGCGTTTGGCAGGGCG‑TCTGCTGGAAGCGAAAAA‑CCCCCCCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGa                                                                                 <  1:220150/238‑1 (MQ=255)
                                                                                                                                                     cgggtccggtcggatgggaGTAAAAAC‑CGTGTCCCGTTGACC‑GCCCTGGTTAGGGGGGGGGCAGAGCG‑GTTGCGGAAAGCGAAAAC‑TCCCACCGCCATGCTGACCACGTTCAACGAAGGCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACg                                                                             <  2:643033/221‑1 (MQ=255)
                                                                                                                                                     ccGGCTTTGGCTGCACGTAGTGAAAAA‑CGTGTCCAGATGACT‑CGCCTGCGTAAGCGTGGGGCAGAGCG‑TCGGCTGGAAGCGAAAAA‑CCCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACg                                                                             <  1:65344/239‑1 (MQ=255)
                                                                                                                                                          gctggcgtgccGTAGTGAAAAA‑CGTGTCCCGAGGCCC‑CGCCTGCGTAGGCGTGTGGCAGAGCG‑TCTGCTGGAAGCGAAAAA‑CTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTAcc                                                                        <  1:305315/229‑1 (MQ=255)
                                                                                                                                                           cgtgCTCCACCTGGTGAAAAA‑CGTGTCTCCGTGACTCCGCCGCCGTAGGCGGGGTGAAGTGCT‑CTTGCTGAAACCAAAAAC‑CCCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTAcc                                                                        <  2:31711/236‑1 (MQ=255)
                                                                                                                                                                ggCACGTA‑TGAAAAAGCGTGCCCCGAGGACT‑CGCCTGCGTAAGCGTGTGGCAGAGCG‑TCTGCTGGAAGCGAAAAA‑CTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGaa                                                                   <  1:184709/237‑1 (MQ=255)
                                                                                                                                                                   aCGTAGTGAAAAA‑CGTGTCCCGATGACT‑CGCCTGCGTAAGCGTGTGGCAGAGCG‑TCTGCTGGAAGCGAAAAA‑CTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCCTCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGGGGTTGAAGCCCCGAAACGTTACCCGGaaagga                                                               >  2:66233/1‑235 (MQ=255)
                                                                                                                                                                   aCGTAGTGAAAAA‑CGTGTCCCGATGACT‑CGCCTGCGTAAGCGTGTGGCAGAGCG‑TCTGCTGGAAGCGAAAAA‑CTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAACGGTGGTTGAAGCCCTGAAActgcctcttataaac                                                                >  1:316263/1‑224 (MQ=255)
                                                                                                                                                                    acggggaaaaag‑CTTGTCCCGTTGTCT‑CGCCTGTGTACGCGGGGGGACGGGCG‑TCTGCTTGAACCAAAAAC‑CCCGCCCGCCATGCTGACCACGTTCAACGAAGGCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGaa                                                              <  2:338410/232‑1 (MQ=255)
                                                                                                                                                                      tAGTGAAAAA‑CGTGTCCCGATGACT‑CGCCTGCGTAAGCGTGTGGCAGAGCG‑TCTGCTGGAAGCGAAAAA‑CTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCATCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAACCGTTACCCGCAAGTGAAc                                                             >  1:58803/1‑238 (MQ=255)
                                                                                                                                                                       gagaaaaac‑CGGGTCCCGATGCCT‑CGCCTGCGTAAGCGTGTGGCAGAGCG‑TCTGCTGGAAGCGAAAAA‑CCCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGaa                                                              <  1:672196/232‑1 (MQ=255)
                                                                                                                                                                        gTGAAAAA‑CGTGTCCCGATGACT‑CGCCTGCGTAAGCGTGTGGCAGAGCG‑TCTGCTGGAAGCGAAAAA‑CTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGGGGTTGAAGCCCGGAAACGTTACCCGGAAGTGAACg                                                            >  2:532025/1‑237 (MQ=255)
                                                                                                                                                                          aaggagaccggcccagATGACT‑CGCTGGGGGAAGCGTGTGAGAAAGCG‑TCGGGGGGAA‑CGAAAAACCTCCGCCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCtt                                                         <  1:482947/224‑1 (MQ=255)
                                                                                                                                                                                 cGTGCCCCGATGACT‑CGCCTGCGTAAGCGTGTGGCAGAGCG‑TCTGCTGGAAGCGAAAAA‑CCCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGa                                                  <  1:972189/239‑1 (MQ=255)
                                                                                                                                                                                    gtCCCGATGACT‑CGCCTGCGTAAGCGTGTGGCAGAGCG‑TCTGCTGGAAGCGAAAAA‑CTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGaa                                                              >  1:302645/1‑224 (MQ=255)
                                                                                                                                                                                    gtCCCGATGACT‑CGCCTGCGTAAGCGTGTGGCAGAGCG‑TCTGCTGGAAGCGAAAAA‑CTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGGTTGTATc                                                    >  1:338410/1‑234 (MQ=255)
                                                                                                                                                                                    gtCCCGATGACT‑CGCCTGCGTAAGCGTGTGGCAGAGCG‑TCTGCTGGAAGCGAAAAA‑CTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGGGGTTGAAGCCCGGAAACGTTACCCGGAAGGGAACCCTCCTATCGACgg                                               >  2:131859/1‑239 (MQ=255)
                                                                                                                                                                                    gtCCCGATGACT‑CGCCTGCGTAAGCGTGTGGCAGAGCG‑TCTGCTGGAAGCGAAAAA‑CTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGACGGGGCTTTATGTACTTCTACTGGAAAGCGGGGGGTGCAGCCCTGAAACGTTACCCGAAAGGGaaagactcaataaaagg                                               >  2:203218/1‑224 (MQ=255)
                                                                                                                                                                                    gtCCCGATGACT‑CGCCTGCGTAAGCGTGTGGCAGAGCG‑TCTGCTGGAAGCGAAAAA‑CTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATGATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGCCTGGGTTTTATGTCCTCCTACGTGAAACGGGGGGTTGAAACCCTAAAACGTACCCCGGaacgaaacgatgaataagaggg                                               >  2:458988/1‑219 (MQ=255)
                                                                                                                                                                                            gACT‑CGCCGGGGGAAGCGTGTGGAAG‑GCGTTCTGCTGGAAGCGAAAAA‑CTCCACCGCCATGCTGAACACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGACGGCGATGACg                                       <  1:1005395/239‑1 (MQ=255)
                                                                                                                                                                                             aCT‑CCCCTGGATAAGGGTGGGGCAGGGGG‑TGTGGTGGAAGCGAAAAC‑CTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGACGGCGATGACg                                       <  1:131859/238‑1 (MQ=255)
                                                                                                                                                                                               t‑CGCCTGCGTAAGCGTGTGGCAGAGCG‑TCTGCTGGAAGCGAAAAA‑CTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGGGGTTGAACCCCTGAAACGTTACCCGAGAGGGAACGCTTCTATCGAGGGcctggccctta                                    >  2:482947/1‑229 (MQ=255)
                                                                                                                                                                                               t‑CGCCAGCGGAGGGGTGTGGC‑GGGCGTTCTGCTGGAAGCAAAAAC‑CCCCACCGCCATGCTGACCAAGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGACGGCGATGACGTgg                                    <  1:62515/239‑1 (MQ=255)
                                                                                                                                                                                                  gCCTGCGAAAGCGTGGGGGAGGGCG‑TCTGCTGGAAGCGAAAAA‑CTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGACGGCGATGACGTGGtt                                  <  1:232959/239‑1 (MQ=255)
                                                                                                                                                                                                                 tgGCAGAGCG‑TCTGCTGGAAGCGAAAAA‑CTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCttt                                                                                                                      <  1:123041/140‑1 (MQ=255)
                                                                                                                                                                                                                 tgGCAGAGCG‑TCTGCTGGAAGCGAAAAA‑CTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCttt                                                                                                                      >  2:123041/1‑140 (MQ=255)
                                                                                                                                                                                                                 tgGCAGAGCG‑TCTGCTGGAAGCGAAAAA‑CTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGACGGCGATTACGTGGTTTACCACAACTATTTc                   >  1:148481/1‑239 (MQ=255)
                                                                                                                                                                                                                 tgGCAGAGCG‑TCTGCTGGAAGCGAAAAA‑CTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGACGGCGATGACGTGGTTTACCACAACTATTTc                   >  1:178470/1‑239 (MQ=255)
                                                                                                                                                                                                                 tgGCAGAGCG‑TCTGCTGGAAGCGAAAAA‑CTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGGGAAAGCGGGGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGGCGGGGGGTACGGGGTTTACACAAACAATTTc                   >  2:874959/1‑239 (MQ=255)
                                                                                                                                                                                                                       agCG‑TCTGCTGGAAGCGAAAAA‑CTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGACGGCGATGACTTGGTTTACCCCAACTATTTCGACGTc             >  1:715360/1‑239 (MQ=255)
                                                                                                                                                                                                                            tCTGCTGGAAGCGAAAAA‑CTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGACGGCGATGACGGGGTTTACCACAACTATTTCGACGTCAGCa         >  1:631313/1‑239 (MQ=255)
                                                                                                                                                                                                                              tgtaggAAGCAAAACA‑CCCCACCGCCATGCTGACCACGTTAAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGACGGCGATGACGTGGTTTACCACAACTATTTCGACGTCAGCATg       <  1:790920/235‑1 (MQ=255)
                                                                                                                                                                                                                                   gAAGCGAAAAA‑CTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGACGGCGATGACGTGGTTTACCACAACTATTTCGACGTCAGCATGGCggg  >  1:214603/1‑238 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                  
CGCCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTGGTCGTCTGACTCGTGAAGATGTGGAAAAACATCTGGCGAAAGCCCCGGCG‑AAAGAGTCTGCTCCGGCAGCGGCTGCTCCGGCGGCGC‑AACCGGCTCTGGCTGCACGTAGTGAAAAA‑CGTGTCCCGATGACT‑CGCCTGCGTAAGCGTGTGGCAGAGCG‑TCTGCTGGAAGCGAAAAA‑CTCCACCGCCATGCTGACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGACGGCGATGACGTGGTTTACCACAACTATTTCGACGTCAGCATGGCGGT  >  NC_000913/761879‑762339

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