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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	762,113	0	T	C	25.9%	67.4 / 2.6	41	S198P (TCC→CCC)	sucB	dihydrolipoyltranssuccinase
Reads supporting (aligned to +/- strand):  ref base T (22/7);  new base C (2/10);  total (24/17)
Fisher's exact test for biased strand distribution p-value = 1.06e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTG‑GTCGTCTGACTCGTGAAGATGTGGAAAAACA‑TCTGGCGAAAGCCCCGGCGAAAGAGTCT‑GCTCCGGCAGC‑GGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGACGGCGATGACGTGGTTTACCACAACTATTTCGACGTCAGCATGGC  >  NC_000913/761881‑762336                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                 ccTGTTGGTTGAACACATTTCGGGCTCCACGGCCTTTAAAGGCCTCGGTGTGGGTG‑GTCGCTTGACCGGGGAAGAGGTGGAAAAACA‑TCGGGCGAAAGCCCCGGCAAAAGAGTCG‑GCTCCGGCAGG‑GGCTGCTCCGGCGGCGCAACCGGCTTTGGCTGCCCGTAGTGAAAACCGTGCCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCg                                                                                                                                                                                                                            <  2:369033/239‑1 (MQ=255) ccTGCTGGCTGAACACAATCTCGACGCCAGGGCAATTAAAGGCACCGGTGTGGGTG‑GTCGTCTGACTCGTGAAGATGTGGAAAAACA‑TCTGGCGAAAGCCCCGGCGAAAGAGTCT‑GCTCCGGCAGC‑GGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCg                                                                                                                                                                                                                            <  1:590688/239‑1 (MQ=255)        gcggtacacAATCTCCGCGCCCGCGCCCTTTAAGGGCCCCGGGTGGGGGTGTCGCTTGAC‑CGTTAAAAGGTGGGAAAACATTCTGGAGAAACCCCCGG‑GAAAGGATCTGGCTCCGGCAGC‑GGTGCCCCCGGGGGCCCACCCGGTTTTGGCTGCCGGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCAAAAAC‑CCCCACCGCCATGc                                                                                                                                                                                                                      <  2:274272/234‑1 (MQ=255)                               cgcCATTAAAGGCACCGGTGTGGGTG‑GTCGTCTGACTCGTGAAGATGTGGAAAAACA‑TCTGGCGAAAGCCCCGGCGAAAGAGTCT‑GCTCCGGCAGC‑GGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACa                                                                                                                                                                                             >  1:432896/1‑239 (MQ=255)                                    ttAAAGGCCCCGGTGTGGGTG‑GGCGTCTGACTCGTGAAGATGTGGAAAAACA‑TCTGGCGAAAGCCCCGGGAAAAGAGTCT‑GCTCCGGGAGC‑GGCCGCTCCGGCGGCGCAACCGGCTCTGGCTGCCCGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAg                                                                                                                                                                                        <  1:626918/239‑1 (MQ=255)                                            aCCGGTGTGGGTG‑GTCGTCTGACTCGTGAAGATGTGGAAAAACA‑TCTGGCGAAAGCCCCGGCGAAAGAGTCT‑GCTCCGGCAGC‑GGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTAt                                                                                                                                                                                >  1:510449/1‑239 (MQ=255)                                            aCCGGTGTGGGTG‑GTCGTCTGACTCGTGAAGATGTGGAAAAACA‑TCTGGCGAAAGCCCCGGCGAAAGAGTCT‑GCTCCGGCAGC‑GGCTGCTCCGGCGACGCAACCGGCTCTGGCTGCACGCAGTGAAAACAGTGTTCCGGTAACTCGCCTGCGGAAGCGTGTGGCAGAGGGTCTGCTTGCA‑GCGCAAAA‑CCCCCCCCCCCTCTT‑GGCCACGTTCAACGAAATCAACAGAAAACCGAtatg                                                                                                                                                                                >  2:80145/1‑236 (MQ=255)                                              cGGTGTGGGTG‑GTCGTCTGACTCGTGAAGATGTGGAAAAACA‑TCTGGCGAAAGCCCCGGCGAAAGAGTCT‑GCTCCGGCAGC‑GGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCACCGAAGTCAACATGAAGCCGATAATga                                                                                                                                                                              >  1:328050/1‑238 (MQ=255)                                              cGGTGTGGGTG‑GTCGTCTGACTCGTGAAGATGTGGAAAAACA‑TCTGGCGAAAGCCCCGGCGAAAGAGTCT‑GCTCCGGCAGC‑GGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATgg                                                                                                                                                                              >  1:423124/1‑239 (MQ=255)                                                gggggggtgtgtcgtcTGAATCGGGAAGATGTGAAAAACCA‑GCTGGCGAAAGCCCCGGCGAAATGGCCT‑GTCCGGGAAGG‑GGCTGCTCCGGCGGCGCAACCGGCTCTGGCGGCGCGTAGTGAAAAACGTGTCCCGATGCCTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGa                                                                                                                                                                             <  2:575493/229‑1 (MQ=255)                                                           tcgtcTGACTCGTGAAGATGTGGAAAAACA‑TCTGGCGAAAGCCCCGGCGAAAGAGTCT‑GCTCCGGCAGC‑GGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGGGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCAGGTTCAACGAAGGCAAAAGAAAGCAGATTATGGAGTTGGGTAAga                                                                                                                                                                  >  2:869791/1‑238 (MQ=255)                                                            cgtcTGACTCGTGAAGATGTGGAAAAACA‑TCTGGCGAAAGCCCCGGCGAAAGAGTCT‑GCTCCGGCAGC‑GGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTGAGGGAAAACGGGTCCCGATGACTCGCCTGCGTAGGCGTGGGCCGGAGCGTCGTCTGGAA‑ACGAAAAA‑CCCCACCCCCCAGCT‑GACACCGTTCAACAACGCCCACCTACAGCCCATTTTGGATCTGCGTCAgaa                                                                                                                                                                 >  2:695151/1‑237 (MQ=255)                                                              tcTGACTCGTGAAGATGTGGAAAAACA‑TCTGGCGAAAGCCCCGGCGAAAGAGTCT‑GCTCCGGCAGC‑GGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAga                                                                                                                                                               >  1:25832/1‑238 (MQ=255)                                                              cttcaCTCGGGAAGTTGTGAAAAACCC‑TCTGGCGAAACCCCCGGGGAAGGAGCCT‑GCCCCGGAGGC‑GGCTCCCCCGGGCGCGCACCCGGCTCTGGCTGCACGTAGGTAAAAACGTGTCCGGAGGCCCCGCCTCCGTAAGCTGGTGGCAGAGCCTCTGCTGGAA‑GGGAAAAA‑CCCAACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGt                                                                                                                                                               <  2:425605/235‑1 (MQ=255)                                                                              gtgGAAAAACA‑TCTGGCGAAAGCCCCGGCGAAAGAGTCT‑GCTCCGGCAGC‑GGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGa                                                                                                                                               >  1:620875/1‑239 (MQ=255)                                                                                ggAAAAACA‑TCTGGCGAAAGCCCCGGCGAAAGAGTCT‑GCTCCGGCAGC‑GGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAAAATGAAGCCGATTATGGATCTGCGTAAGAAGTACGGGGAAAGGTTTGaaa                                                                                                                                             >  2:519536/1‑239 (MQ=255)                                                                                ggAAAAACA‑TCTGGCGAAAGCCCCGGCGAAAGAGTCT‑GCTCCGGCAGC‑GGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGCCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTTTGGATCTGCGAAAGCGTTACGTTGAAGCGTTTGaaa                                                                                                                                             >  1:892852/1‑239 (MQ=255)                                                                                          tCTGGCGAAAGCCCCGGCGAAAGAGTCT‑GCTCCGGCAGC‑GGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGGGAAGCGTTTGcaaaaacccccg                                                                                                                                    >  1:298971/1‑232 (MQ=255)                                                                                                         gggaaaaGGGTCT‑GCTCCGGCAGC‑GGCTGCTCCGGCGGCGCAACCGGCTCGGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGGGCGTCTGCTGGAA‑GCGAAAAA‑CCCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCt                                                                                                                     <  1:148783/235‑1 (MQ=255)                                                                                                                   gtgttgcccGGGAGCGGGCTGCCCCGGCGGCCCAACGGGCTTGGGCGGCCGGTAGGGAAAAACGTGTCCCGATGACCCGCCTGGGTAAGCGTGGGGAAGGGCGTGTGCTGGAA‑GCAAAAAC‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGt                                                                                                               <  2:510449/230‑1 (MQ=255)                                                                                                                                    gctgctCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCCCGTTCAACGAAGGCAACATGAACCCGATTATGGATCCGGGGAAGCAGTACGGTGAAGCATTTGAAAAAAGCAACGGCACCCGTCTGGGCGTTATGTCCTTCTACGGGAAAGCggg                                                                                            >  2:404031/1‑238 (MQ=255)                                                                                                                                                  gcAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGGGGTTGAAGCCCTGa                                                                              >  1:335790/1‑239 (MQ=255)                                                                                                                                                                cTGAACGAAGTAAAAATCGGCTCCAGATGACTCGCCTGCGTAAGCGTGGGGCAGAGCGTCTGCTGAAA‑GCGAAAAA‑CCCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGa                                                                 <  1:695151/238‑1 (MQ=255)                                                                                                                                                                                    gtCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACACGGAAGTGAACGCTTCTATCGACg                                             >  1:88050/1‑238 (MQ=255)                                                                                                                                                                                    gtCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCCGAAACGTTACCCGGAAGTGAACGCTTCTATCGACgg                                            >  1:468572/1‑239 (MQ=255)                                                                                                                                                                                    gtCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGGGGGTGAAGCCCTGAAACGTTACCCGGAAGGGAACGCTTCTATCGACgg                                            >  2:617982/1‑239 (MQ=255)                                                                                                                                                                                             ccTCGCCGGCGTAAGCGTGGGGCAGGGCGTCTGTGGGAA‑GCGAAAAC‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGACGGCGATGACg                                    <  1:519536/237‑1 (MQ=255)                                                                                                                                                                                             aCTCGCCTGCTTAAGCTGGGGGCAGGGCGTCGGCTGAAA‑GCGAAAAC‑CCCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGACGGCGATGACg                                    <  1:80145/238‑1 (MQ=255)                                                                                                                                                                                                    ttgtgatgCGTTTGGCAGGGCTTTTGTTGG‑A‑GAGAAAAACCCCCACCTCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGACGGCGATGACGTGGTTTAc                            <  2:328050/232‑1 (MQ=255)                                                                                                                                                                                                         aggCGTGTGGCATAGCGTTGGCTGGAA‑GCAAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTAAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGACGGCGATGACGTGGTTTACcaca                        <  1:879023/236‑1 (MQ=255)                                                                                                                                                                                                         aaGCGTGGGGCAAAGCGTTCGGGGGAA‑GCGAAAAA‑CCCAGCCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGACGGCGATGACGTGGTTTACCACaa                       <  1:678520/239‑1 (MQ=255)                                                                                                                                                                                                            cGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGGTCAACGAAGTCAACATGAAGCCGATTATGGATATGCGTAAGCAGAACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCGTTATGTCCTTCTACGCGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGACGGATTCTATCGACGGCGATGACGTGGGTTACCACAACTa                    >  1:831781/1‑239 (MQ=255)                                                                                                                                                                                                                tgGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCTGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGACGGCGCAGACGGGGTTTACAACAAATATTTc                >  2:368308/1‑239 (MQ=255)                                                                                                                                                                                                                tgGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTAGGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGACGGCGATGACGTGGTTTACCAAAACTATTTc                >  1:256386/1‑239 (MQ=255)                                                                                                                                                                                                                tgGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGCAAGCGGTGGTTGAAGCCCTGAAAGGTTAACCGGAAGTGAACGATTCTTTCGACGGCGAGGACCGGGTTTACCACAACAAATTc                >  2:1656/1‑239 (MQ=255)                                                                                                                                                                                                                tgGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCTACGGCGATGACGCGGTTTACCACAACTATTTc                >  1:89471/1‑239 (MQ=255)                                                                                                                                                                                                                        cGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGACGGCGATGACGTGGTTTACCACAACTATTTCGACGTCAg        >  1:386768/1‑239 (MQ=255)                                                                                                                                                                                                                         caccgtgggAA‑ACAAAAAA‑CCCCCCCGCCATGTTGGACCA‑GTTCAACGAAGTCAAAAAGAAGTCGATTATGGATTTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCCCGGTATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGACGGCGATGACGTGGTTTACCACAACTATTTCGACGTCAGc       <  2:89471/232‑1 (MQ=255)                                                                                                                                                                                                                            tgctgGCA‑GCGAAAAC‑CTCCACCGCCATGCT‑GACCACGTTCCACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGACGTGAACGCTTCAATCGACGGCGGTGACGTGTTTTACCACCACCATTTCCACCTCAgacgg    >  2:676380/1‑235 (MQ=255)                                                                                                                                                                                                                             ggtgGAAGGCGAAACA‑CCCCACCGTCATGCT‑GGCCGCGTTCAAGGGAGTCAACATGAACCCGGTTATGGAGCGGCGTAAGCGGTAAGGTGAAGCGTGTGAAAAACGCCACGGCATCCGTCTTGTCTTTAGGGTCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTAGCGACGGCGATGACGTGGTTTGCCACAACTATTTCGACGTCAGCATg    <  2:3404/237‑1 (MQ=255)                                                                                                                                                                                                                                tggaacgaaaaac‑CCCCCCCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGACGGCGATGACGTGGTTTACCACAACTATTTCGACGTCAGCATGGc  <  1:617982/230‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                 CCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGTGTGGGTG‑GTCGTCTGACTCGTGAAGATGTGGAAAAACA‑TCTGGCGAAAGCCCCGGCGAAAGAGTCT‑GCTCCGGCAGC‑GGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGACGGCGATGACGTGGTTTACCACAACTATTTCGACGTCAGCATGGC  >  NC_000913/761881‑762336

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