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%	175.7 / NA	55	*105* (TAG→TAA)	priB	primosomal protein N
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/29);  total (26/29)

GGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGC  >  NC_000913/4425699‑4425964                                                                                                                                        |                                                                                                                                   ggCCATTACTCACAGTATAACGGTCTGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGATCGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCa                                                                                                                                 >  1:1487048/1‑139 (MQ=255)  gCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCat                                                                                                                                >  2:913431/1‑139 (MQ=255)      ttACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATgg                                                                                                                            <  2:1267019/139‑1 (MQ=255)      ttACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATgg                                                                                                                            <  1:913431/139‑1 (MQ=255)                   aaCGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCcg                                                                                                               <  1:690570/139‑1 (MQ=255)                        tCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCgtcgtc                                                                                                          >  1:91529/1‑139 (MQ=255)                          ggCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGc                                                                                                        >  2:1547900/1‑139 (MQ=255)                                         gTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGc                                                                                                        >  2:156121/1‑124 (MQ=255)                                         gTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGc                                                                                                        <  1:156121/124‑1 (MQ=255)                                          ttCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATgg                                                                                                                            <  1:1035556/103‑1 (MQ=255)                                          ttCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATgg                                                                                                                            >  2:1035556/1‑103 (MQ=255)                                          ttCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCa                                                                                        <  1:1547900/139‑1 (MQ=255)                                          ttCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCa                                                                                        >  1:332895/1‑139 (MQ=255)                                          ttCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCa                                                                                        <  2:332895/139‑1 (MQ=255)                                          ttCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCa                                                                                        >  1:554629/1‑139 (MQ=255)                                            cAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTg                                                                                                >  2:720669/1‑129 (MQ=255)                                            cAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTg                                                                                                <  1:720669/129‑1 (MQ=255)                                                 gTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGa                                                                                 <  2:91529/139‑1 (MQ=255)                                                 gTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGa                                                                                 <  2:1487048/139‑1 (MQ=255)                                                                   ggCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAt                                                                  >  1:581466/1‑136 (MQ=255)                                                                   ggCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAt                                                                  <  2:581466/136‑1 (MQ=255)                                                                       aaGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACtat                                                           <  2:554629/139‑1 (MQ=255)                                                                               aCTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCat                                                                                                                                >  1:1149754/1‑62 (MQ=255)                                                                               aCTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCat                                                                                                                                <  2:1149754/62‑1 (MQ=255)                                                                               aCTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAt                                                                  <  1:1332940/124‑1 (MQ=255)                                                                               aCTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAt                                                                  >  2:1332940/1‑124 (MQ=255)                                                                                  gAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGa                                                                                 >  1:799979/1‑106 (MQ=255)                                                                                  gAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGa                                                                                 <  2:799979/106‑1 (MQ=255)                                                                                   aGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTc                                                                                         <  1:563921/97‑1 (MQ=255)                                                                                   aGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTc                                                                                         >  2:563921/1‑97 (MQ=255)                                                                                   aGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCt                                               >  1:336380/1‑139 (MQ=255)                                                                                    gCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTa                                              <  2:336380/139‑1 (MQ=255)                                                                                    gCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTa                                              >  1:1262858/1‑139 (MQ=255)                                                                                            gTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGaaa                                      <  2:1262858/139‑1 (MQ=255)                                                                                             ttttGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCa                                                                                                                                 >  2:1359177/1‑47 (MQ=255)                                                                                             ttttGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCa                                                                                                                                 <  1:1359177/47‑1 (MQ=255)                                                                                             ttttGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATc                                                                 <  1:149692/111‑1 (MQ=255)                                                                                             ttttGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATc                                                                 >  2:149692/1‑111 (MQ=255)                                                                                                       cGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGa                                                                                 <  1:1008715/85‑1 (MQ=255)                                                                                                       cGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGa                                                                                 >  2:1008715/1‑85 (MQ=255)                                                                                                           cAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTg                                                                                                >  1:455013/1‑66 (MQ=255)                                                                                                           cAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTg                                                                                                <  2:455013/66‑1 (MQ=255)                                                                                                                tgaattgaTAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACtat                                                           >  2:1438316/1‑98 (MQ=255)                                                                                                                tgaattgaTAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACtat                                                           <  1:1438316/98‑1 (MQ=255)                                                                                                                   attgaTAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAg                                                                                                     >  1:874506/1‑53 (MQ=255)                                                                                                                   attgaTAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAg                                                                                                     <  2:874506/53‑1 (MQ=255)                                                                                                                        tAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGa                                                                                 >  1:1275440/1‑68 (MQ=255)                                                                                                                        tAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGa                                                                                 <  2:1275440/68‑1 (MQ=255)                                                                                                                           aTTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTcc       <  1:380699/139‑1 (MQ=255)                                                                                                                            ttCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGtt                                                                                           <  2:323277/54‑1 (MQ=255)                                                                                                                            ttCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGtt                                                                                           >  1:323277/1‑54 (MQ=255)                                                                                                                            ttCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTccc      >  1:649498/1‑139 (MQ=255)                                                                                                                              cTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTc                                                                                                  <  1:1250639/45‑1 (MQ=255)                                                                                                                              cTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTc                                                                                                  >  2:1250639/1‑45 (MQ=255)                                                                                                                                ggAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGc  <  2:649498/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                   GGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGC  >  NC_000913/4425699‑4425964

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: ‑