Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	105,107	T→C	Y376H (TAT→CAT)	ftsA →	ATP‑binding cell division FtsK recruitment protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	105,107	0	T	C	100.0%	105.3 / NA	28	Y376H (TAT→CAT)	ftsA	ATP‑binding cell division FtsK recruitment protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/16);  total (12/16)

TGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACGTGTTACAGCATCAGTT  >  NC_000913/104977‑105196                                                                                                                                   |                                                                                          tGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATCATTCGAcg                                                                                   >  2:83517/1‑139 (MQ=255)        aCCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATCATTCGACGGCGGTgg                                                                            <  2:165745/139‑1 (MQ=255)        aCCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATCATTCGACGGCGGTgg                                                                            >  1:165745/1‑139 (MQ=255)        aCCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATCATTCGACGGCGGTgg                                                                            <  2:15461/139‑1 (MQ=255)          cGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATACTCAGGAGCCGTATCATTCGACGGCGGTGGGa                                                                          <  2:190268/139‑1 (MQ=255)           ggTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATCATTCGACGGCGGTGGGAt                                                                         >  2:472113/1‑139 (MQ=255)                cgcagcgcagATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATCATTCGACGGCGGTGGGATTGCtt                                                                    <  1:173834/139‑1 (MQ=255)                cgcagcgcagATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATCATTCGACGGCGGTGGGATTGCtt                                                                    <  1:83517/139‑1 (MQ=255)                          aTCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATCATTCGACGGCGGTGGGATTGCTTCACTATGGGa                                                          <  2:280133/139‑1 (MQ=255)                           tCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATCATTCGACGGCGGTGGGATTGCTTCACTATGGGaa                                                         <  2:1904/139‑1 (MQ=255)                                                     cgcgTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATCATTCGa                                                                                     >  1:366191/1‑85 (MQ=255)                                                     cgcgTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATCATTCGa                                                                                     <  2:366191/85‑1 (MQ=255)                                                     cgcgTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATCATTCGAcggcgg                                                                               <  1:273946/91‑1 (MQ=255)                                                     cgcgTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATCATTCGAcggcgg                                                                               >  2:273946/1‑91 (MQ=255)                                                     cgcgTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATCATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAg                                  >  2:285257/1‑136 (MQ=255)                                                     cgcgTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATCATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAg                                  <  1:285257/136‑1 (MQ=255)                                                        gtgtTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATCATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAg                            >  1:367295/1‑139 (MQ=255)                                                             tCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATCATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGaa                       <  1:472113/139‑1 (MQ=255)                                                                     aaGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATCATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACgtgt               <  2:367295/139‑1 (MQ=255)                                                                       gTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATCATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACGTGTTa             <  1:443879/139‑1 (MQ=255)                                                                         gCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATCATTCGACGGCGGt                                                                              <  2:346381/72‑1 (MQ=255)                                                                         gCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATCATTCGACGGCGGt                                                                              >  1:346381/1‑72 (MQ=255)                                                                                  cgcgcCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATCATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACGTGTTACAGCATCAGtt  >  1:455765/1‑139 (MQ=255)                                                                                  cgcgcCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATCATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACGTGTTACAGCATCAGtt  >  2:497407/1‑139 (MQ=255)                                                                                                 tACCGGTTTAACGGATTATGCTCAGGAGCCGTATCATTCGACGGCGGTgg                                                                            <  2:252736/50‑1 (MQ=255)                                                                                                 tACCGGTTTAACGGATTATGCTCAGGAGCCGTATCATTCGACGGCGGTgg                                                                            >  1:252736/1‑50 (MQ=255)                                                                                                 tACCGGTTTAACGGATTATGCTCAGGAGCCGTATCATTCGACGGCGGTGGGATTGCTTCa                                                                  <  1:449257/60‑1 (MQ=255)                                                                                                 tACCGGTTTAACGGATTATGCTCAGGAGCCGTATCATTCGACGGCGGTGGGATTGCTTCa                                                                  >  2:449257/1‑60 (MQ=255)                                                                                                                                   |                                                                                          TGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACGTGTTACAGCATCAGTT  >  NC_000913/104977‑105196

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