Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,984,689	T→C	23.1%	E301G (GAG→GGG)	yqeF ←	short chain acyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,984,689	0	T	C	23.1%	50.4 / 2.9	41	E301G (GAG→GGG)	yqeF	short chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base T (5/23);  new base C (6/4);  total (13/28)
Fisher's exact test for biased strand distribution p-value = 1.93e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTG‑CAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAG‑CCAACTGCCAGCCT‑ACACGCTCCAGGCAACGGCGGG‑TCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACA‑CCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCT  >  NC_000913/2984465‑2984907                                                                                                                                                                                                                                    |                                                                                                                                                                                                                             ttGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTG‑CAGCAAACGCTTCCTTAGCCTCGATAAGATCGCCTTCAG‑CCAACTGCCAGCCT‑ACACGCTCCAGGCAAcggcgg                                                                                                                                                                                                                >  2:104976/1‑239 (MQ=255)  tGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGCCATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTG‑CAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAG‑CCACCTCCCACCCT‑ACACACACCCGGCACcggcgg                                                                                                                                                                                                                >  1:257710/1‑238 (MQ=255)  tGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGCCATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTG‑CAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTAAG‑CCAACCGCCAACCA‑AAAACCCCCAGGCAACGGCggg                                                                                                                                                                                                               >  1:775259/1‑239 (MQ=255)  tGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGCCGCTCATCCCACCCAAGCATCTTTCCGACCGAAAGCGCCGTTC‑CAGCAACCGCTTCATTAACCTCGTTAAGATACAATTCCC‑CAACCACCCACCCC‑ACACACCCCAGCCAAcggcgg                                                                                                                                                                                                                >  1:384773/1‑238 (MQ=255)  tGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTG‑CAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTAAG‑CCAACTGCCAACCT‑AACCGCACCAGGCAACGCCggg                                                                                                                                                                                                               >  2:401262/1‑239 (MQ=255)  tGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTG‑CAGCAAACGCCTCATTACCCTCGATAAGATCGCCTTCAG‑CCAACTGCCACCCC‑AAACCCCCCAGGAAACgggggg                                                                                                                                                                                                               >  1:848948/1‑239 (MQ=255)    ccAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAACGCCTGGC‑CAGCAAACCCTTCATTAGCCCCGATAAGAACGCATTCAA‑CCAACGCCAAGCCTAACAC‑CCCCAGGCAAAGGCggggg                                                                                                                                                                                                             >  1:422554/1‑237 (MQ=255)         ccTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTG‑CAGCAACCGCTTCATTAGCCTCGATAAGATCGACTTCAG‑CCAACTGCCAGCCT‑ACCCGCTCCAGGCAACGGCggggccccat                                                                                                                                                                                                        >  2:28228/1‑232 (MQ=255)             tGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCTACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTG‑CAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAG‑CCAACTCCCAGCCT‑ACACGCTCCAGGCAAaggcgagctgcccccac                                                                                                                                                                                                     >  1:792379/1‑226 (MQ=255)             tGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTG‑CAGCAACCGCTTCATTAGCCTCGATAAGATCGACTTCAG‑CCAACTGCCAGCCT‑ACACGCCCCAGGCAACGGCGGG‑TCGCacaacc                                                                                                                                                                                                    >  2:457753/1‑233 (MQ=255)              gCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTTTG‑CAGCAAACGCTTCATTAGCCTCGATAAGATCGCCTTCAG‑CCACCTCCCAGCCT‑AAACCCTCCAGGCACCGACGGG‑CCGCAtgcccc                                                                                                                                                                                                   >  2:137820/1‑233 (MQ=255)                gggCCTTACGTTTCACCATTTCATGAACCAGAGAACCCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTG‑CAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAG‑CCAACTGCCAGCCT‑ACACGCTCCAGGCAACGGCGGG‑TCGCATACACCgg                                                                                                                                                                                                 <  2:752339/238‑1 (MQ=255)                  ggcttggttttCCCCTTTTCATGACCACGGGAACCCAAGATTGGGACACCGGAAGTGGCTTCGGGGTTGCCGCGGTGGATCGCGCCAACATTGACATTGACCCTACTCTCCTCCCACTCAAGCATCTAGCAAACCGAAAGCGCCTGGG‑CCGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAG‑CCAACTGCCAGCCT‑ACACGCTCCAGGCAACAGCGGG‑TCGCATACACCGgcg                                                                                                                                                                                               <  1:137820/230‑1 (MQ=255)                  gctttACGTTTACCCATTTCAAGGACCCGAGAAACCAGGATTCGGCAACCGGGAGCGCCTATCCGGGGTCCGCGTGCGATCCCGCCACCATTGGCATTGGCCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTG‑CAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAG‑CCAACTGCCAGCCT‑ACACGCTCCAGGCAACGGCGGG‑TCGCATACACCGgcg                                                                                                                                                                                               <  1:96124/235‑1 (MQ=255)                  gCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTG‑CAGCAAACGCTTCATTAGCCTCGATAAGATCGCCTTCAG‑CCACCTCCCAGCCT‑CCACCCCCCCGCCAACGGCGGG‑CCGCATGCAcccaccc                                                                                                                                                                                              >  1:17366/1‑234 (MQ=255)                                 aTTTCAAGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTG‑CAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAG‑CCAACTGCCAGCCT‑ACACGCTCCAGGCAACGGCGGG‑TCGCATACACCGGCGCAATTCCCATCAATGc                                                                                                                                                                               <  2:280072/239‑1 (MQ=255)                                                               acaCCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTG‑CAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAG‑CCAACTGCCAGCCT‑ACACGCTCCAGGCAACGGCGGG‑TCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACA‑CCAACGCTGGCAAATgcgc                                                                                                                                                 <  2:74211/236‑1 (MQ=255)                                                                                        ccggggCGATCGCGCCACCATTGACATTGACCCGACGCTCAACCCACTCAAGCATCTTGACGACCGAAAGCGCCTGTG‑CAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAG‑CCAACTGCCAGCCT‑ACACGCTCCAGGCAACGGCGGG‑TCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACA‑CCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTc                                                                                                                       <  2:384773/234‑1 (MQ=255)                                                                                                              gacattggcCCGACGCCCAACCCACCCAAGCATCTTGCCAACCGAAAGCGCCTTTG‑CAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAG‑CCAACTGCCAGCCT‑ACACGCTCCAGGCAACGGCGGG‑TCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACA‑CCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGc                                                                                                 <  2:872502/239‑1 (MQ=255)                                                                                                                                  ccccCTCAAGCATCTTGCCAACCGAAAGCGCCTGTG‑CAGCAAACGCTTCATTAGCCTCGATAAGATCGCCTTCAG‑CCAACTGCCAGCCT‑ACACGCTCCAGGCAACGGCGGG‑TCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACA‑CCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCt                                                                             <  2:817347/239‑1 (MQ=255)                                                                                                                                       tctaacaccTTTCCACCCGAAAGCGCCTGTG‑CAGCAAACGCTTCATTCGCCTAGATAAGATCGACCTCAG‑CCAACGGCCAGCCT‑ACACGCTCCAGGCAACTGCGGG‑TCGCATACACCGGCGTAATTCCCCTCAACGCCGGATCTACA‑CCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGcc                                                                        <  1:145374/231‑1 (MQ=255)                                                                                                                                       cacgtctctttcccaacc‑aaAGCGCCTCTG‑CACACAACGCTTCACTCGCCTCGATAAGGACTACTTCAG‑CCCACTGCCAGACT‑ACCCGTCCCAGGCAACGGCGGG‑TCGCATACACCGGAGCAATTCCCAGCCATGCCGGATCTACA‑CCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCa                                                                       <  1:578542/227‑1 (MQ=255)                                                                                                                                        ccagcctctctccAACCGAAATCGCCGTTC‑CACCAAACGCTTCCTGAGCCCCAATAAAATCGCCTTACG‑CCACCTGCCACCCA‑ACACGCCCCAGGCAACGGCGGG‑TCGCATACGCCGGCGCAATTCCCCTCCATGCCGGATCTACA‑CCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCa                                                                       <  1:852179/228‑1 (MQ=255)                                                                                                                                        caaaccacTTGCACACCAAAAGCCGCTGTG‑CACGCAACGCTTCTATAGCCTCGATACAACCGGCTCCAG‑ACAACTCCCATCCT‑ACACGCTCCAGGCAACGGCGGG‑TCGCATGCGCCGGCGCAATTCCCAGCAATGCCGGATCTACA‑CCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCa                                                                       <  1:451487/232‑1 (MQ=255)                                                                                                                                        cAAGCATCTTTCCAACCGAAAGCGCCTGTG‑CAGCAAACGCTTCCTTAGCCTCGATAAGATCGACTTCAG‑CCAACTGCCAGCCT‑ACACGCTCCAGGCAACGGCGGG‑TCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACA‑CCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCa                                                                       <  2:53051/239‑1 (MQ=255)                                                                                                                                        cAAGCATCTTGCCAACCGAAAGCGCCTGTG‑CACCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAG‑CCAACTGCCAGCCT‑ACACGCTCCAGGCAACGGCGGG‑TCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACA‑CCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCa                                                                       <  2:313214/239‑1 (MQ=255)                                                                                                                                               cTTGCCAACCGAAATCGCCTGTG‑CAGACAACGGTTCATTAGCCTGGATAAGAGCGACTTCAG‑CCAACTGCCAGCCT‑ACACGCTCCAGGCAACGGCGGG‑CCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACA‑CCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTa                                                                <  2:848931/239‑1 (MQ=255)                                                                                                                                                ttGCCAAACCAAAGCCCCTGGG‑CAGCAAACGCTTCATTTGCCCCCATAAGACCCACCTCAG‑CCAACTGCCAGCCT‑ACCCGCTCCGAGCAACGGCGGG‑TCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACA‑CCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTAt                                                               <  2:534461/239‑1 (MQ=255)                                                                                                                                                 tGCCAACCGAAAGCGCCTGTG‑CAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAG‑CCAACTGCCAGCCT‑ACACGCTCCAGGCAACGGCGGGCCCGCA‑ACAACGCCCCCATTTCCCTCCACTCCCGCGCTACACCCAACCCGGGCAAA‑CCCCTGATGCTGGCCCGCACCGCTACAAACAACCCCCGGTCTTTGGCTTCGCCCATCATCATTACCCCaccagcccacaactgtatg                                                              >  1:313214/1‑221 (MQ=255)                                                                                                                                                    cAACCGAAAGCGCCTGTG‑CAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAG‑CCAACTGCCAGCCT‑ACACGCTCCAGGCAACGGCGGG‑TCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACA‑CCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGga                                                            <  2:399316/238‑1 (MQ=255)                                                                                                                                                     acccGAAAGCCCCTGTG‑CAGCAACCCCTTAATTCGCCTCGACAAGACCGCCTTCAG‑CCAACGGCCAGCCT‑ACCCGCCCCAGGCAACGGCGGG‑TCGCATACGCCGGCGCAATTCCCATCAATGCCGGATCTACA‑CCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGgatg                                                          <  1:620637/237‑1 (MQ=255)                                                                                                                                                     aaCCGAAAGCGCC‑GTGTCAAAAAACGCTTGATTAGCCCCGATAAGATCGCCTTCAG‑CCAACTGCCAGCCT‑ACACGCTTCAGGCAACGGCGGG‑TCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACA‑CCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGgatg                                                          <  2:482405/239‑1 (MQ=255)                                                                                                                                                               ctcTGGG‑CAGCAAAAGCTTGTTTAGCCTCGATAAGATCG‑CCTCAGCCCACCTGCCAGCCT‑ACACGCTCCAGGCAACGGCGGG‑TCGCATACACCGGCGCAATGCCCATCAATGCCGGATCTACA‑CCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTAcc                                                <  2:646245/237‑1 (MQ=255)                                                                                                                                                                ccTGTG‑CAGCAAACGCTTCATTAGCCTCGAGAAGATCGACTTCAG‑CCAACTGCCAGCCT‑ACCCTCTCCAGGCAACGGCGGG‑TCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACA‑CCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTAcc                                                <  2:583694/238‑1 (MQ=255)                                                                                                                                                                           aaaCGCTTCATTAGCCTCGATAAGATCGACTTCAG‑CCAACTGCCAGCCT‑ACACGCTCCAGGCAACGGCGGG‑TCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACA‑CCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTca                                        <  2:257710/236‑1 (MQ=255)                                                                                                                                                                                   cATTAGCCTCGATAAGATCGACTTCAG‑CCAACTGCCAGCCT‑ACACGCTCCAGGCAACGGCGGG‑TCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACA‑CCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCt                                                                                                     <  2:258821/167‑1 (MQ=255)                                                                                                                                                                                             caTAAGCTCGCCTTCAC‑CCACCCGCCATCCT‑ACACGCGCCAGGCAACGGCGGG‑TCGCATACCCCGGCGCAATTCCCATCAAGGCCGGATCTACC‑CCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATg                   <  1:104976/238‑1 (MQ=255)                                                                                                                                                                                              aaaaGATCGACTTCAG‑CCAACTGCCAGCCT‑ACCCGCTCCAGGCAACGGCGGG‑TCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACA‑CCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGa                  <  2:609931/237‑1 (MQ=255)                                                                                                                                                                                                        ctcccg‑CCACCTGCCCCCCT‑ACCCGCTCCAGGCACCGGGGGG‑CCCCATACCCCGGCGTCATTCCCACGAATGCCGTATCTACA‑CCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTaaa        <  1:457753/234‑1 (MQ=255)                                                                                                                                                                                                         ttCAG‑CCAACTCCCCGCCT‑ACCCGCCCCAGGCAACGGCGGG‑CCGCATAACCCGGCGCAATTCCCATCAATGCCGGATCCACA‑CCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTaaa        <  2:268541/238‑1 (MQ=255)                                                                                                                                                                                                               ccAACCGCCAGCCT‑AACCGCTCCAGGCAACGGCGGG‑TCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACC‑CCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCt  <  2:775259/239‑1 (MQ=255)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                             TTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTG‑CAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAG‑CCAACTGCCAGCCT‑ACACGCTCCAGGCAACGGCGGG‑TCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACA‑CCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCT  >  NC_000913/2984465‑2984907

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