Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,395,394	T→C	T289A (ACG→GCG)	yhdP ←	DUF3971‑AsmA2 domains protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,395,394	0	T	C	100.0%	130.7 / NA	38	T289A (ACG→GCG)	yhdP	DUF3971‑AsmA2 domains protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/18);  total (20/18)

GTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGTAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGTTTGCTTCTCGCCCAACCAGCTGGCACCGCCCTGTTTCAGCCAGACGTCACCGCCGGTTACATCGCCTTTGTCGATCGTCATC  >  NC_000913/3395269‑3395515                                                                                                                              |                                                                                                                          gTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGTCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGCAATATGCGCGGTc                                                                                                              <  1:1485986/139‑1 (MQ=255)       gccgccAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGCAATATGCGCGGTCAGATTa                                                                                                        <  1:586338/139‑1 (MQ=255)       gccgccAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGCAATATGCGCGGTCAGATTa                                                                                                        <  1:442794/139‑1 (MQ=255)               cGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGCAATATGCGCGGTCAGATTATCCACCGa                                                                                                <  2:907100/139‑1 (MQ=255)                           gTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGCAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGt                                                                                    >  2:1265076/1‑139 (MQ=255)                                 aGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGCAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGTTTGCtt                                                                              <  1:1265076/139‑1 (MQ=255)                                   gCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGCAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGTTTGCTtct                                                                            >  1:1148981/1‑139 (MQ=255)                                   gCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGCAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGTTTGCTtct                                                                            >  2:732438/1‑139 (MQ=255)                                      aGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGCAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGTTTGCTTCTCGc                                                                         >  2:64039/1‑139 (MQ=255)                                        cGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGCAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGTTTGCTTCTCGccc                                                                       >  1:1226149/1‑139 (MQ=255)                                        cGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGCAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGTTTGCTTCTCGccc                                                                       >  2:1512985/1‑139 (MQ=255)                                        cGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGCAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGTTTGCTTCTCGccc                                                                       <  1:942255/139‑1 (MQ=255)                                         gTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGCAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGTTTGCTTCTCGCCCa                                                                      <  2:535854/139‑1 (MQ=255)                                               gCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGCAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGTTTGCTTCTCGCCCAACCAGc                                                                <  1:1291477/139‑1 (MQ=255)                                                  ccGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGCAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGTTTGCTTCTCGCCCAACCAGCTgg                                                             >  1:332602/1‑139 (MQ=255)                                                     cccGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGCAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGTGTGCTTCTCGCCCAACCAGCTGGCAc                                                          <  1:732438/137‑1 (MQ=255)                                                     cTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGCAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGTTTGCTTCTCGCCCAACCAGCTGGCAc                                                          >  1:441636/1‑139 (MQ=255)                                                                    ccGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGCAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGTTTGCt                                                                               <  2:21926/103‑1 (MQ=255)                                                                    ccGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGCAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGTTTGCt                                                                               >  1:21926/1‑103 (MQ=255)                                                                                 tCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGCAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGTTTGCTTCTCGCCCAACCAGCTGGCACCGCCCTGTTTCAGCCAGACGTCAccgcc                              >  1:1428014/1‑139 (MQ=255)                                                                                     tgtATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGCAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGTTTGCTTCTCGCCCAACCAGCTGGCACCGCCCTGTTTCAGCCAGACGTCACCGCCGGtt                          >  2:1103762/1‑139 (MQ=255)                                                                                            ggAATAGAGAACTGCCAGCCCGGATTTTCACGCGCAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGTTTGCTTCTCGCCCaa                                                                     >  1:607470/1‑89 (MQ=255)                                                                                            ggAATAGAGAACTGCCAGCCCGGATTTTCACGCGCAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGTTTGCTTCTCGCCCaa                                                                     <  2:607470/89‑1 (MQ=255)                                                                                                   agaACTGCCAGCCCGGATTTTCACGCGCAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGTTTGCTTCTCGCCCAACCAGCTGGCACCGCCCTGTTTCAGCCa                                          <  1:496081/109‑1 (MQ=255)                                                                                                   agaACTGCCAGCCCGGATTTTCACGCGCAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGTTTGCTTCTCGCCCAACCAGCTGGCACCGCCCTGTTTCAGCCa                                          >  2:496081/1‑109 (MQ=255)                                                                                                   agaACTGCCAGCCCGGATTTTCACGCGCAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGTTTGCTTCTCGCCCAACCAGCTGGCACCGCCCTGTTTCAGCCAGACGTCACCGCCGGTTACATCGCCTTTGTc            >  2:299973/1‑139 (MQ=255)                                                                                                      aCTGCCAGCCCGGATTTTCACGCGCAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGTTTGCTTCTCGCCCAACCAGCTGGCACCGCCCTGTTTCAGCCAGACGTCACCGCCGGTTACATCGCCTTTGTCGAt         <  2:1226149/139‑1 (MQ=255)                                                                                                        tGCCAGCCCGGATTTTCACGCGCAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGTTTGCTTCTCGCCCAACCAGCTGGCAc                                                          <  2:1254948/88‑1 (MQ=255)                                                                                                        tGCCAGCCCGGATTTTCACGCGCAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGTTTGCTTCTCGCCCAACCAGCTGGCAc                                                          >  1:1254948/1‑88 (MQ=255)                                                                                                        tGCCAGCCCGGATTTTCACGCGCAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGTTTGCTTCTCGCCCAACCAGCTGGCACCGCCCTGTTTCAGCCAGACGTCACCGCCGGTTACATCGCCTTTGTCGATCg       >  1:1554004/1‑139 (MQ=255)                                                                                                        tGCCAGCCCGGATTTTCACGCGCAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGTTTGCTTCTCGCCCAACCAGCTGGCACCGCCCTGTTTCAGCCAGACGTCACCGCCGGTTACATCGCCTTTGTCGATCg       <  1:64039/139‑1 (MQ=255)                                                                                                            aGCCCGGATTTTCACGCGCAATATGCGCGGTCAGATTATc                                                                                                      >  2:157119/1‑40 (MQ=255)                                                                                                            aGCCCGGATTTTCACGCGCAATATGCGCGGTCAGATTATc                                                                                                      <  1:157119/40‑1 (MQ=255)                                                                                                             gCCCGGATTTTCACGCGCAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGTTTGCTTCTCGCCCAACCAGCTGGCACCGCCCTGTTTCAGCCAGACGTCACCGCCGGTTACATCGCCTTTGTCGATCGTCATc  >  1:1296558/1‑139 (MQ=255)                                                                                                               ccGGATTTTCACGCGCAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGTTTGCTTCTCGCCCAACCAGCTGGCACCGCCCTGTTTCAGCCAGACGTCAc                                  <  2:42779/105‑1 (MQ=255)                                                                                                               ccGGATTTTCACGCGCAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGTTTGCTTCTCGCCCAACCAGCTGGCACCGCCCTGTTTCAGCCAGACGTCAc                                  >  1:42779/1‑105 (MQ=255)                                                                                                                  gATTTTCACGCGCAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGTTTGCTTCTCGCCCAACCAGCTGGCACCGCCCTGTTTCAGCCAGACGTCAccgccg                             <  2:1610696/107‑1 (MQ=255)                                                                                                                  gATTTTCACGCGCAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGTTTGCTTCTCGCCCAACCAGCTGGCACCGCCCTGTTTCAGCCAGACGTCAccgccg                             >  1:1610696/1‑107 (MQ=255)                                                                                                                              |                                                                                                                          GTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGTAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGTTTGCTTCTCGCCCAACCAGCTGGCACCGCCCTGTTTCAGCCAGACGTCACCGCCGGTTACATCGCCTTTGTCGATCGTCATC  >  NC_000913/3395269‑3395515

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