Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	225,004	G→T	E31* (GAA→TAA)	murQ →	possible sugar phosphate isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	225,004	0	G	T	100.0%	78.5 / NA	28	E31* (GAA→TAA)	murQ	possible sugar phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (16/12);  total (16/12)

GCCGAATGATACATTTGATTTTCAAAAAACTAGGGAGTGTAAGAAGTGATGGAAAATAGTACGACCGAAGCGCGTAATGAAGCGACGATGCATCTTGATGAAATGACTGTGGAAGAGGCTTTAATTACGATGAATAAAGAAGATCAGCAAGTCCCGTTAGCAGTTCGAAAGGCAATACCACAATTGACAAAAGTAATTAAAAAAACAATTGCACAGTATAAAAAGGGTGGACGATTGATTTATATCGGTGCAGGTACAAGTGGAAGGTTGGGTGTC  >  CP000730/224866‑225141                                                                                                                                           |                                                                                                                                          gCCGAATGATACATTTGATTTTCAAAAAACTAGGGAGTGTAAGAAGTGATGGAAAATAGTACGACCGAAGCGCGTAATGAAGCGACGATGCATCTTGATGAAATGACTGTGGAAGAGGCTTTAATTACGATGAataaataa                                                                                                                                         >  1:354461/1‑141 (MQ=255)      aTGATACATTTGATTTTCAAAAAACTAGGGAGTGTAAGAAGTGATGGAAAATAGTACGACCGAAGCGCGTAATGAAGCGACGATGCATCTTGATGAAATGACTGTGGAAGAGGCTTTAATTACGATGAATAAATAAGATCa                                                                                                                                    >  1:42586/1‑141 (MQ=255)        gATACATTTGATTTTCAAAAAACTAGGGAGTGTAAGAAGTGATGGAAAATAGTACGACCGAAGCGCATAATTAAGCTACGATGCAACTTGATGAAATGACTGTGGAAGAGGCTTTAATTACGATGAATAAATAAGATCAGc                                                                                                                                  <  1:363555/141‑1 (MQ=255)         aTACATTTGATTTTCAAAAAACTAGGGAGTGTAAGAAGTGATGGAAAATAGTACGACCGAAGCGCGTAATGAAGCGACGATGCATCTTGATGAAATGACTGTGGAAGAGGCTTTAATTACGATGAATAAATAAGATCAGCa                                                                                                                                 <  1:448641/141‑1 (MQ=255)         aTACATTTGATTTTCAAAAAACTAGGGAGTGTAAGAAGTGATGGAAAATAGTACGACCGAAGCGCGTAATGAAGCGACGATGCATCTTGATGAAATGACTGTGGAAGAGGCTTTAATTACGATGAATAAATAAGATCAGCa                                                                                                                                 <  1:398800/141‑1 (MQ=255)          tACATTTGATTTTCAAAAAACTAGGGAGTGTAAGAAGTGATGGAAAATAGTACGACCGAAGCGCGTAATGAAGCGACGATGCATCTTGATGAAATGACTGTGGAAGAGGCTTTAATTACGATGAATAAATAAGATCAGCaa                                                                                                                                <  1:467190/141‑1 (MQ=255)                  aTTTTCAAAAAACTAGGGAGTGTAAGAAGTGATGGAAAATAGTACGACCGAAGCGCGTAATGAAGCGACGATGCATCTTGATGAAATGACTGTGGAAGAGGCTTTAATTACGATGAAAAAATAAGACCGGCAAGGCCCGtt                                                                                                                        >  1:172683/1‑141 (MQ=255)                     ttCAAAAACCTTGGGAGTGTAAGAAGTGATGGAAAATAGTACGACCGAAGCGCGTAATGAAGCGACGATGCATCTTGATGAAATGACTGTGGAAGAGGCTTTAATTACGATGAATAAATAAGATCAGCAAGTCCCGTTAGc                                                                                                                     <  1:198957/141‑1 (MQ=255)                          aaaaCTAGGGAGTGTAAGAAGTGATGGAAAATAGTACGACCGAAGCGCGTAATGAAGCGACGATGCATCTTGATGAAATGACTGTGGAAGCGGCATTAATTACGATGAATAAATAAGATCAGCAAGTCCCGTTAGCAGTTc                                                                                                                <  1:95750/141‑1 (MQ=255)                            aaCTAGGGAGTGTAAGAAGTGATGGAAAATAGTACGACCGAAGCGCGTAATGAAGCGACGATGCATCTTGATGAAATGACTGTGGAAGAGGCTTTAATTACGATGAATAAATAAGATCAGCAAGTCCCGTTAGCAGTTCGa                                                                                                              <  1:170792/141‑1 (MQ=255)                             aCTAGGGAGAGTAAGAAGTGATGGAAAATAGTACGACCGAAGCGCGTGATGAAGCGACGATGCATCTTGATGAGATGACTGTGGAAGAGGCTTTAATTACGATGAATAAATAAAATCAGCAAGTCCCGTTAGCAGTTCGaa                                                                                                             >  1:171456/1‑141 (MQ=255)                                       gtAAGAAGTGATGGAAAATAGTACGACCGAAGCGCGTAATGAAGCGACGATGCATCTTGATGAAATGACTGTGGAAGAGGCTTTAATTACGATGAATAAATAAGATCAGCAAGTCCCGTTAGCAGTTCGAAAGGCAATAcc                                                                                                   >  1:433878/1‑141 (MQ=255)                                            aagTGATGGAAAATAGTCAGAGACTAACGTTTAATGAAGCGACGATGCATCTTGATGAAATGACTGTGGACGAGGCTTTAATTACAATGAATAAATAAGATCAGCAAGTCCCGTTAGCAGTTCGAAAGGCAATACCACAAt                                                                                              >  1:118683/1‑141 (MQ=255)                                                  tGGAAAATAGTACGACCGAAGCGCGTAATGAGGCGACGATGCATCTTGATGAAATGACTGTGGAAGAGGCTTTAATTACGATGAATAAATAAGATCAGCAAGTCCCGTTAGCAGTTAGAAAGGCAATACCACAATTGACaa                                                                                        >  1:98640/1‑141 (MQ=255)                                                   ggAAAATAGTACGACCGAAGCGCGTAATGAAGCGACGATGCATCTTGATGAAATGACTGTGGAAGAGGCTTTAATTACGATGAATAAATAAGATCAGCAAGTCCCGTTAGCAGTTCGAAAGGCAATACCACAATTCACaaa                                                                                       >  1:414313/1‑141 (MQ=255)                                                                      gcgcGTAATGAAGCGACGATGCATCTTGATGAAATGACTGTGGAAGAGGCTTTAATTACGATGAATAAATAAGATCAGCAAGTCCCGTTAGCAGTTAGAAAGGCAATACCACAATTGACAAAAGTAATTAAAAAAACAAtt                                                                    >  1:467847/1‑141 (MQ=255)                                                                          gTAATGAAGCGACGATGCAACTTGATGAAATGACTGTGGAAGAGGCTTTAATTACGATGAATAAATAAGATCAGCAAGTCCCGTTAGCAGTTCGAAAGGCAATACCACAATTGACAAAAGTAATTAAAAAAACAATTGcac                                                                >  1:297125/1‑141 (MQ=255)                                                                             aTGAAGCGACGATGCATCTTGATGAAATGACTGTGGAAGAGGCTTTAATTACGATGAATAAATAAGATCAGCAAGTCCCGTTAGCAGTTCGACAGGCAATTCCTTAATTGACAAAAGTAATTAAAAAAACAATTGCACAGt                                                             <  1:217836/141‑1 (MQ=255)                                                                                    gacgaTGCATCTTGATGAAATGACTGTGGAAGAGGCTTTAATTACGATGAATAAATAAGATCAGCAAGTCCCGTTAGCAGTTCGAAAGGCAATACCACAATTGACAAAAGTAATTAAAAAAACAATTGCCAAGTCaaaaaa                                                      >  1:236596/1‑141 (MQ=255)                                                                                      cgaTGCATCTTGATGAAATGACTGTGGAAGAGGCTTTAATTACGATGAATAAATAAGATCAGCAAGTCCCGTTAGCAGTTCGAAAGGCAATACCACAATTGACAAAAGTAATTAAAAAAACAATTGCACAGTATAAAAAgg                                                    <  1:492179/141‑1 (MQ=255)                                                                                         tgcatttcgcggaAATGACTGTGGAAGAGTCTTTAATTACGATGAATCAATAAGACAAGCAAGTCACGTTAGCAGTGCGAAAGGCAATACCACAATTGACAAAAGTAATTAAAAAAACAATTGCACAGTATAAAAAGGGTg                                                 <  1:211690/130‑1 (MQ=255)                                                                                         tGCACCGTGATGAAATGACTGGGGAAGAGGCTTTAATTACGATGAATAAATAAGATCAGCAAGTCCCGTTAGCAGTTCGAAAGGCAAGACCACAATTGCCAAAAGTAATTAAAAAAACAATTGCACAGTATAAAAAGGGTg                                                 <  1:375419/141‑1 (MQ=255)                                                                                          gCACCTGGATGGAATGACTGTGGAAGAGGCTTTAATTACAATGAATAAATAAGATCAGCAAGTCCCGTTAGCAGTTCGAAAGGCAATACCACAACTGACAAAAGTAATTAAAAAAACAATTGCACAGTATAAAAAGGGTgg                                                <  1:485357/141‑1 (MQ=255)                                                                                            aTCTTGATGAAATGACTGTGGAAGAGGCTTTAATTACGATGAATAAATAAGATCAGCAAGTCCCGTTAGCAGTTCGAAAGGCAATACCACAATTGACAAAAGTAATTAAAAAAACAATTGCACACAATAAAAAGGGGGGAc                                              >  1:234844/1‑141 (MQ=255)                                                                                                              tgGAAGAGGCTTTAATTACGATGAATAAATAAGATCAGCAAGTCCCGTTAGCAGTTCGAAAGGCAATACCACAATTGACAAAAGTAATTAAAAAAACAATTGCACAGTATAAAAAGGGTGGACGATTGATTTATATCGGTg                            >  1:286655/1‑141 (MQ=255)                                                                                                                     ggCTTTAATTACGATGAATAAATAAGATCAGCAAGTCCCGTTAGCAGTTCGGAAGGCAATACCACTATTGACAAAAGTAATTAAAAAAACAATTGCACAGTATAAAAAGGGTCGACGATTGATTTATATCGGTGCAGGTAc                     >  1:345162/1‑141 (MQ=255)                                                                                                                        tttAATTACGATGAATAAATAAAAACCACTAGGCCCGTTAGCAGTTCGAAAGGCAATACCACAATTGACAAAAGTAATTAAAAAAACAATTGCACAGTATAAAAAGGGTGGACGATTGATTTATATCGGTGCAGGTACGAg                  >  1:365830/1‑141 (MQ=255)                                                                                                                               aCGATGAATAAATAAGATCAGCAAGTCCCGTTAGCAGTTCGAAAGGCAATACCACAATTGACAAAAGTAATTAAAAAAACAATTGCACAGTATAAAAAGGGTGGACGATTGATTTATATCGGTGCAGGTACAAGTGGAAgg           >  1:374072/1‑141 (MQ=255)                                                                                                                                      ataaataaGATCAGCAAGTCCCGTTAGCAGTTCGAAAGGCAATACCACAATTGACAAAAGTAATTAAAAAAACAATTGCACAGTATAAAAAGGGTGGACGATTTATTTATATCGGTGCAGGTACAAGTGGAAGGTTGGgtg    >  1:11885/1‑141 (MQ=255)                                                                                                                                      ataaataaGATCAGCAAGTCCCGTTAGCAGTTCGAAAGGCAATACCACAATTGACAAAAGTAATTAAAAAAACAATTGCACAGTATAAAAAGGGTGGACGATTGATTTATATCGGTGCAGGTACAAGTGGAGGGTTGGgtg    >  1:450399/1‑141 (MQ=255)                                                                                                                                        aaaTAAGATCAGCAAGTCCCGTTAGCAGTTCGGAAGGCAATACCACAATTGACAAAAGTAATTAAAAAAACAATTGCACAGTATAAAAAGGGTGGGCGATTGAGTTATATCGGTGCAGGTACAAGTGGAAGGTTGGgtgta  >  1:104929/1‑140 (MQ=255)                                                                                                                                           |                                                                                                                                          GCCGAATGATACATTTGATTTTCAAAAAACTAGGGAGTGTAAGAAGTGATGGAAAATAGTACGACCGAAGCGCGTAATGAAGCGACGATGCATCTTGATGAAATGACTGTGGAAGAGGCTTTAATTACGATGAATAAAGAAGATCAGCAAGTCCCGTTAGCAGTTCGAAAGGCAATACCACAATTGACAAAAGTAATTAAAAAAACAATTGCACAGTATAAAAAGGGTGGACGATTGATTTATATCGGTGCAGGTACAAGTGGAAGGTTGGGTGTC  >  CP000730/224866‑225141

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑