Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,601,428	A→G	G30G (GGT→GGC)	dnaT ←	DNA biosynthesis protein (primosomal protein I)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,601,428	0	A	G	100.0%	152.5 / NA	44	G30G (GGT→GGC)	dnaT	DNA biosynthesis protein (primosomal protein I)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/29);  total (15/29)

TGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACA  >  NC_000913/4601285‑4601575                                                                                                                                                |                                                                                                                                                    tGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACGCCGCCtt                                                                                                                                              >  2:28220/1‑151 (MQ=255)   cGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACGCCGCCTTc                                                                                                                                             >  2:595987/1‑150 (MQ=255)    ggTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACGCCGCCTTc                                                                                                                                             >  1:375917/1‑149 (MQ=255)                agTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACGCCGCCTTCAGCTTTTGCCAGa                                                                                                                                >  1:1057032/1‑150 (MQ=255)                   tGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACGCCGCCTTCAGCTTTTGCCAGAACg                                                                                                                             <  2:362829/150‑1 (MQ=255)                     aTCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACGCCGCCTTCAGCTTTTGCCAGAACGGtt                                                                                                                          <  2:699843/151‑1 (MQ=255)                       cGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACGCCGCCTTCAGCTTTTGCCAGAACGGTTTg                                                                                                                        <  2:589433/151‑1 (MQ=255)                                      cGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACGCCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCa                                                                                                          <  1:637094/150‑1 (MQ=255)                                      cGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACGCCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCa                                                                                                          <  2:381281/150‑1 (MQ=255)                                         ttCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACGCCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCAGGGc                                                                                                      >  1:265224/1‑151 (MQ=255)                                          tCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACGCCGCCTTCAGCTTTTGCCAGATCGGTTTGGTGATCGTGTTCCAGGGc                                                                                                      <  2:503756/150‑1 (MQ=255)                                          tCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACGCCGCCTTCAGCTTTTGCCAGAACGGATTGGTGATCGTGTACCAGGGc                                                                                                      >  1:809967/1‑150 (MQ=255)                                               gACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACGCCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCAGGGCGTCAAt                                                                                                <  2:54578/151‑1 (MQ=255)                                               gACGCGCCAGCTTTTCTTCCAGCGCCAGCACTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACGCCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCAGGGCGTCAAt                                                                                                >  1:381737/1‑151 (MQ=255)                                                   cgcCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACGCCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCAGGGCGTCAATACCa                                                                                            <  2:330947/151‑1 (MQ=255)                                                                tcCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACGCCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCggg                                                                                <  2:408625/150‑1 (MQ=255)                                                                tcCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACGCCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCgg                                                                                 <  2:513386/149‑1 (MQ=255)                                                                 ccAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACGCCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCgggg                                                                               <  2:146029/150‑1 (MQ=255)                                                                    gcgcCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCTCGCCGCCTTCAGCTTGTGCCAGAACGGTTTGGTGATCGTGTACCAGGGCGGCAATCCCAACGACGACCGGGGTCaa                                                                           <  2:158643/151‑1 (MQ=255)                                                                        cagcagTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACGCCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAAct                                                                       <  2:44378/151‑1 (MQ=255)                                                                               tCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACGCCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGa                                                                  >  2:657271/1‑149 (MQ=255)                                                                                 aGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACGCCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGAc                                                              >  1:881617/1‑151 (MQ=255)                                                                                         gtgcaggcgtgACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACGCCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGtt                                                      <  2:273379/151‑1 (MQ=255)                                                                                            caggcgtgACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACGCCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGTTtc                                                    >  2:993152/1‑150 (MQ=255)                                                                                              ggcgtgACGGCGTAAAACGCCGGGGCATTGTTAGCAAATACGGCAGCCACGCCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCAGGTCGGCAATACCAACGACGTCCGGGGTGAAAACTCTGGAAGACATACTGTTtctc                                                  <  2:540326/150‑1 (MQ=255)                                                                                                 gtgGCGGCATAAAACGCCGGGGGATTGTTAACAAAGACGGCAACCACGGCGCCTTCAGCTTTTGCCAGAACGTTTTGGTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGTTTCTCaa                                                <  2:520244/149‑1 (MQ=255)                                                                                                   gACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACGCCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGTTTCTCAAata                                             <  1:294265/150‑1 (MQ=255)                                                                                                       gCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACGCCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGTTTCTCAAATATGGa                                         >  1:872481/1‑150 (MQ=255)                                                                                                        cATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACGCCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGTTTCTCAAATATGGAAc                                       <  2:667029/151‑1 (MQ=255)                                                                                                                gCCGGGGCATTGTTAGCAAATACGGCAACCACGCCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGTTTCTCAAATATGGAACGATGAACg                               <  2:275059/151‑1 (MQ=255)                                                                                                                 ccGGGGCATTGTTAGCAAATACGGCAACCACGCCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGTTTCTCAAATATGGAACGATGAACga                              <  1:790329/151‑1 (MQ=255)                                                                                                                  cGGGGCATTGTTAGCAAATACGGCAACCACGCCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGTTTCTCAAATATGGAACGATGAACgag                             <  2:737462/151‑1 (MQ=255)                                                                                                                  cGGGGCATTGTTAGCAAATACGGCAACCACGCCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGTTTCTCAAATATGGAACGATGAACga                              <  2:729658/150‑1 (MQ=255)                                                                                                                    gggCATTGTTAGCAAATACGGCAACCACGCCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAGCTCTGGAAGACAGACTGTTTCTCAAATATGGAACGATGAACgaga                            <  2:218186/150‑1 (MQ=255)                                                                                                                       cATTGTTAGCAAATACGGCAACCACGCCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGTTTCTCAAATATGGAACGATGAACGAGAAggg                        >  1:743330/1‑151 (MQ=255)                                                                                                                            ttAGCAAATACGGCAACCACGCCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCAGGGCGTCAATACCAACGACGTCTGGGGTCAAAACTCTGGAAGACATACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGAt                    >  1:1022593/1‑150 (MQ=255)                                                                                                                              aGCAAATACGGCAACCACGCCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGTGATAAc                 <  1:411989/151‑1 (MQ=255)                                                                                                                                 aaaTACGGCAACCACGCCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACttt              <  2:188185/151‑1 (MQ=255)                                                                                                                                     aCGGCAACCACGCCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTtata          <  1:868606/151‑1 (MQ=255)                                                                                                                                     aCGGCAACCACGCCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTtata          <  1:892213/151‑1 (MQ=255)                                                                                                                                     aCGGCAACCACGCCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTtat           >  1:873587/1‑150 (MQ=255)                                                                                                                                       ggCAACCACGCCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGt        <  1:56101/151‑1 (MQ=255)                                                                                                                                       ggCAACCACGCCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGt        <  2:786121/151‑1 (MQ=255)                                                                                                                                         cAACCACGCCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTa       >  2:65024/1‑150 (MQ=255)                                                                                                                                             cacgCCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAacaca  <  2:207998/151‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                    TGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCAGGGCGTCAATACCAACGACGTCCGGGGTCAAAACTCTGGAAGACATACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACA  >  NC_000913/4601285‑4601575

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