Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,244,074	G→A	T223I (ACC→ATC)	uxaC ←	uronate isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,244,074	0	G	A	100.0%	131.8 / NA	42	T223I (ACC→ATC)	uxaC	uronate isomerase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/19);  total (23/19)

CTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGGTCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCGCCATGTAGTCGTTAAAGGTCGCCTGTTCGATGTTGAAGGCTTTGTCCGGACGCCAGCTCGGCAGCACTTTGA  >  NC_000913/3243939‑3244209                                                                                                                                        |                                                                                                                                        ctcAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCa                                                                                                                                      >  2:645015/1‑139 (MQ=255)    aGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAggg                                                                                                                                   <  1:376150/139‑1 (MQ=255)                                  aTGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGt                                                                                                     <  2:729483/139‑1 (MQ=255)                                           aGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCgg                                                                                            <  1:433818/139‑1 (MQ=255)                                             ttCCACTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGaa                                                                                          >  2:812901/1‑139 (MQ=255)                                                  cTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTc                                                                                     >  1:922726/1‑139 (MQ=255)                                                  cTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTc                                                                                     >  2:1093773/1‑139 (MQ=255)                                                  cTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTc                                                                                     >  1:627368/1‑139 (MQ=255)                                                           cTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCtt                                                                            <  2:819162/139‑1 (MQ=255)                                                             tCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCg                                                                          <  1:335564/139‑1 (MQ=255)                                                                        caccacATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAg                                                                               <  1:1323646/123‑1 (MQ=255)                                                                        caccacATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAg                                                                               >  2:1323646/1‑123 (MQ=255)                                                                           cacaTCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAggg                                                                                                                                   >  2:824334/1‑68 (MQ=255)                                                                           cacaTCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAggg                                                                                                                                   <  1:824334/68‑1 (MQ=255)                                                                                ccAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCa                                                                                                                          >  1:546845/1‑72 (MQ=255)                                                                                ccAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCa                                                                                                                          <  2:546845/72‑1 (MQ=255)                                                                                   gcgcGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCGCCATGTAGTCGTTAAAGGTCGc                                                    <  2:922726/139‑1 (MQ=255)                                                                                     gcgTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACt                                                                                       <  2:96094/102‑1 (MQ=255)                                                                                     gcgTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACt                                                                                       >  1:96094/1‑102 (MQ=255)                                                                                     gcgTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCGCCATGTAGTCGTTAAAGGTCGCCt                                                  >  2:232200/1‑139 (MQ=255)                                                                                      cgTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGc                                                                                   <  1:325842/105‑1 (MQ=255)                                                                                      cgTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGc                                                                                   >  2:325842/1‑105 (MQ=255)                                                                                      cgTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCGCCATGTAGTCGTTAAAGGTCGCCTg                                                 >  2:1469404/1‑139 (MQ=255)                                                                                        tGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCGCCAt                                                                      <  2:863675/116‑1 (MQ=255)                                                                                        tGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCGCCAt                                                                      >  1:863675/1‑116 (MQ=255)                                                                                         ggTCAGACACTTTACAGCCGTTAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCGCCATGTAGTCGTTAAAGGTCGCCTGTTc                                              >  1:1090910/1‑139 (MQ=255)                                                                                         ggTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCg                                                                          <  1:645696/111‑1 (MQ=255)                                                                                         ggTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCg                                                                          >  2:645696/1‑111 (MQ=255)                                                                                         ggTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCGCCATGTAGTCGTTAAAGGTCGCCTGTTc                                              >  2:553611/1‑139 (MQ=255)                                                                                            cAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCGCCATGTAg                                                                  <  2:290515/116‑1 (MQ=255)                                                                                            cAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCGCCATGTAg                                                                  >  1:290515/1‑116 (MQ=255)                                                                                              gACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCGCCATGTAGTCGTTAAAGGTCGCCTGTTCGATGt                                         >  2:1244352/1‑139 (MQ=255)                                                                                                   tttACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCGCCATGTAGTCGTTAAAGGTCGCCTGTTCGATGTTGAAg                                    >  2:1266122/1‑139 (MQ=255)                                                                                                     tACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGcc                                                                                  >  1:122615/1‑91 (MQ=255)                                                                                                     tACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGcc                                                                                  <  2:122615/91‑1 (MQ=255)                                                                                                                           aTCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCg                                                                          >  1:668326/1‑77 (MQ=255)                                                                                                                           aTCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCg                                                                          <  2:668326/77‑1 (MQ=255)                                                                                                                             cccgACGTTTGCTCAGGGCCGTTTGCAGGTCAGCAAAGAGGTGAATGACGGTATCGGAAACTTCGCCCAGCTTCGCCATGTAGTCGTTAAAGGTCGCCTGTTCGATGTTGAAGGCTTTGTCCGGACGCCAGCTCGgcag          <  1:812901/136‑1 (MQ=255)                                                                                                                             ccAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCGCCATGTAGTCGTTAAAGGTCGCCTGTTCGATGTTGAAGGCTTTGTCCGGACGCCAGCTCGgcag          <  2:1090910/139‑1 (MQ=255)                                                                                                                               aGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCGCCATGTAGTCGTTAAAGGTCGCCTGTTCGATGTTGAAGGCTTTGTCCGGACGCCAGCTCGgcagca        <  1:553611/139‑1 (MQ=255)                                                                                                                                 aCGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAACTTCGCCATGTAGTCGTTAAAGGTCGCCTGTTCGATGTTGAAGGCTTTGTCCGGACGCCAGCTCGGCAGCACt      >  1:876666/1‑139 (MQ=255)                                                                                                                                   gTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCGCCATGTAGTCGTTAAAGGTCGCCTGTTCGATGTTGAAGGCTTTGTCCGGACGCCAGCTCGGCAGCACttt    <  1:1244352/139‑1 (MQ=255)                                                                                                                                     ttGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCGCCATGTAGTCGTTAAAGGTCGCCTGTTCGATGTTGAAGGCTTTGTCCGGACGCCAGCTCGGCAGCACTTTGa  >  2:1285593/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                        CTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGGTCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCGCCATGTAGTCGTTAAAGGTCGCCTGTTCGATGTTGAAGGCTTTGTCCGGACGCCAGCTCGGCAGCACTTTGA  >  NC_000913/3243939‑3244209

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