Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,226,817	A→G	Y73C (TAC→TGC)	ycgJ →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,226,817	0	A	G	100.0%	144.8 / NA	39	Y73C (TAC→TGC)	ycgJ	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/21);  total (18/21)

GACGCTGGCGGACTTCAATCTCCGGCGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCG  >  NC_000913/1226675‑1226961                                                                                                                                               |                                                                                                                                                 gACGCTGGCGGACTTCAATCTCCGGCGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTGCGATcc                                                                                                                                            >  1:140327/1‑149 (MQ=255)         cGGACTTCAATCTCCGGCGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTGCGATCCCAGCGAAttt                                                                                                                                  <  1:237751/151‑1 (MQ=255)                  aTCTCCGGCGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTGCGATCCCAGCGAATTTACATTCGCt                                                                                                                         >  2:183782/1‑151 (MQ=255)                      ccGGCGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTGCGATCCCAGCGAATTTACATTCGCTAACg                                                                                                                     <  1:556977/151‑1 (MQ=255)                                   ttgttTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTGCGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCg                                                                                                         >  2:161595/1‑150 (MQ=255)                                       ttGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTGCGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGtt                                                                                                    >  1:38644/1‑151 (MQ=255)                                           gACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTGCGATCCCAGCGAATTTACAtt                                                                                                                             >  2:397076/1‑122 (MQ=255)                                           gACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTGCGATCCCAGCGAATTTACAtt                                                                                                                             <  1:397072/122‑1 (MQ=255)                                              ccGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTGCGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGaaa                                                                                             <  2:488935/151‑1 (MQ=255)                                                 tataTTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTGCGATCCCAGCGAATTTACATTCGCTAACGGTATATTTTGCGATGTTAAAGAAAAAt                                                                                          >  1:194538/1‑151 (MQ=255)                                                          gtgAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTGCGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGtg                                                                                 <  1:86874/151‑1 (MQ=255)                                                            gAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTGCGATCCCAGCGAATTTACACTCGCTAACGGTGTATTTTGCGATGTTAAACAAAAATTATGTCGtgat                                                                               <  1:222357/151‑1 (MQ=255)                                                             aaCTCAGATGGCATTTCCCCAGAGTTAATAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTGCGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGtgatg                                                                              <  2:447318/151‑1 (MQ=255)                                                              aCTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTGCGATCCCAGCGAATTTACATTCGCTAACGGTGTAtttt                                                                                                            <  1:156935/120‑1 (MQ=255)                                                              aCTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTGCGATCCCAGCGAATTTACATTCGCTAACGGTGTAtttt                                                                                                            >  2:156938/1‑120 (MQ=255)                                                              aCTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTGCGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGtgatga                                                                             <  2:399290/151‑1 (MQ=255)                                                               cTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTGCGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTgatgat                                                                            >  2:247873/1‑151 (MQ=255)                                                                  aGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTGCGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGt                                                                         >  1:42306/1‑151 (MQ=255)                                                                   gATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTGCGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGtt                                                                        <  2:276541/151‑1 (MQ=255)                                                                       gCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTGCGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTAttt                                                                    <  2:223569/151‑1 (MQ=255)                                                                         aTTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTGCGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTg                                                                  <  2:103866/151‑1 (MQ=255)                                                                         aTTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTGCGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTg                                                                  >  2:132014/1‑151 (MQ=255)                                                                            tCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTGCGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGgt                                                                >  2:469676/1‑150 (MQ=255)                                                                             ccccAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTGCGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGgtgt                                                              <  2:620637/151‑1 (MQ=255)                                                                              cccAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTGCGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGtgtg                                                             >  2:430809/1‑151 (MQ=255)                                                                                   agTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTGCGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATgg                                                        <  2:650184/151‑1 (MQ=255)                                                                                         cAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTGCGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCg                                                  >  1:342794/1‑151 (MQ=255)                                                                                         cAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTGCGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCg                                                  >  1:328255/1‑151 (MQ=255)                                                                                          aaGGAAATATCTCGGCGAAAAAGCCGCTGAAATCTTACAATCATTACAAGGCTGCGATCCCAGCGAATTTACATTCGGTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCg                                                  <  1:382270/150‑1 (MQ=255)                                                                                             gAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTGCGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGt                                              <  2:617619/151‑1 (MQ=255)                                                                                                   tctcGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTGCGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGaaa                                         >  2:538553/1‑150 (MQ=255)                                                                                                       ggtgAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTGCGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAaatca                                    <  2:338085/148‑1 (MQ=255)                                                                                                                      gAAAACTTACAATCATTACAAGGCTGCGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACcaca                      >  1:422582/1‑150 (MQ=255)                                                                                                                         aaCTTACAATCATTACAAGGCTGCGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACaacg                   >  1:473452/1‑148 (MQ=255)                                                                                                                                   cATTACAAGGCTGCGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAATCACCACAAAGATGttatttat        <  1:262850/151‑1 (MQ=255)                                                                                                                                       aCAAGGCTGCGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATgtgt    >  1:237675/1‑151 (MQ=255)                                                                                                                                        cAAGGCTGCGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTc   <  1:433274/151‑1 (MQ=255)                                                                                                                                        cAAGGCTGCGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTc   <  2:623692/151‑1 (MQ=255)                                                                                                                                         aaGGCTGCGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCg  <  1:316007/151‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                 GACGCTGGCGGACTTCAATCTCCGGCGTGCGGCGTTGTTTGCGACCCGTATATTTGTGTGAACTCAGATGGCATTTCCCCAGAGTTAACAAGGAAATATCTCGGCGAAAAAGCCGCTGAAAACTTACAATCATTACAAGGCTACGATCCCAGCGAATTTACATTCGCTAACGGTGTATTTTGCGATGTTAAAGAAAAATTATGTCGTGATGATCGTTATTTTGGTGTGGATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCG  >  NC_000913/1226675‑1226961

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