Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,874,078	C→T	*35* (TAG→TAA)	yoaI ←	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,874,078	0	C	T	100.0%	118.6 / NA	37	*35* (TAG→TAA)	yoaI	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/21);  total (16/21)

CTCTTCGAACGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAACTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACGTGATAATCAGTGTCTCGACAAACATTTGATCGTTCATAGCATGCCCCCGGGAGATAGCTGTTTTCAC  >  NC_000913/1873952‑1874212                                                                                                                               |                                                                                                                                       ctctTCGAACGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGt                                                                                                                            <  1:754398/139‑1 (MQ=255)     tCGAACGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGa                                                                                                                        >  1:909954/1‑139 (MQ=255)        aaCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATaa                                                                                                                     <  2:222055/139‑1 (MQ=255)                 atgaAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACgg                                                                                                            <  2:674366/139‑1 (MQ=255)                      aTCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACaaa                                                                                                       <  1:50609/139‑1 (MQ=255)                       tCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTAGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACaaaa                                                                                                      >  1:654791/1‑139 (MQ=255)                           cTTGAATGCAGCCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAg                                                                                                  <  2:654791/139‑1 (MQ=255)                                           cTGGATAAATCCGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGCCAAAACCAGTACAACAGCAATGGCa                                                                                  >  1:290430/1‑139 (MQ=255)                                           cTGGATAAATCCGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCa                                                                                  >  1:360977/1‑139 (MQ=255)                                             ggATAAATCCGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCaaa                                                                                <  2:909954/139‑1 (MQ=255)                                                     ccGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATg                                                                        >  2:121352/1‑139 (MQ=255)                                                       gTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGAc                                                                      <  2:360977/139‑1 (MQ=255)                                                        ttATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACg                                                                     >  1:873452/1‑139 (MQ=255)                                                        ttATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTACCGCCGCACTGCGTAATTAGCCGGTCCGTACGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACg                                                                     >  1:664532/1‑139 (MQ=255)                                                                     cgcgGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACGTGATAATCAgtgt                                                        >  1:626655/1‑139 (MQ=255)                                                                       cgGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATaag                                                                                                                    >  2:889272/1‑77 (MQ=255)                                                                       cgGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATaag                                                                                                                    <  1:889272/77‑1 (MQ=255)                                                                         gATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACGTGATAATCAGTGTCTCg                                                    <  2:321361/139‑1 (MQ=255)                                                                           ttattTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACGTGATAATCAGTGTCTCGAc                                                  <  1:794842/139‑1 (MQ=255)                                                                            tattTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCa                                                                                        <  1:224507/100‑1 (MQ=255)                                                                            tattTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCa                                                                                        >  2:224507/1‑100 (MQ=255)                                                                                      ttgcgcaACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACGTGATAATCAGTGTCTCGACAAACATTTGa                                        <  1:758849/138‑1 (MQ=255)                                                                                      ttgcgcaACCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACGTGATAATCAGTGTCTCGACAAACATTTGa                                        >  2:758849/1‑138 (MQ=255)                                                                                             aCCAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACGTGATAATCAGTGTCTCGATAAACATTTGATCGTTCAt                                <  2:664532/139‑1 (MQ=255)                                                                                              ccAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACGTGATAATCAGTGTCTCGAc                                                  <  1:853030/120‑1 (MQ=255)                                                                                              ccAGAACTGGTCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACGTGATAATCAGTGTCTCGAc                                                  >  2:853030/1‑120 (MQ=255)                                                                                                        tCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGAc                                                                                                          >  1:933835/1‑54 (MQ=255)                                                                                                        tCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGAc                                                                                                          <  2:933835/54‑1 (MQ=255)                                                                                                        tCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACGTGATAATCAGTGTCTCGACAAACATTTGATCGTTc                                  >  2:431442/1‑126 (MQ=255)                                                                                                        tCAATTTCCGCCGCACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACGTGATAATCAGTGTCTCGACAAACATTTGATCGTTc                                  <  1:431442/126‑1 (MQ=255)                                                                                                                    gcACTGCGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACGTGATAATCAGTGTCTCGACAAACATTTGATCGTTCATAGCATGCCCCCGGGAGATAGCTg         <  2:588980/139‑1 (MQ=255)                                                                                                                          cGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAAcc                                                                                                    >  1:90899/1‑42 (MQ=255)                                                                                                                          cGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAAcc                                                                                                    <  2:90899/42‑1 (MQ=255)                                                                                                                          cGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAAc                                                                            <  2:189577/66‑1 (MQ=255)                                                                                                                          cGTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAAc                                                                            >  1:189577/1‑66 (MQ=255)                                                                                                                           gTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACGTGATAATCAGTGTCTCGACAAACATTTGATCGTTCATAGCATGCCCCCGGGAGATAGCTGTTTTCAc  <  2:290430/139‑1 (MQ=255)                                                                                                                           gTAATTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACGTGATAATCAGTGTCTCGACAAACATTTGATCGTTCATAGCATGCCCCCGGGAGATAGCTGTTTTCAc  <  2:626655/139‑1 (MQ=255)                                                                                                                               |                                                                                                                                       CTCTTCGAACGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGAACTGGTCAATTTCCGCCGCACTGCGTAACTAGCCGGTCCGTTCGATAAGAAGAACGGACAAAACCAGTACAACAGCAATGGCAAAAAACGATGACGTGATAATCAGTGTCTCGACAAACATTTGATCGTTCATAGCATGCCCCCGGGAGATAGCTGTTTTCAC  >  NC_000913/1873952‑1874212

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