Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,425,834	G→A	*105* (TAG→TAA)	priB →	primosomal protein N

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,425,834	0	G	A	100.0%	180.6 / NA	56	*105* (TAG→TAA)	priB	primosomal protein N
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/26);  total (30/26)

GGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGCCGT  >  NC_000913/4425699‑4425967                                                                                                                                        |                                                                                                                                      ggCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCa                                                                                                                                    >  2:337684/1‑139 (MQ=255)                 ataACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTc                                                                                                                    >  2:146424/1‑139 (MQ=255)                   aaCGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCcg                                                                                                                  <  1:337684/139‑1 (MQ=255)                     cGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCcgtc                                                                                                                >  1:822796/1‑139 (MQ=255)                        tCGGCAGTTGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGACAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATCTGGCACGTTATTTCCgccgtc                                                                                                             >  1:1031532/1‑139 (MQ=255)                          ggCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGc                                                                                                           >  2:1752551/1‑139 (MQ=255)                          ggCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGc                                                                                                           >  1:1554330/1‑139 (MQ=255)                                cGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTc                                                                                                     <  1:1877026/139‑1 (MQ=255)                                          ttCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGt                                                                                                                          <  2:1678569/108‑1 (MQ=255)                                          ttCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGt                                                                                                                          >  1:1678569/1‑108 (MQ=255)                                          ttCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCa                                                                                           >  1:847652/1‑139 (MQ=255)                                          ttCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCa                                                                                           >  2:31916/1‑139 (MQ=255)                                          ttCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCa                                                                                           >  1:1674413/1‑139 (MQ=255)                                           tCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCAc                                                                                          <  1:489367/139‑1 (MQ=255)                                                 gTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGa                                                                                    >  2:1220428/1‑139 (MQ=255)                                                 gTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGa                                                                                    >  1:1469324/1‑139 (MQ=255)                                                         catGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTc                                                                            <  2:1469324/139‑1 (MQ=255)                                                         catGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTc                                                                            <  1:1220428/139‑1 (MQ=255)                                                              cacaAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAagag                                                                       <  2:1674413/139‑1 (MQ=255)                                                                caAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAt                                                                     <  2:822796/139‑1 (MQ=255)                                                                 aaGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATc                                                                    >  2:253812/1‑139 (MQ=255)                                                                 aaGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATc                                                                    <  1:253812/139‑1 (MQ=255)                                                                    gCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGAc                                                                 <  2:847652/139‑1 (MQ=255)                                                                        aGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACtata                                                             >  1:1647575/1‑139 (MQ=255)                                                                                cTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGa                                                                  >  2:1746790/1‑126 (MQ=255)                                                                                cTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGa                                                                  <  1:1746790/126‑1 (MQ=255)                                                                                  gAGCAAAATGGTTTTGCATGCCGAGTAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGc                                                   >  1:978963/1‑139 (MQ=255)                                                                                  gAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGc                                                   >  2:1027812/1‑139 (MQ=255)                                                                                   aGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAt                                                                     >  1:1799460/1‑120 (MQ=255)                                                                                   aGCAAAATGGTTGTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAt                                                                     <  2:1799460/120‑1 (MQ=255)                                                                                    gCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCg                                                                               >  1:1555528/1‑109 (MQ=255)                                                                                    gCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCg                                                                               <  2:1555528/109‑1 (MQ=255)                                                                                       aaaTGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGc                                              <  1:146424/139‑1 (MQ=255)                                                                                                      ccGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATc                                                                    <  1:554596/102‑1 (MQ=255)                                                                                                      ccGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATc                                                                    >  2:554596/1‑102 (MQ=255)                                                                                                           cAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAt                                                                     >  2:942667/1‑96 (MQ=255)                                                                                                           cAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAt                                                                     <  1:942667/96‑1 (MQ=255)                                                                                                           cAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGaa                          <  2:1756258/139‑1 (MQ=255)                                                                                                            aGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAt                                                                     <  2:1446206/95‑1 (MQ=255)                                                                                                            aGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAt                                                                     >  1:1446206/1‑95 (MQ=255)                                                                                                                tgaattgaTAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCgg                     <  2:1647575/139‑1 (MQ=255)                                                                                                                      gaTAGATTCTTGAGACTAACCATATGGTACGTTATTTCCGGGGTCGCAAGTTCTGCCGTTTCACCGCGGAATGCGTTCAAGAGATCCACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGAt               <  1:1142758/139‑1 (MQ=255)                                                                                                                        tAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTc                                                                            >  1:736873/1‑76 (MQ=255)                                                                                                                        tAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTc                                                                            <  2:736873/76‑1 (MQ=255)                                                                                                                        tAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTg             >  1:477302/1‑139 (MQ=255)                                                                                                                           aTTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTcc          <  1:1248716/139‑1 (MQ=255)                                                                                                                            ttCTGGAGACTCACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGa                                                                                    <  2:1269508/64‑1 (MQ=255)                                                                                                                            ttCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGa                                                                                    >  1:1269508/1‑64 (MQ=255)                                                                                                                            ttCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAt                                                                     >  1:730560/1‑79 (MQ=255)                                                                                                                            ttCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAt                                                                     <  2:730560/79‑1 (MQ=255)                                                                                                                             tCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCa        >  1:1700243/1‑139 (MQ=255)                                                                                                                                 gagaCTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGcc    >  1:509479/1‑139 (MQ=255)                                                                                                                                 gagaCTAACCATATGGCACATTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAagag                                                                       <  2:378821/72‑1 (MQ=255)                                                                                                                                 gagaCTAACCATATGGCACATTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAagag                                                                       >  1:378821/1‑72 (MQ=255)                                                                                                                                   gaCTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGCCGt  <  1:1027812/139‑1 (MQ=255)                                                                                                                                     cTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATc                                                                    >  2:1126468/1‑71 (MQ=255)                                                                                                                                        |                                                                                                                                      GGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGCCGT  >  NC_000913/4425699‑4425967

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