Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,743,164	C→T	E27K (GAG→AAG)	ribC ←	riboflavin synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,743,164	0	C	T	100.0%	86.4 / NA	27	E27K (GAG→AAG)	ribC	riboflavin synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (10/17);  total (10/17)

TGGTCATAATATGACCTGACATTAAGTGTCCGCCAATTTCATCACTGAATTTCGCCGCACGCTCAACGTTTACCCAATCCCCCACTTTTAAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGAAAACGCCCTTACCAGTGAAGGGGTAAGAATGGCTATTTTTTCACTGGAGAATTAATAAATCCTCGCTACAATAGACTGAATTTCCCCTGCTTCTTCTTTTTGCT  >  NC_000913/1742927‑1743401                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                              tGGTCATAATATGACCTGACATTAAGTGTCCGCCAATTTCATCACTGAATTTCGCCGCACGCTCAACGTTTACCCAATCCCCCACTTTTAAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTTCACCACATGCGt                                                                                                                                                                                                                                   >  1:461751/1‑250 (MQ=255)     cATAATATGACCTGACATTAAGTGTCCGCCAATTTCATCACTGAATTTCGCCGCACGCTCAACGTTTACCCAATCCCCCACTTTTAAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTTCACCACATGCGTACGa                                                                                                                                                                                                                               <  2:104112/250‑1 (MQ=255)             gACCTGACATTAAGTGTCCGCCAATTTCATCACTGAATTTCGCCGCACGCTCAACGTTTACCCAATCCCCCACTTTTAAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTTCACCACATGCGTACGAAAATTTGGt                                                                                                                                                                                                                      >  2:509194/1‑251 (MQ=255)                     aTTAAGTGTCCGCCAATTTCATCACTGAATTTCGCCGCACGCTCAACGTTTACCCAATCCCCCACTTTTAAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTTCACCACATGCGTACGAAAATTTGGTTTCTCGTc                                                                                                                                                                                                              >  2:460663/1‑251 (MQ=255)                         aGTGTCCGCCAATTTCATCACTGAATTTCGCCGCACGCTCAACGTTTACCCAATCCCCCACTTTTAAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTTCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATc                                                                                                                                                                                                          <  1:509194/251‑1 (MQ=255)                             tCCGCCAATTTCATCACTGAATTTCGCCGCACGCTCAACGTTTACCCAATCCCCCACTTTTAAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGACGTCCAGCATGTGGTCGGGTAACTTCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGac                                                                                                                                                                                                       >  1:528043/1‑250 (MQ=255)                               cGCCAATTTCATCACTGAATTTCGCCGCACGCTCAACGTTTACCCAATCCCCCACTTTTAAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTTCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACAcc                                                                                                                                                                                                    <  1:475986/251‑1 (MQ=255)                                           caCTGAATTTCGCCGCTCGCTCAACGTTTACCCAATCCCCCACTTTTAAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGGATCTTCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATGGACACCAGTTTTGCGGTg                                                                                                                                                                                        <  2:511129/251‑1 (MQ=255)                                                         gcACGCTCAACGTTTACCCAATCCCCCACTTTTAAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTGCCGTCACGGTCAGGCAGCATCCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTTCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATccc                                                                                                                                                                          <  1:447168/251‑1 (MQ=255)                                                                 aaCGTTTACCCAATCCCCCACTTTTAAATCGACAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTTCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACa                                                                                                                                                                  <  2:131370/251‑1 (MQ=255)                                                                  aCGTTTACCCAATCCCCCATTGTTAAATCACCAAGATTGGCAATGCGTAACTTTTCTTTCAACAACGCACAACTGACATGGTTGCCGTTAATTTCCGTCACGTTCAGGCAGCTACCGTTATGCGCCACCGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTTCACCTCATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGTTTCCCTGTACAATCCCCGTAAACAt                                                                                                                                                                 <  2:481495/251‑1 (MQ=255)                                                                         cccAATCCCCCACTTTTAAATCGCCAAGATTGGTAATGCGTAACGTTTCTTGCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCCCCGGTTTCCAGGCCGGCCAGCATGTGGTCGGGTAACTTCACCACATGCGTACGAAAAGTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACa                                                                                                                                                          <  2:529783/251‑1 (MQ=255)                                                                          ccAATCCCCCACTTTTAAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTTCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACaa                                                                                                                                                         <  2:427775/251‑1 (MQ=255)                                                                                            aTCGCCACGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTTCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTaa                                                                                                                                        >  2:375478/1‑250 (MQ=255)                                                                                            aTCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTTCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAg                                                                                                                                       <  1:460663/251‑1 (MQ=255)                                                                                               gCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTTCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACa                                                                                                                                    <  1:192991/251‑1 (MQ=255)                                                                                                   aGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTTCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCt                                                                                                                                >  2:150307/1‑251 (MQ=255)                                                                                                     aTTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTTCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGt                                                                                                                              <  1:150307/251‑1 (MQ=255)                                                                                                          tAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTTCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGc                                                                                                                         >  1:424387/1‑251 (MQ=255)                                                                                                           aattcGTAACGTTTCGTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTGCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTGTCCAGGCCGTCCAGCATGTGGTCGGGTAACTTCACCACAGGCGTACGAAAATTTGTTTTCTCGTCAATCGACGCCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACGCCGGAAAGCAGTTAAGACATTATGCTCAGca                                                                                                                        <  2:55667/247‑1 (MQ=255)                                                                                                            aTGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTTCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTTCAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGcac                                                                                                                       >  1:314402/1‑251 (MQ=255)                                                                                                                             catcaGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTTCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGAAAAc                                                                                                      <  1:70542/251‑1 (MQ=255)                                                                                                                                               aTGGTTTCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTTCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACGTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGAAAACGCCCTTACCAGTGAAggg                                                                                    <  1:509538/251‑1 (MQ=255)                                                                                                                                                                          cAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTTCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGAAAACGCCCTTACCAGTGAAGGGGTAAGAATGGCTATTTTTTCACTGga                                                          <  1:502930/250‑1 (MQ=255)                                                                                                                                                                                                  ggAAGCACCGGGTTCCAGGCAGTCCAGCATGTGGTCGGGTAACTTCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGAAAACGACCTTACCAGTGACGGGGCAAGAATGGCTATTTTTTCACTGGAGAATTAATAAATCCTCGCTACAAt                                  >  1:133334/1‑250 (MQ=255)                                                                                                                                                                                                        aCCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTTCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGAAAACGCCCTTACCAGTGAAGGGGTAAGAATGGCTATTTTTTCACTGGAGAATTAATAAATCCTCGCTACAATAGACTg                            >  2:425695/1‑250 (MQ=255)                                                                                                                                                                                                                                  ggTCGGGTAACTTCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTTCAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGAAAACGCCCTTACCAGTGAAGGGGTTAGAATGGCTATTTTTTCACTGGAGAATTAATAAATCCTCGCTACAATAGACTGAATTTCCCCTGCTTCTTCTTTTtgct  <  1:47997/250‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                              TGGTCATAATATGACCTGACATTAAGTGTCCGCCAATTTCATCACTGAATTTCGCCGCACGCTCAACGTTTACCCAATCCCCCACTTTTAAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGAAAACGCCCTTACCAGTGAAGGGGTAAGAATGGCTATTTTTTCACTGGAGAATTAATAAATCCTCGCTACAATAGACTGAATTTCCCCTGCTTCTTCTTTTTGCT  >  NC_000913/1742927‑1743401

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑