Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	762,113	T→C	25.5%	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.5%	78.1 / 3.5	49	S198P (TCC→CCC)	sucB	dihydrolipoyltranssuccinase
Reads supporting (aligned to +/- strand):  ref base T (24/10);  new base C (2/12);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 4.83e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGT‑GTGGGTG‑GTCGTC‑TGACTCGTGAAGATGT‑GGAAAAA‑CA‑T‑CTGGCGAAAGCCCCGGCGAAAGAGTCT‑GCTCCGGCAGC‑GGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGACGGCGATGACGTGGTTTACCACAACTATTTCGACGTCAGCATGGC  >  NC_000913/761881‑762336                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                 ccTGTTGGTTGAACACATTTCGGGCTCCACGGCCTTTAAAGGCCTCGGT‑GTGGGTG‑GTCGCT‑TGACCGGGGAAGAGGT‑GGAAAAA‑CA‑T‑CGGGCGAAAGCCCCGGCAAAAGAGTCG‑GCTCCGGCAGG‑GGCTGCTCCGGCGGCGCAACCGGCTTTGGCTGCCCGTAGTGAAAACCGTGCCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCg                                                                                                                                                                                                                            <  1:369033/239‑1 (MQ=255) ccTGCTGGCTGAACACAATCTCGACGCCAGGGCAATTAAAGGCACCGGT‑GTGGGTG‑GTCGTC‑TGACTCGTGAAGATGT‑GGAAAAA‑CA‑T‑CTGGCGAAAGCCCCGGCGAAAGAGTCT‑GCTCCGGCAGC‑GGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCg                                                                                                                                                                                                                            <  2:590688/239‑1 (MQ=255)        gcggtacacAATCTCCGCGCCCGCGCCCTTTAAGGGCCCCGG‑GTGGGGGTGTCGCT‑TGAC‑CGTTAAAAGGT‑GGGAAAA‑CATT‑CTGGAGAAACCCCCGG‑GAAAGGATCTGGCTCCGGCAGC‑GGTGCCCCCGGGGGCCCACCCGGTTTTGGCTGCCGGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCAAAAAC‑CCCCACCGCCATGc                                                                                                                                                                                                                      <  1:274272/234‑1 (MQ=255)           gaaccaaatgccgccGCCCGGGCCATTAAGGGC‑CCGGTGGGGGGTG‑GT‑GTCTTGACTCGGGAAGATGT‑GAAAAAC‑CA‑T‑CGGGCGAAACCCCCGGCGAAAGGTGCT‑GCTCCGGCGGC‑GGCTGCTCCGGCGGCGCAACCGGCTGTGGCGGCACTTAGTGAAAAACGTTTCCCGAGGACGCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCGGGAA‑GGGAAAAA‑CTCCGCCGCCATGCT‑GAc                                                                                                                                                                                                                 <  1:811230/225‑1 (MQ=255)                           cccgggcagtaaaagtggcgggT‑GTGGGTG‑GGGTGA‑TGACGCGTGAACGGGT‑GGAAACA‑CA‑G‑GTGGGGAAAGCCCCGGGGAAGGGTTTT‑GCCCCGGCGGC‑GGCTGCTCGGGCGGCGAAACCGGCTCGGGCTGCCGTTGGTGAAAAAGGTGTCACGATAACCCGCCGGGGTAAGGGTGGGGCAGGGTGTCTGCGGGAG‑GCAAAAAA‑CGCCACCGCCATGCT‑GACCACGTTCAACGAAGt                                                                                                                                                                                                  <  1:262265/218‑1 (MQ=255)                               cgcCATTAAAGGCACCGGT‑GTGGGTG‑GTCGTC‑TGACTCGTGAAGATGT‑GGAAAAA‑CA‑T‑CTGGCGAAAGCCCCGGCGAAAGAGTCT‑GCTCCGGCAGC‑GGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACa                                                                                                                                                                                             >  2:432896/1‑239 (MQ=255)                                    ttAAAGGCCCCGGT‑GTGGGTG‑GGCGTC‑TGACTCGTGAAGATGT‑GGAAAAA‑CA‑T‑CTGGCGAAAGCCCCGGGAAAAGAGTCT‑GCTCCGGGAGC‑GGCCGCTCCGGCGGCGCAACCGGCTCTGGCTGCCCGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAg                                                                                                                                                                                        <  2:626918/239‑1 (MQ=255)                                    taaagagcccgtgt‑gtGGGGG‑GTCGTT‑TAACTGGGAAGGATGG‑GGAAAAC‑CC‑T‑CGGGCGAAACCCCCGGCGAAAGATTCT‑CCCCCGGCGGC‑GGCTGCTCGGGCGGCGCAACCGGCTCTGGACGCCCGGAGTGAAAAACGGGTCCGGATGACTCGCCTGCGTAAGCGTGCGGCAGGGCGTCGGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAg                                                                                                                                                                                        <  1:72726/227‑1 (MQ=255)                                            aCCGGT‑GTGGGTG‑GTCGTC‑TGACTCGTGAAGATGT‑GGAAAAA‑CA‑T‑CTGGCGAAAGCCCCGGCGAAAGAGTCT‑GCTCCGGCAGC‑GGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTAt                                                                                                                                                                                >  2:510449/1‑239 (MQ=255)                                            aCCGGT‑GTGGGTG‑GTCGTC‑TGACTCGTGAAGATGT‑GGAAAAA‑CA‑T‑CTGGCGAAAGCCCCGGCGAAAGAGTCT‑GCTCCGGCAGC‑GGCTGCTCCGGCGACGCAACCGGCTCTGGCTGCACGCAGTGAAAACAGTGTTCCGGTAACTCGCCTGCGGAAGCGTGTGGCAGAGGGTCTGCTTGCA‑GCGCAAAA‑CCCCCCCCCCCTCTT‑GGCCACGTTCAACGAAATCAACAGAAAACCGAtatg                                                                                                                                                                                >  1:80145/1‑236 (MQ=255)                                              cGGT‑GTGGGTG‑GTCGTC‑TGACTCGTGAAGATGT‑GGAAAAA‑CA‑T‑CTGGCGAAAGCCCCGGCGAAAGAGTCT‑GCTCCGGCAGC‑GGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCACCGAAGTCAACATGAAGCCGATAATga                                                                                                                                                                              >  2:328050/1‑238 (MQ=255)                                              cGGT‑GTGGGTG‑GTCGTC‑TGACTCGTGAAGATGT‑GGAAAAA‑CA‑T‑CTGGCGAAAGCCCCGGCGAAAGAGTCT‑GCTCCGGCAGC‑GGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATgg                                                                                                                                                                              >  2:423124/1‑239 (MQ=255)                                                 gggggggtgtgtcgtc‑TGAATCGGGAAGATGT‑GAAAAAC‑CA‑G‑CTGGCGAAAGCCCCGGCGAAATGGCCT‑GTCCGGGAAGG‑GGCTGCTCCGGCGGCGCAACCGGCTCTGGCGGCGCGTAGTGAAAAACGTGTCCCGATGCCTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGa                                                                                                                                                                             <  1:575493/229‑1 (MQ=255)                                                            tcgtc‑TGACTCGTGAAGATGT‑GGAAAAA‑CA‑T‑CTGGCGAAAGCCCCGGCGAAAGAGTCT‑GCTCCGGCAGC‑GGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGGGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCAGGTTCAACGAAGGCAAAAGAAAGCAGATTATGGAGTTGGGTAAga                                                                                                                                                                  >  1:869791/1‑238 (MQ=255)                                                             cgtc‑TGACTCGTGAAGATGT‑GGAAAAA‑CA‑T‑CTGGCGAAAGCCCCGGCGAAAGAGTCT‑GCTCCGGCAGC‑GGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTGAGGGAAAACGGGTCCCGATGACTCGCCTGCGTAGGCGTGGGCCGGAGCGTCGTCTGGAA‑ACGAAAAA‑CCCCACCCCCCAGCT‑GACACCGTTCAACAACGCCCACCTACAGCCCATTTTGGATCTGCGTCAgaa                                                                                                                                                                 >  1:695151/1‑237 (MQ=255)                                                               tc‑TGACTCGTGAAGATGT‑GGAAAAA‑CA‑T‑CTGGCGAAAGCCCCGGCGAAAGAGTCT‑GCTCCGGCAGC‑GGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAga                                                                                                                                                               >  2:25832/1‑238 (MQ=255)                                                                ctttgaactagtaGATGTGGGAAAAACCA‑T‑GTGGAGAAGCCCCCGGCAAAAGGGCTC‑GCCCCGGTAGC‑GGCTGCTCGGGCGGCGAACCCGGCTTGGGCCGCAGGTAGGGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCTTGTGGCAGCGCGTTTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCa                                                                                                                                                                 <  1:499353/227‑1 (MQ=255)                                                                cttcaCTCGGGAAGTTGT‑GAAAAAC‑CC‑T‑CTGGCGAAACCCCCGGGGAAGGAGCCT‑GCCCCGGAGGC‑GGCTCCCCCGGGCGCGCACCCGGCTCTGGCTGCACGTAGGTAAAAACGTGTCCGGAGGCCCCGCCTCCGTAAGCTGGTGGCAGAGCCTCTGCTGGAA‑GGGAAAAA‑CCCAACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGt                                                                                                                                                               <  1:425605/235‑1 (MQ=255)                                                                cgtctgccctggagagttggggaaaa‑CA‑TGCGGGCGAAACCCCCGGCGAAAGAGTCT‑GCTCCGGCAGG‑GGCGGCCCCGGCGGCGAACCGGGCTCTGGCTGCCCGTGGTGAAAAACGTGCCCCGATGACTCGCCTGCGTAAGCGGGTGGCAGAGCGTCTGCTGGAA‑GAGAAAAA‑CCCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCa                                                                                                                                                                 <  1:423124/217‑1 (MQ=255)                                                                                gt‑gGAAAAA‑CA‑T‑CTGGCGAAAGCCCCGGCGAAAGAGTCT‑GCTCCGGCAGC‑GGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGa                                                                                                                                               >  2:620875/1‑239 (MQ=255)                                                                                   ggAAAAA‑CA‑T‑CTGGCGAAAGCCCCGGCGAAAGAGTCT‑GCTCCGGCAGC‑GGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAAAATGAAGCCGATTATGGATCTGCGTAAGAAGTACGGGGAAAGGTTTGaaa                                                                                                                                             >  1:519536/1‑239 (MQ=255)                                                                                   ggAAAAA‑CA‑T‑CTGGCGAAAGCCCCGGCGAAAGAGTCT‑GCTCCGGCAGC‑GGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGCCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTTTGGATCTGCGAAAGCGTTACGTTGAAGCGTTTGaaa                                                                                                                                             >  2:892852/1‑239 (MQ=255)                                                                                              t‑CTGGCGAAAGCCCCGGCGAAAGAGTCT‑GCTCCGGCAGC‑GGCTGCTCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGGGAAGCGTTTGcaaaaacccccg                                                                                                                                    >  2:298971/1‑232 (MQ=255)                                                                                                              gggaaaaGGGTCT‑GCTCCGGCAGC‑GGCTGCTCCGGCGGCGCAACCGGCTCGGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGGGCGTCTGCTGGAA‑GCGAAAAA‑CCCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCt                                                                                                                     <  2:148783/235‑1 (MQ=255)                                                                                                                        gtgttgcccGGGAGCGGGCTGCCCCGGCGGCCCAACGGGCTTGGGCGGCCGGTAGGGAAAAACGTGTCCCGATGACCCGCCTGGGTAAGCGTGGGGAAGGGCGTGTGCTGGAA‑GCAAAAAC‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGt                                                                                                               <  1:510449/230‑1 (MQ=255)                                                                                                                                         gctgctCCGGCGGCGCAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCCCGTTCAACGAAGGCAACATGAACCCGATTATGGATCCGGGGAAGCAGTACGGTGAAGCATTTGAAAAAAGCAACGGCACCCGTCTGGGCGTTATGTCCTTCTACGGGAAAGCggg                                                                                            >  1:404031/1‑238 (MQ=255)                                                                                                                                                       gcAACCGGCTCTGGCTGCACGTAGTGAAAAACGTGTCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGGGGTTGAAGCCCTGa                                                                              >  2:335790/1‑239 (MQ=255)                                                                                                                                                                     cTGAACGAAGTAAAAATCGGCTCCAGATGACTCGCCTGCGTAAGCGTGGGGCAGAGCGTCTGCTGAAA‑GCGAAAAA‑CCCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGa                                                                 <  2:695151/238‑1 (MQ=255)                                                                                                                                                                                      cttctccctctgccccccctcggTAGGCTTGTGGCAGCGTGTCTGGCGGAA‑GCAAAAAC‑CCCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGa                                               <  1:432896/217‑1 (MQ=255)                                                                                                                                                                                         gtCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACACGGAAGTGAACGCTTCTATCGACg                                             >  2:88050/1‑238 (MQ=255)                                                                                                                                                                                         gtCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCCGAAACGTTACCCGGAAGTGAACGCTTCTATCGACgg                                            >  2:468572/1‑239 (MQ=255)                                                                                                                                                                                         gtCCCGATGACTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGGGGGTGAAGCCCTGAAACGTTACCCGGAAGGGAACGCTTCTATCGACgg                                            >  1:617982/1‑239 (MQ=255)                                                                                                                                                                                                  ccTCGCCGGCGTAAGCGTGGGGCAGGGCGTCTGTGGGAA‑GCGAAAAC‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGACGGCGATGACg                                    <  2:519536/237‑1 (MQ=255)                                                                                                                                                                                                  aCTCGCCTGCTTAAGCTGGGGGCAGGGCGTCGGCTGAAA‑GCGAAAAC‑CCCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGACGGCGATGACg                                    <  2:80145/238‑1 (MQ=255)                                                                                                                                                                                                  aCTCGCCTGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGGGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTACCGagcgcaaagccgg                                   >  1:99279/1‑227 (MQ=255)                                                                                                                                                                                                         ttgtgatgCGTTTGGCAGGGCTTTTGTTGG‑A‑GAGAAAAACCCCCACCTCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGACGGCGATGACGTGGTTTAc                            <  1:328050/232‑1 (MQ=255)                                                                                                                                                                                                         tGCGTAAGCGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGACGCCCTGAACCGTTCCCCGGAAGGGAACCCTTCTATCGACggggggggctgggttta                             >  1:796206/1‑223 (MQ=255)                                                                                                                                                                                                              aggCGTGTGGCATAGCGTTGGCTGGAA‑GCAAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTAAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGACGGCGATGACGTGGTTTACcaca                        <  2:879023/236‑1 (MQ=255)                                                                                                                                                                                                              aaGCGTGGGGCAAAGCGTTCGGGGGAA‑GCGAAAAA‑CCCAGCCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGACGGCGATGACGTGGTTTACCACaa                       <  2:678520/239‑1 (MQ=255)                                                                                                                                                                                                                 cGTGTGGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGGTCAACGAAGTCAACATGAAGCCGATTATGGATATGCGTAAGCAGAACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCGTTATGTCCTTCTACGCGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGACGGATTCTATCGACGGCGATGACGTGGGTTACCACAACTa                    >  2:831781/1‑239 (MQ=255)                                                                                                                                                                                                                     tgGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCTGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGACGGCGCAGACGGGGTTTACAACAAATATTTc                >  1:368308/1‑239 (MQ=255)                                                                                                                                                                                                                     tgGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTAGGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGACGGCGATGACGTGGTTTACCAAAACTATTTc                >  2:256386/1‑239 (MQ=255)                                                                                                                                                                                                                     tgGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGCAAGCGGTGGTTGAAGCCCTGAAAGGTTAACCGGAAGTGAACGATTCTTTCGACGGCGAGGACCGGGTTTACCACAACAAATTc                >  1:1656/1‑239 (MQ=255)                                                                                                                                                                                                                     tgGCAGAGCGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCTACGGCGATGACGCGGTTTACCACAACTATTTc                >  2:89471/1‑239 (MQ=255)                                                                                                                                                                                                                             cGTCTGCTGGAA‑GCGAAAAA‑CTCCACCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGACGGCGATGACGTGGTTTACCACAACTATTTCGACGTCAg        >  2:386768/1‑239 (MQ=255)                                                                                                                                                                                                                              caccgtgggAA‑ACAAAAAA‑CCCCCCCGCCATGTTGGACCA‑GTTCAACGAAGTCAAAAAGAAGTCGATTATGGATTTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCCCGGTATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGACGGCGATGACGTGGTTTACCACAACTATTTCGACGTCAGc       <  1:89471/232‑1 (MQ=255)                                                                                                                                                                                                                                 tgctgGCA‑GCGAAAAC‑CTCCACCGCCATGCT‑GACCACGTTCCACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGACGTGAACGCTTCAATCGACGGCGGTGACGTGTTTTACCACCACCATTTCCACCTCAgacgg    >  1:676380/1‑235 (MQ=255)                                                                                                                                                                                                                                  ggtgGAAGGCGAAACA‑CCCCACCGTCATGCT‑GGCCGCGTTCAAGGGAGTCAACATGAACCCGGTTATGGAGCGGCGTAAGCGGTAAGGTGAAGCGTGTGAAAAACGCCACGGCATCCGTCTTGTCTTTAGGGTCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTAGCGACGGCGATGACGTGGTTTGCCACAACTATTTCGACGTCAGCATg    <  1:3404/237‑1 (MQ=255)                                                                                                                                                                                                                                     tggaacgaaaaac‑CCCCCCCGCCATGCT‑GACCACGTTCAACGAAGTCAACATGAAGCCGATTATGGATCTGCGTAAGCAGTACGGTGAAGCGTTTGAAAAACGCCACGGCATCCGTCTGGGCTTTATGTCCTTCTACGTGAAAGCGGTGGTTGAAGCCCTGAAACGTTACCCGGAAGTGAACGCTTCTATCGACGGCGATGACGTGGTTTACCACAACTATTTCGACGTCAGCATGGc  <  2:617982/230‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                 CCTGCTGGCTGAACACAATCTCGACGCCAGCGCCATTAAAGGCACCGGT‑GTGGGTG‑GTCGTC‑TGACTCGTGAAGATGT‑GGAAAAA‑CA‑T‑CTGGCGAAAGCCCCGGCGAAAGAGTCT‑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: ‑