Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,130,929	T→C	17.7%	S366P (TCA→CCA)	metL →	Bifunctional aspartokinase/homoserine dehydrogenase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,130,929	0	T	C	17.7%	120.6 / 2.2	63	S366P (TCA→CCA)	metL	Bifunctional aspartokinase/homoserine dehydrogenase 2
Reads supporting (aligned to +/- strand):  ref base T (16/33);  new base C (11/2);  total (27/36)
Fisher's exact test for biased strand distribution p-value = 1.19e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCC‑GGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCG‑CCAGTTGCT‑GCAATTTTGCTACACCTCAGAAGTGGCCG‑ACAGTGCGCTG‑AAAATCCTCGACGAAG‑CGGGATTACCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCGGTCGAATTTACCTGGCAGTCCGATGACGGCATCAGCCTGGTGGCAGTACTGCGCACCGGCCCGACCGAAAGCCTGA  >  NC_000913/4130694‑4131166                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                  tttccccccccgcgATCACGGTCCCCCGCGCTTTGAACCCGTG‑GGCCCTCCGGGTACGGGTGCGCGTTTTGGCCCCGGCCAGGATGTGGTCTGTTTGATTGAGTTTCGGGGTCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCG‑CCAGTTGCT‑GCAATTTTGCTACACCTCAg                                                                                                                                                                                                                                                <  1:745466/226‑1 (MQ=255)         cGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCC‑GGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTTGGGGACATAACGATCG‑CCAGTTGCTCGCAATTTT‑TTACACCCCCGAAATGGc                                                                                                                                                                                                                                         >  2:79427/1‑238 (MQ=255)              gATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCC‑GGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGCACATAAAAAACGCCCCGTTG‑T‑GCAACTTTGTTCCACCCCAGAAGAGGccacaag                                                                                                                                                                                                                                   >  1:583580/1‑234 (MQ=255)               aTCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCC‑GGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACCTAACGCTTC‑CCAGTTGTA‑GCCATTTTGCTATAACTCAAAAGAGGCCc‑acaag                                                                                                                                                                                                                                 >  2:381264/1‑237 (MQ=255)                        ccACGCGCATTGAACGCGTGCTGGCCTCC‑GGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGAAAAAAACAATCC‑CCAATTGCG‑GCAATTTTGTTAAACCCCAGAAAGGGCCG‑CCaaggggcgtaaa                                                                                                                                                                                                                        >  1:875877/1‑228 (MQ=255)                                attgacccggttcTGGCCTTC‑GGTATTGGGGAGAGTATGGCCCCCAGCCAGGATGATTTCTGTTTGATTGAGTTCCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCG‑CCAGTTGCT‑GCAATTTTGCTACACCTCAGAAGTGGCCG‑ACAGTGCGCTG‑AAAATCCTcga                                                                                                                                                                                                               <  2:496960/227‑1 (MQ=255)                                     ccggttGCTGGCCCCG‑GGTACGGGTGCGCGTATTGTCCCCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCG‑CCAGTTGCT‑GCAATTTTGCTACACCTCAGAAGTGGCCG‑ACAGTGCGCTG‑AAAATCCTCGACGAAg                                                                                                                                                                                                          <  2:780548/234‑1 (MQ=255)                                          tGCTGGCCTCC‑GGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGGACATAACGATCG‑CCAGGTGCT‑GCAATTTTGCTACACCCCCGAAGGGGCCG‑ACCGTGGCCTG‑AAAATCCTCGACGAAA‑Cgggg                                                                                                                                                                                                    >  2:273826/1‑238 (MQ=255)                                          tGCTGGCCTCC‑GGCACTGGTGCGCGTATTGTCACCAGACACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGGTGGCGGACATAACCATCG‑CCACTTGCT‑GCCATTTTGCTACCACCCCGAACTGGCCGCACCGT‑CCCCG‑GAAATCCCCCACGAAA‑CGGGa                                                                                                                                                                                                    >  1:550588/1‑239 (MQ=255)                                           tccggCCTCC‑GGTACTGGTCGGCGTATTGCCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCG‑CCAGTTGCT‑GCAATTTTGCTACACCTCAGAAGTGGCCG‑ACAGTGCGCTG‑AAAATCCTCGACGAAG‑CGGGa                                                                                                                                                                                                    <  2:165948/235‑1 (MQ=255)                                                  ccgcgtgtttgtggtggtttttGCCACCCGCCAGGTTGTTGTCTGTTTGTTGGGGTTCCGGGGGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGCCCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCG‑CCAGTTGCT‑GCAATTTTGCTACACCTCAGAAGTGGCCG‑ACAGTGCGCTG‑AAAATCCTCGACGAAG‑CGGGATTACCTg                                                                                                                                                                                             <  1:850426/219‑1 (MQ=255)                                                               gcgcGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCG‑CCAGTTGCT‑GCAATTTTGCTACACCTCAGAAGTGGCCG‑ACAGTGCGCTG‑AAAATCCTCGCCGAAG‑CGGGATTACCTGGCGAACTGCGcc                                                                                                                                                                                 >  2:72581/1‑238 (MQ=255)                                                                  cgTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCG‑CCAGTTGCT‑GCAATTTTGCTACCCCCCAGAAGTGGCCG‑ACAGTGCGCTG‑GAAATCCTCGCGCAAA‑CGGGGTTACCTGGCGAACtcggccccg                                                                                                                                                                              >  2:656285/1‑230 (MQ=255)                                                                  cgTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGGACAAAACGATCG‑CCCACTGCC‑GCAATTTTTCTACACCCCAGAAAGGGCCG‑GCAGTGGCCGT‑AAAAACCCCGCCGAAA‑CGGGAGTACCCGGCGAACAGCGCCcgcg                                                                                                                                                                             >  1:877718/1‑239 (MQ=255)                                                                     aTTGTCCCCAGCCAGTATGTTGCTTTTTTAGTTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAGCGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCG‑CCAGTTGCT‑GCAATTTTGCTACATCTCAGAAGTGGCCG‑ACAGTGCGCTG‑AAAATCCTCGACGAAG‑CGGGATTACCTGGCGAACTGCGCCTGCGT‑Ca                                                                                                                                                                         <  2:112326/239‑1 (MQ=255)                                                                       tGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCG‑CCCGTTGTC‑GCCACTTTGCTCCACCCCGAGAGTGGCCG‑CCAGGGGGCCG‑GAAAACCTCG‑CGCAGACGGGATTACCCGGCggaacagcccccgcccaag                                                                                                                                                                        >  1:806037/1‑221 (MQ=255)                                                                               gCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCG‑CCAGTTGCT‑GCAATTTTGCCAAACCCCAGAAGTGGCCG‑GCAGGGGGCGG‑AAAATCCCCGACGAAG‑CGGGATTACCTGGGCAACTGGCCCTGCGC‑CAGGGGCTggg                                                                                                                                                                >  2:70246/1‑237 (MQ=255)                                                                                          gTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGACCC‑CCAGTTGTC‑GCAATTTTGCCACACCCCAAAAGGGGCCG‑ACAGGACGCTG‑AAAAACCTCTACGAAA‑CAGGGTTACCTGGCGAAACGCCCCCCCGG‑CAgggggggggggcggggggggg                                                                                                                                                    >  1:808518/1‑222 (MQ=255)                                                                                             tgtttgttgtgtgtcgtggTCCCCGTCAGTCAGGATTTAACACTGGCGCATAAAGAGATCGACCAAAGTCTGAAACGCGCGCAGGTACGGCCGCGCGCGGTGGGCGTCCATAACGATCG‑CCAGTTGCT‑GCAATTTTGATACGCCTCAGAAGGTGCCG‑TCAGTGCGCTG‑AAAATCCTCGACGAAG‑CGGGATTACCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCGCTGGTGGATATGGt                                                                                                                                                 <  1:653165/222‑1 (MQ=255)                                                                                             tGTTTGATTGAGTTTCAGGTGCCCGCCAGTAAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGGGTACATAACAAAAG‑CCAGTTGCT‑GCAAATTTGCTACACCTCAGAACGGGGCG‑AAAGGGCGCTG‑AAAAAACTCTCCGAAG‑CGGGGATAACAGGCTGAAAGCGCCTGCGTACAGGGGGTgggcacggggggggggg                                                                                                                                                  >  1:145046/1‑224 (MQ=255)                                                                                             tGTTTGATGGATTTTCATGTGCCCGCCAGGCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAACCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCG‑CCAGTTGCT‑GCAATTTTGCTACACCTCAGAAGTGGCCG‑ACAGTGCGCTG‑AAAATCCTCGACGAAG‑CGGGATTACCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCGCTGGTGGCGATGGt                                                                                                                                                 <  2:534411/239‑1 (MQ=255)                                                                                                  gtttttgttctgggTGCCCCCTAGTCAGGATTTCAAACTGGCGCATAAAGGGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCG‑CCAGTTGCT‑GCAATTTTGCTACACCTCAGAAGTGGCCG‑ACAGTGCGCTG‑AAAATCCTCGACGAAG‑CGGGATTACCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCGCTGGTGGCGATGGTCGGTg                                                                                                                                            <  1:72581/227‑1 (MQ=255)                                                                                                       aGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCG‑CCCGTTGCT‑GCAATTTTTCTACACCTCCGCAGTGGCCC‑GCCGAGCGCTG‑AAAAACCCCCACGGAG‑CCGGGTTACCTGGCGAACTGCGCCTGCGG‑CCGGGGCCGGCGCTGGGGGGGAGGGTCGGGGCAggg                                                                                                                                       >  1:549461/1‑238 (MQ=255)                                                                                                         tttCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGGGGACACAACGATCG‑CCAGGTGCT‑GGCATTTTTCTACACCTCAAAAGGGGCCG‑ACAGGGCGCCGCAAAAACCCCG‑CGCAGAGGGGTTTACCT‑GCGAACATCCCCCGCGG‑CACGGGCCGGCGCCGGGGGGGGGGGTCGGGGCCGGCg                                                                                                                                      >  1:105152/1‑238 (MQ=255)                                                                                                         tttCAGGTCCCCGCCAGTCAGGATTTCAAACGGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCG‑CCAGTTGCT‑GCAATTTTGCTACACCTCAGAAGTGGCCG‑ACAGTGCGCTG‑AAAATCCTCGACGAAG‑CGGGATTACCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCg                                                                                                                                      <  2:715759/238‑1 (MQ=255)                                                                                                         ttccAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCG‑CCAGTTGCT‑GCAATTTTGCTACACCTCAGAAGTGGCCG‑ACAGTGCGCTG‑AAAATCCTCGACGAAG‑CGGGATTACCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCg                                                                                                                                      <  2:19272/235‑1 (MQ=255)                                                                                                         ttcagggTGCCCGCCAGTGGGGATTTCAAACTGGCGCAGAAAGAGAACGGACAAATCCGGAAACGGGCGGAGGTAGGCCCGCTGGCGGGTGGCGGACATAACGATCG‑CAGGTTGGT‑GCAATTTTGCTACACCTGAGAAGTGGCCG‑ACAGTGCGCTG‑AAAATCCTCGACGAAG‑CGGGATGAGCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGt                                                                                                                                     <  1:65143/234‑1 (MQ=255)                                                                                                          ttCAGGGGCCCGGCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCG‑CCAGTTGCT‑GCAATTTTGCTACACCTCAGAAGTGGCCG‑ACAGTGCGCTG‑AAAATCCTCGACGAAG‑CGGGATTACCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGt                                                                                                                                     <  2:8738/238‑1 (MQ=255)                                                                                                                       caggcaggATTTCAAACGGGCACATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCG‑CCAGTTGCT‑GCAATTTTGCTACACCTCAGAAGTGGCCG‑ACAGTGCGCTG‑AAAATCCTCGACGAAG‑CGGGATTACCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGc                                                                                                                       <  2:806037/239‑1 (MQ=255)                                                                                                                                 ttCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGGGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCG‑CCAGTTGCT‑GCAATTTTGCTACACCTCAGAAGTGGCCG‑ACAGTGCGCTG‑AAAATCCTCGACGAAG‑CGGGATTACCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCa                                                                                                             <  2:877718/239‑1 (MQ=255)                                                                                                                                                aaaggggtcccccAAACCCTGAACCGGGGGGAGGTTCCGCCGGTGGCGGTTGGCGTACATACCGGTCG‑CCAGTTGCT‑GCAATTTTGGTACACCTCAGAAGTGGCCG‑ACAGTGCGCTG‑AAAAGCCTCGCCGAAG‑CGGGATTACCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGgcagca                                                                                              <  1:478525/228‑1 (MQ=255)                                                                                                                                                     tgtCGCCCAACCCCCGAAACGGGGGAAGGTCATCCCGGTGGGGGTTGGCTGACATAACGATCG‑CAAGTTGAT‑GCAATTTGGCTACACCCCAGAAGTGGCCG‑ACAGTGCGCTG‑AAAATCCTCGACGAAG‑CGGGATTACCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTg                                                                                          <  1:543052/236‑1 (MQ=255)                                                                                                                                                          aCCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTTGCGTACATAACGATCG‑CCAGTTGCT‑GCAATTTTGCTACACCTCCGAAGTGGCCG‑CCCGTGCGCTG‑AAAATCCTCGACCAAG‑CGGGATTAACTGGCGAACTGCGCCCGCGT‑CCAGGGCGGGGGCCGGTGGCGGTGGGCGGGGCAGGGGTCCCCCGCAACCCCCAGCATTTGCACCCCTTCCGGGGGAAACAGAAAgg                                                                                     >  1:293655/1‑238 (MQ=255)                                                                                                                                                              aaTCCTGAAACGCGCGCAGGGTCGCCCGCTGGCGTTTGGCGTACATAACGATCG‑CCAGTTGCT‑GCAATTTTGCTACACCTCAGAAGTGGCCG‑ACAGTGCGCTG‑AAAATCCTCGACGAAG‑CGGGATTACCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGc                                                                                <  2:86853/239‑1 (MQ=255)                                                                                                                                                                tCCTGAAACGCGCGCAGGTCCGCCCGCTGGCGGTTGGCGTACATAACGATCG‑CCAGTTGCT‑GCAATTTTGCTACACCTCAGAAGTGGCCG‑ACAGTGCGCTG‑AAAATCCTCGACGAAG‑CGGGATTACCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCg                                                                              <  2:790307/239‑1 (MQ=255)                                                                                                                                                                 ccTGAAACGCTTGCAGGTTCGCCCCGTGGCGGTTGGCGTACGTAACGATCT‑CCGGTTGCG‑GTAATTTTGCTACGCCTGAGAAGTGGCCG‑ACAGTGCGCTG‑AAAATCCTGGACGAAG‑CGGGATTGCCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCgg                                                                             <  1:70246/239‑1 (MQ=255)                                                                                                                                                                  cagaaaccggcgggggTACGCCCGCTGGGGGTTGGCGTACAATACGATCG‑CCAGTTGCT‑TCAATTTTGCTACCCCTGAGAAGTGGCCG‑ACAGTGCGCTG‑AAAATCCTCGACGAAG‑CGGGATTACCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCGGt                                                                            <  1:342980/225‑1 (MQ=255)                                                                                                                                                                     aaaccgcggcAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCG‑CCAGTTGCT‑GCAATTTTGCTACACCCCAGAAGTGGCCG‑ACAGTGCGCTG‑AAAATCCTCGACGAAG‑CGGGATTACCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCGGTCGa                                                                         <  2:424570/231‑1 (MQ=255)                                                                                                                                                                     aaaCGCGCGCAGGTACGCCCGCTGGCGGTTTGCGTACATAAAGATCG‑CCAGTTGCT‑GCAATTTTGCTACAACCCAGAAGTGGCCG‑ACAGAGCGCTG‑AAAATCCTCGACGTAG‑CGGGATTACCCGGGGAAATGCGCCGCGGT‑CCGGGGCGGGCGGCGGGGGGGATGGGTGGGGCGAGGGGCACCCTGAGACCGCCGCACTGTCCCCGCTCCTGGGAGGAACCGAAAAGGCAGACgcggcg                                                                         >  1:217478/1‑234 (MQ=255)                                                                                                                                                                      acccgtcgccgtatcGCCCTCTGGGGTTTTGCTTACAGTACGAGCG‑CCAGGTGCT‑GCAATTTTGCGACCCCTCAGAAGTGGCCG‑ACAGTGCGCTG‑AAAATCCTCGACGAAG‑CGGGATTACCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCGGTCGa                                                                         <  1:381264/224‑1 (MQ=255)                                                                                                                                                                          cgcgcAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCG‑CCAGTTGCT‑GCAATTTTGCTACACCTCAGAAGTGGCCG‑ACAGTGCGCTG‑AAAATCCTCGACGAAG‑CGGGATTACCTGGCGAACTGCGCCTTCGT‑CAGGGGCTGGCGCTGGGGGCGATGGTCGGGGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCGGTCGAAttt                                                                     >  2:232821/1‑238 (MQ=255)                                                                                                                                                                             gcAGGTACGCCCGGTGGGGGTTGGCGTACATAACGATCG‑CCAGTTGCT‑GCAATTTTGCTACACCGCAGAAGTGGCCG‑ACAGTGCGCTG‑AAAATCCTCGACGAAG‑CGGGATTACCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCGGTCGAATTTAcc                                                                  <  2:875877/238‑1 (MQ=255)                                                                                                                                                                             gcAGGTACGCCCGGTGGCGGTTGGCGTACATAACGGTCG‑CCAGTTGCT‑GCAATTTTGCTACCCCTCAGAAGTGGCCG‑ACAGGTCGCTG‑AAAATCCTCGACGAAG‑CGGGATTACCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCGGTCGAATTTAcc                                                                  <  2:36048/238‑1 (MQ=255)                                                                                                                                                                                  tgcGCCCGCTGGCGGTTGGCGTACATAACGATCG‑CCAGTTGCT‑GCAATTTTGCTACACCTCAGAAGTGGCCG‑ACAGTGCGCTG‑AAAATCCTCGACGAAG‑CGGGATTACCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCGGTCGAATTTACCTgg                                                               <  2:3406/234‑1 (MQ=255)                                                                                                                                                                                          tgggccgtttgctgtattAACGATCG‑CCCGTTCGC‑GCAATTTTGCTACTCCTCAGAAGGGGCCG‑ACGTGGCGCTG‑AAAATCCTCGGCGAAG‑CGGGATTCCCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCGGTCGAATTTACCTGGCAGTCCGATGa                                                    <  1:93600/222‑1 (MQ=255)                                                                                                                                                                                                     gttACATACCGATCG‑CCAGTTGCT‑GCAATTTGGCTACACCTCGGAAGTGGCCG‑ACAGTGCGCTG‑AAACTCCTCGACGAAG‑CGGGATTACCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCGGTCGAATTTACCTGGCAGTCCGATGACGGCATCAGc                                          <  2:72181/236‑1 (MQ=255)                                                                                                                                                                                                      cttcATAACGTCCG‑CCTGTTGCC‑TCAATTTTGCCTCCCCTCAGAAGGGGCCG‑ACAGTGCGCTG‑AAAATCCTCGACGGAG‑CGGGATTACCTGGCGAACTGCGCCGGCGT‑CAGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCGGTCGAATTTACCTGGCAGTCCGATGACGGCATCAGCCt                                        <  1:273826/236‑1 (MQ=255)                                                                                                                                                                                                       tacagcaCGATCG‑CCAGTTGCT‑GCAATTTTGCTACCCCTCAGAAGGGGCCG‑ACAGTGCGCTG‑AAAATCCTCGACGAAG‑CGGGATTACCTGGCGAACTGCGCCTGCGT‑CAGGGGCGGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCGGTCGAATTTACCTGGCAGTCCGATGACGGCATCAGCCt                                        <  2:288660/232‑1 (MQ=255)                                                                                                                                                                                                       tACATAACGATCG‑CCAGTTGCT‑GCAATTTTGCTACCCCTCAGAAGTGGCCG‑ACAGTGCGCTG‑AAAATCCTCGACGAAG‑CGGGATTACCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCGGTCGAATTTACCTGGCAGTCCGATGACGGCATCAGCCt                                        <  2:549461/238‑1 (MQ=255)                                                                                                                                                                                                       tACATAACGATCG‑CCAGTTGCT‑GCAATTTTGCTACACCTCAGAAGTGGCCG‑ACAGTGCGCTG‑AAAATCCTCGACGAAG‑CGGGATTACCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCGGTCGAATTTACCTGGCAGTCCGATGACGGCATCAGCCtg                                       >  2:613425/1‑239 (MQ=255)                                                                                                                                                                                                        aCATAACGTTCG‑CAAGTTGCT‑GCAATTTTGCTACACCTAAGAAGTGGCCG‑ACAGGGCGCTG‑AAAATCCTCGACGAAG‑CGGGATTACCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCGGTCGAATTTACCTGGCAGTCCGATGACGGCATCAGCCtgg                                      <  2:52838/239‑1 (MQ=255)                                                                                                                                                                                                              cGATCG‑CCAGTTGCT‑GCAATTTTGCTACACCTCAGAAGTGGCCG‑ACAGTGCGCTG‑AAAATCCTCGACGAAG‑CGGGATTACCTGGCGAACTGCGCCTGCGT‑CATGGGCTGGCGCTGGTGGCGCTGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAAGCCACCCGGTCGCATTTTCCTGGCCGTCCGCTGACGCCCTCAcccctgggtgag                                >  1:849162/1‑230 (MQ=255)                                                                                                                                                                                                              cGATCG‑CCAGTTGCT‑GCAATTTTGCTACACCTCAGAAGTGGCCG‑ACAGTGCGCTG‑AAAATCCTCGACGAAG‑CGGGATTACCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCa                                                                                                             >  1:122824/1‑162 (MQ=255)                                                                                                                                                                                                              cGATCG‑CCAGTTGCT‑GCAATTTTGCTACACCTCAGAAGTGGCCG‑ACAGTGCGCTG‑AAAATCCTCGACGAAG‑CGGGATTACCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCa                                                                                                             <  2:122824/162‑1 (MQ=255)                                                                                                                                                                                                              cGATCG‑CCAGTTGCT‑GCAATTTTGCTACACCTCAGAAGTGGCCG‑ACAGTGCGCTG‑AAAATCCTCGACGAAG‑CGGGATTACCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCGGTCGAATTTACCTGGCAGTCGGATGACGGCATCAGCCtgtgggcgg                                >  2:617426/1‑232 (MQ=255)                                                                                                                                                                                                                 tCG‑CCAGTTGCT‑GCAATTTTGCTACACCTCAGAAGTGGCCG‑ACAGCGCGCTG‑AAAATCCTCGACGACG‑CGGGCTTACCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCCCTGGTTGCCATGGTCGGGGCAGGCGTGACCCGTAACCCGCTGCATTGCCGCCGCTTCTGGCAGCCACTTGAAGGCCAGCCGGGCGAGTTTACCCGGCGTGCCGATGACGGCCTAAGCCtggtgggcgaat                             >  1:544432/1‑233 (MQ=255)                                                                                                                                                                                                                 gtcccaggTTGCT‑GCAATTTTGCTACACCTCAGAAGGGCCCG‑ACAGTGCGCTG‑AAAATCCTCGGCGAAG‑CGGGATTCCCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCGGTCGAATTTACCTGGCAGTCCGATGACGGCATCAGCCTGGTGGCAGTa                              <  1:521471/232‑1 (MQ=255)                                                                                                                                                                                                                  cG‑CCAGTTGCT‑GCAATTTTGCTACACCTCAGAAGTGGCCG‑ACAGTGCGCTG‑AAAATCCTCGACGAAG‑CGGGATTACCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCCCCCTTAACCCGCTGCATTGCCACCGCTTCTGGAAGCAACTGAAAGGCCAGCCGGCCGAATTTACCTGGCATGCCGATGAcaggaacacccgggggcgcgtacg                            >  1:691267/1‑216 (MQ=255)                                                                                                                                                                                                                     ccAGTTGCT‑GCAATTTTGCTACACCTCAGAAGTGGCCG‑ACAGTGCGCTG‑AAAATCCTCGACGAAG‑CGGGATTACCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCGGTCGAATTTACCTGGCAGCCCGATTACGGCCATAGCCTGGTGGCAGTACTg                           >  2:308241/1‑238 (MQ=255)                                                                                                                                                                                                                            ct‑gcAATTTTGCTACACCTCAGAAGTGGCCG‑ACAGTGCGCTG‑AAAATCCTCGACGAAG‑CGGGATTACCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGACAGGCCAGCCGGTCGAATTTACCTGGCAGTCCGATGACGGCATCAGCCTGGTGGCAGTACTGCGCACCg                    >  2:428715/1‑238 (MQ=255)                                                                                                                                                                                                                                      ggCTACCCCTCAGAAGGGGCCG‑GCAGTGCGCTG‑AAAATCCTCGACGAAG‑CGGGATTACCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCGGTCGAATTTACCTGGCAGTCCGATGACGGCATCAGCCTGGTGGCAGTACTGCGCACCGGCccgaccga          <  2:379580/238‑1 (MQ=255)                                                                                                                                                                                                                                              cTCAGAAGTGGCCG‑ACAGTGCGCTG‑AAAATCCTCGACGAAG‑CGGGATTACCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCGGTCGAATTTACCTGGCAGTCCGATGAc                                                   <  2:890410/190‑1 (MQ=255)                                                                                                                                                                                                                                              cTCAGAAGTGGCCG‑ACAGTGCGCTG‑AAAATCCTCGACGAAG‑CGGGATTACCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCGGGCGCATTTACCTGGCAGTCCGATGCCGGCATACGCCTGGTGGCAGTACTGCGAAACGGcccacacgaaaacccga  >  1:25618/1‑225 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                  TAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCC‑GGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCG‑CCAGTTGCT‑GCAATTTTGCTACACCTCAGAAGTGGCCG‑ACAGTGCGCTG‑AAAATCCTCGACGAAG‑CGGGATTACCTGGCGAACTGCGCCTGCGT‑CAGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCGGTCGAATTTACCTGGCAGTCCGATGACGGCATCAGCCTGGTGGCAGTACTGCGCACCGGCCCGACCGAAAGCCTGA  >  NC_000913/4130694‑4131166

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