Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,094,234	T→C	T283A (ACA→GCA)	yggR ←	putative PilT family AAA+ ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,094,234	0	T	C	100.0%	136.7 / NA	39	T283A (ACA→GCA)	yggR	putative PilT family AAA+ ATPase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/18);  total (21/18)

GTCACAAACGCCCTTCCCCCACCCGGTGCTGATAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGTGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCACTGCCCGTAAACTACCTGCCAGTTGATTACGCACGGGGTCTTTTTCCTGCGCCGGAAATGAATC  >  NC_000913/3094099‑3094366                                                                                                                                        |                                                                                                                                     gTCACAAACGCCCTTCCCCCACCCGGTGCTGATAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGtt                                                                                                                                   >  2:855864/1‑139 (MQ=255)                      cccGGTGCTGATAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAgcgc                                                                                                              >  1:998583/1‑139 (MQ=255)                        cGGTGCTGATAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCa                                                                                                            <  1:855864/139‑1 (MQ=255)                          gTGCTGATAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCAcg                                                                                                          <  1:533865/139‑1 (MQ=255)                            gCTGATAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCAcgcg                                                                                                        <  1:836691/139‑1 (MQ=255)                                aTAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGtcct                                                                                                    <  2:998583/139‑1 (MQ=255)                                     tctGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTg                                                                                               >  2:899804/1‑139 (MQ=255)                                       tgctgAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGc                                                                                                                            >  2:918513/1‑108 (MQ=255)                                       tgctgAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGc                                                                                                                            <  1:918513/108‑1 (MQ=255)                                            aaaCGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATcc                                                                                        >  2:442310/1‑139 (MQ=255)                                                 ttATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCATCCTTCCTGTTTATCCACTTc                                                                                   >  1:67481/1‑139 (MQ=255)                                                   atcatcCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCAACTTCCa                                                                                 >  1:598363/1‑139 (MQ=255)                                                      atcCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCAc                                                                                                           >  2:891155/1‑110 (MQ=255)                                                      atcCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCAc                                                                                                           <  1:891155/110‑1 (MQ=255)                                                         cccACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCtttt                                                                           >  1:411652/1‑139 (MQ=255)                                                         cccACCTGCTGCCCGGTTTGAAT‑ACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTtg                                                                          >  2:860631/1‑139 (MQ=255)                                                            aCCTGCTGCCCGGTTTGAATAACATGCGGTAACTAGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAACCAGCCATTCAAATAGCGCCACGCGTCCTTCCTGTTTAACCAATTCCAGCTTTtgtg                                                                        >  2:26640/1‑139 (MQ=255)                                                                 ctgcCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCa                                                                                                                              >  2:178820/1‑80 (MQ=255)                                                                 ctgcCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCa                                                                                                                              >  2:147640/1‑80 (MQ=255)                                                                 ctgcCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCa                                                                                                                              <  1:147640/80‑1 (MQ=255)                                                                 ctgcCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCa                                                                                                                              <  1:178820/80‑1 (MQ=255)                                                                 ctgcCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGc                                                                   >  2:728968/1‑139 (MQ=255)                                                                   gcCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCAc                                                                 <  2:700945/139‑1 (MQ=255)                                                                   gcCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCAc                                                                 <  1:280313/139‑1 (MQ=255)                                                                             aaTAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCACTGCCCGTaaa                                                       <  2:910802/139‑1 (MQ=255)                                                                                       ggTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCATCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCACTGCCCGTAAACTACCTGCCa                                             <  2:67481/139‑1 (MQ=255)                                                                                         tAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAAt                                                                                                                         >  2:817918/1‑61 (MQ=255)                                                                                         tAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAAt                                                                                                                         <  1:817918/61‑1 (MQ=255)                                                                                                  gggtttCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCa                                                                                                                           <  2:620958/47‑1 (MQ=255)                                                                                                  gggtttCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCa                                                                                                                           >  1:620958/4‑50 (MQ=255)                                                                                                       tCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCAcg                                                                                                          <  2:517601/62‑1 (MQ=255)                                                                                                       tCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCAcg                                                                                                          >  1:517601/1‑62 (MQ=255)                                                                                                         ccTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGtccttcc                                                                                                 >  1:252445/1‑69 (MQ=255)                                                                                                         ccTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGtccttcc                                                                                                 <  2:252445/69‑1 (MQ=255)                                                                                                                 aaTCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCACTGc                                                              <  1:572103/96‑1 (MQ=255)                                                                                                                 aaTCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCACTGc                                                              >  2:572103/1‑96 (MQ=255)                                                                                                                              aCCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATc                                                                                         <  2:299434/56‑1 (MQ=255)                                                                                                                              aCCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATc                                                                                         >  1:299434/1‑56 (MQ=255)                                                                                                                                  cgGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCACTGCCCGTAAACTACCTGCCAGTTGATTACGCACGGGGTCTTTTTCCTGCGCCGGAAATGAATc  >  1:108966/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                     GTCACAAACGCCCTTCCCCCACCCGGTGCTGATAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGTGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCACTGCCCGTAAACTACCTGCCAGTTGATTACGCACGGGGTCTTTTTCCTGCGCCGGAAATGAATC  >  NC_000913/3094099‑3094366

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