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%	133.7 / NA	42	T223I (ACC→ATC)	uxaC	uronate isomerase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/21);  total (21/21)

CGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGGTCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCGCCATGTAGTCGTTAAAGGTCGCCTGTTCGATGTTGAAGGCTTTGTCCGGACGCCAGCT  >  NC_000913/3243949‑3244195                                                                                                                              |                                                                                                                          cGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTg                                                                                                              <  2:791561/139‑1 (MQ=255)      gccagACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGt                                                                                                         >  2:672261/1‑139 (MQ=255)      gccagACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGt                                                                                                         >  2:234191/1‑139 (MQ=255)      gccagACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGt                                                                                                         >  2:1004/1‑139 (MQ=255)           aCGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCa                                                                                                    <  1:362345/139‑1 (MQ=255)                gcgcCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAgcg                                                                                               <  1:1004/139‑1 (MQ=255)                 cgcCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAgcgg                                                                                              >  2:389193/1‑139 (MQ=255)                              tCCAGTTCTGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTAt                                                                                 >  1:721434/1‑139 (MQ=255)                                        cTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTc                                                                       >  2:1031055/1‑139 (MQ=255)                                        cTTCGTTCGCTTCAGCAAACATCACCACATCCAGCCCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTc                                                                       >  2:398974/1‑139 (MQ=255)                                             ttcgCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCa                                                                  <  1:389193/139‑1 (MQ=255)                                                     aGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCGcc                                                          <  2:136152/139‑1 (MQ=255)                                                           cATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCGCCATGTAg                                                    <  1:1031055/139‑1 (MQ=255)                                                            aTCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCGCCATGTAGt                                                   >  1:1062852/1‑139 (MQ=255)                                                                 cacaTCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCa                                                                                                                   >  2:482532/1‑70 (MQ=255)                                                                 cacaTCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCa                                                                                                                   <  1:482532/70‑1 (MQ=255)                                                                 cacaTCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCGCCATGTAGTCGTTa                                              >  1:267330/1‑139 (MQ=255)                                                                   caTCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTg                                                                                                              <  1:201373/73‑1 (MQ=255)                                                                   caTCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTg                                                                                                              >  2:201373/1‑73 (MQ=255)                                                                    aTCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGCGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCGCCATGTAGTCGTTAAAg                                           <  2:721434/139‑1 (MQ=255)                                                                     tCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCGCCATGTAGTCGTTAAAgg                                          >  1:601382/1‑139 (MQ=255)                                                                      ccAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAg                                                                                                 <  1:840165/83‑1 (MQ=255)                                                                      ccAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAg                                                                                                 >  2:840165/1‑83 (MQ=255)                                                                        aGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCGCCATGTAGTCGTTAAAGGTCg                                       <  1:396000/139‑1 (MQ=255)                                                                         gcgcGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCGCCATGTAGTCGTTAAAGGTCGc                                      <  1:234191/139‑1 (MQ=255)                                                                         gcgcGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCGCCATGTAGTCGTTAAAGGTCGc                                      <  1:672261/139‑1 (MQ=255)                                                                              tGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCGCCATGTAGTCGTTAAAGGTCGCCTGtt                                 >  1:288205/1‑139 (MQ=255)                                                                                gTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCGCCATGTAGTCGTTAAAGGTCGCCTGTTCg                               <  1:528872/139‑1 (MQ=255)                                                                                gTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCGCCATGTAGTCGTTAAAGGTCGCCTGTTCg                               <  2:1062852/139‑1 (MQ=255)                                                                                gTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCGCCATGTAGTCGTTAAAGGTCGCCTGTTCg                               >  1:463362/1‑139 (MQ=255)                                                                                    gACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCGCCATGTAGTCGTTAAAGGTCGCCTGTTCGATGt                           <  2:288205/139‑1 (MQ=255)                                                                                           tACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCACCCAGCTTCGCCATGTAGTCGTTAAAGGTCGCCTGTTCGATGTTGAAGGc                    >  2:164826/1‑139 (MQ=255)                                                                                                   tGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAg                                                                                                 <  1:220650/54‑1 (MQ=255)                                                                                                   tGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAg                                                                                                 >  2:220650/1‑54 (MQ=255)                                                                                                   tGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAg                                                                                                 <  1:876923/54‑1 (MQ=255)                                                                                                   tGAGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAg                                                                                                 >  2:876923/1‑54 (MQ=255)                                                                                                     aGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTg                                                                                                              >  1:840813/1‑39 (MQ=255)                                                                                                     aGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTg                                                                                                              <  2:840813/39‑1 (MQ=255)                                                                                                     aGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCtt                                                              <  2:1073662/87‑1 (MQ=255)                                                                                                     aGCGGCGAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCtt                                                              >  1:1073662/1‑87 (MQ=255)                                                                                                          cgAAGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCGCCATGTAGTCGTTAAAGGTCGCCTGTTCGATGTTGAAGGCTTTGTCCGGACGCCa     >  2:68819/1‑139 (MQ=255)                                                                                                             aGTGATCCAGACGTTTGATCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACGTCGCACAGCTTCGCCATGTAGTCGTTAAAGGTCGCCTGTTCGATGTTGAAGGCTTTGTCCGGACGCCAggt  <  2:1043352/139‑3 (MQ=255)                                                                                                                              |                                                                                                                          CGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGGTCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTCGCCATGTAGTCGTTAAAGGTCGCCTGTTCGATGTTGAAGGCTTTGTCCGGACGCCAGCT  >  NC_000913/3243949‑3244195

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: ‑