Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	344,674	T→C	Y167H (TAC→CAC)	yahL →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	344,674	0	T	C	100.0%	133.5 / NA	38	Y167H (TAC→CAC)	yahL	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/20);  total (18/20)

CAGGATGATGCGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGAT  >  NC_000913/344554‑344800                                                                                                                         |                                                                                                                               cAGGATGATGCGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTga                                                                                                              >  1:1158546/1‑139 (MQ=255)        atgCGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAAtttt                                                                                                       >  1:386399/1‑139 (MQ=255)         tgCGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTc                                                                                                      >  2:417106/1‑139 (MQ=255)                   gATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGa                                                                                            <  2:240176/139‑1 (MQ=255)                    atatTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGaa                                                                                           <  2:1158546/139‑1 (MQ=255)                           aGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAaga                                                                                        >  2:460788/1‑135 (MQ=255)                           aGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAaga                                                                                        <  1:460788/135‑1 (MQ=255)                                aCGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCa                                                                               >  2:70405/1‑139 (MQ=255)                                            gACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGaaa                                                                                          <  1:357963/116‑1 (MQ=255)                                            gACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGaaa                                                                                          >  2:357963/1‑116 (MQ=255)                                                  ttACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCt                                                                                                                <  1:130799/88‑1 (MQ=255)                                                  ttACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCt                                                                                                                >  2:130799/1‑88 (MQ=255)                                                        ggCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAAt                                                       <  2:386399/139‑1 (MQ=255)                                                         gCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATg                                                      <  1:70405/139‑1 (MQ=255)                                                                   ttAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTc                                            >  2:892183/1‑139 (MQ=255)                                                                       tGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGaa                                                         <  2:719806/122‑1 (MQ=255)                                                                       tGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGaa                                                         >  1:719806/1‑122 (MQ=255)                                                                               aTCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGAcg                                <  2:1138612/139‑1 (MQ=255)                                                                                    ttCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGccc                           >  1:1230421/1‑139 (MQ=255)                                                                                      cAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAAt                                                       >  1:139318/1‑109 (MQ=255)                                                                                      cAGTGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAAt                                                       <  2:139318/109‑1 (MQ=255)                                                                                         tGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGa                                                                                            <  1:242148/69‑1 (MQ=255)                                                                                         tGCGAACAGATTGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGa                                                                                            >  2:242148/1‑69 (MQ=255)                                                                                                   ttGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAaga                                                                                        >  1:1094673/1‑63 (MQ=255)                                                                                                   ttGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAaga                                                                                        <  2:1094673/63‑1 (MQ=255)                                                                                                   ttGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATgtcgt                                                  <  2:983198/101‑1 (MQ=255)                                                                                                   ttGCAAATGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGAAAAAATGGCTTCGCGAAATgtcgt                                                  >  1:983198/1‑101 (MQ=255)                                                                                                       aaaTGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAag                                                                                         <  2:11821/58‑1 (MQ=255)                                                                                                       aaaTGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAag                                                                                         >  1:11821/1‑58 (MQ=255)                                                                                                        aaTGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTGTGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCa       <  1:981374/139‑1 (MQ=255)                                                                                                         aTGTATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCAt      <  1:892183/139‑1 (MQ=255)                                                                                                            tATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATgtcgtc                                                 >  2:174116/1‑93 (MQ=255)                                                                                                            tATTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATgtcgtc                                                 <  1:174116/93‑1 (MQ=255)                                                                                                             aTTCTTCTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGat  <  2:1230421/139‑1 (MQ=255)                                                                                                              ttcttcTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCt                                                >  1:649336/1‑92 (MQ=255)                                                                                                              ttcttcTTGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCt                                                <  2:649336/92‑1 (MQ=255)                                                                                                                     tGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATgg    >  2:1165764/1‑129 (MQ=255)                                                                                                                     tGAACACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATgg    <  1:1165764/129‑1 (MQ=255)                                                                                                                         |                                                                                                                               CAGGATGATGCGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGAT  >  NC_000913/344554‑344800

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