Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,757,772	A→C	L19W (TTG→TGG)	yiaY ←	L‑threonine dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,757,772	0	A	C	100.0%	126.9 / NA	36	L19W (TTG→TGG)	yiaY	L‑threonine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (17/19);  total (17/19)

GCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAG  >  NC_000913/3757642‑3757896                                                                                                                                   |                                                                                                                             gCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCCAt                                                                                                                            >  1:510682/1‑133 (MQ=255)  cTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCCATg                                                                                                                           >  1:449598/1‑133 (MQ=255)     aaaaTATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAACTCCATAATCTGCAATCCTATTCATTACATCGGTCCATGAAt                                                                                                                        >  2:201372/1‑133 (MQ=255)     aaaaTATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCCATGAAt                                                                                                                        <  2:597606/133‑1 (MQ=255)               gTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCCGCCATCATCTTCCTTGAATCAGTCCATGAATGAGCGCCGAt                                                                                                              >  2:295464/1‑133 (MQ=255)                  tttcCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCCATGAATCAGCGCCGATGAc                                                                                                           <  2:4808/132‑1 (MQ=255)                          tcTTTTTGCCCATCGTCCGCCATACCTAGTTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCCATGAATCAGCGCCGATGACATTCACAg                                                                                                   <  2:299744/132‑1 (MQ=255)                          gCTTTTTGCACTCCCCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCCATGAATCAGCGCCGATGACATTCACAg                                                                                                   <  2:450493/133‑1 (MQ=255)                           tttttgGCACCTCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCCATGAATCAGCGCCGATGACATTCACAGa                                                                                                  <  2:553654/132‑1 (MQ=255)                              tttGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCCATGAATCAGCGCCGATGACATTCACAg                                                                                                   >  2:497959/1‑129 (MQ=255)                              tttGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCCATGAATCAGCGCCGATGACATTCACAGAAgg                                                                                               >  1:44376/1‑133 (MQ=255)                              tttCCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCCATGAATCAGCGCCGATGACATTCACAg                                                                                                   <  1:497959/129‑1 (MQ=255)                                     aTCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCCATGAATCAGCGCCGATGACATTCACAGAAgg                                                                                               <  2:118023/126‑1 (MQ=255)                                     aTCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCCATGAATCAGCGCCGATGACATTCACAGAAgg                                                                                               >  1:118023/1‑126 (MQ=255)                                         cccgccATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCCATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACg                                                                                    <  2:549375/133‑1 (MQ=255)                                                  ccTAATTTCGGTAACCTATTGTCAGTGATAATAAAGGTACGTGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCCATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTg                                                                           <  2:252240/133‑1 (MQ=255)                                                  ccTAATTTCGCTAACTTATTGTCAATGAAAATAAAGGTACGTGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCCATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAGGCTg                                                                           <  2:196491/133‑1 (MQ=255)                                                          cGTTAACACATTGTCGGTGACACTTAAGGTTCGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCCATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAaca                                                                   <  1:215530/133‑1 (MQ=255)                                                            ttAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCCATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAacaca                                                                 >  2:460013/1‑133 (MQ=255)                                                             tAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCCATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAacacac                                                                >  1:50894/1‑133 (MQ=255)                                                             tAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCCATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAacacac                                                                <  1:460013/133‑1 (MQ=255)                                                                                      gTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCCATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAgggcggg                                       <  2:449598/133‑1 (MQ=255)                                                                                         cGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCCATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGGAAAAAAGGGCGGGCAg                                    >  1:421711/1‑133 (MQ=255)                                                                                           ggTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCCATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACAATCCTGAAAGTGGAAAAAAGGGCGGGCAGCt                                  >  1:237129/1‑133 (MQ=255)                                                                                                     tAATCTGCCATCATATTCATTGCATCAGTCCATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGAt                        >  1:622561/1‑133 (MQ=255)                                                                                                      aaTCTGCCATCATATTCATTGCATCAGTCCATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGAtt                       >  1:48724/1‑133 (MQ=255)                                                                                                            ccATCATATTCATTGCATCAGTCCATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCAtt                 >  2:86995/1‑133 (MQ=255)                                                                                                              atcatATTCATTGCATCAGTCCATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCg               <  2:313599/133‑1 (MQ=255)                                                                                                              atcatATTCATTGCATCAGTCCATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGGAAAAAAGGGCGCGCAGCTTACTGAGGATTTTCATTCg               <  2:421711/133‑1 (MQ=255)                                                                                                                   atTCATTGCATCAGTCCATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAat          <  2:1608/133‑1 (MQ=255)                                                                                                                      cATTGCATCAGTCCATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATg       <  1:86995/133‑1 (MQ=255)                                                                                                                      cATTGCATCAGTCCATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATg       <  2:514701/133‑1 (MQ=255)                                                                                                                        ttGCATCAGTCCATGAATAAGCTCCGATGATATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGAt     <  1:197386/133‑1 (MQ=255)                                                                                                                         tGCATCAGTCCATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAACGGCGGGCAACTTACTGAGGATTTTCATTCGGAAATATGATc    >  2:197386/1‑133 (MQ=255)                                                                                                                          gCATCAGTCCATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGCGAAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCa   >  1:361564/1‑133 (MQ=255)                                                                                                                           cATCAGTCCATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAg  >  1:112903/1‑133 (MQ=255)                                                                                                                                   |                                                                                                                             GCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAG  >  NC_000913/3757642‑3757896

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