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

GCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGCCGTA  >  NC_000913/4425700‑4425968                                                                                                                                       |                                                                                                                                       gCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCat                                                                                                                                    <  1:662721/139‑1 (MQ=255)               tataACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGACACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTAttt                                                                                                                      >  2:1251838/1‑139 (MQ=255)                       tCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGc                                                                                                                               <  2:862356/122‑1 (MQ=255)                       tCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGc                                                                                                                               >  1:862356/1‑122 (MQ=255)                       tCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCgtcgtc                                                                                                              >  1:993844/1‑139 (MQ=255)                       tCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCgtcgtc                                                                                                              >  1:1256445/1‑139 (MQ=255)                       tCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCgtcgtc                                                                                                              >  2:16685/1‑139 (MQ=255)                         ggCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGc                                                                                                            >  1:90073/1‑139 (MQ=255)                                           cAGGGGTTGATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCGGCCGTTTCAcc                                                                                          >  1:853667/1‑139 (MQ=255)                                                   catttcatGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAgg                                                                                  >  2:375496/1‑139 (MQ=255)                                                        catGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTc                                                                             >  2:647447/1‑139 (MQ=255)                                                               caAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAt                                                                      <  1:16685/139‑1 (MQ=255)                                                                aaGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATc                                                                     <  1:1221646/139‑1 (MQ=255)                                                                aaGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATc                                                                     <  1:375496/139‑1 (MQ=255)                                                                aaGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATc                                                                     >  2:605343/1‑139 (MQ=255)                                                                  ggCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGa                                                                   <  2:90073/139‑1 (MQ=255)                                                                       aGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACtata                                                              >  1:194082/1‑139 (MQ=255)                                                                             gACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACg                                                                                                                            <  1:612053/71‑1 (MQ=255)                                                                             gACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACg                                                                                                                            <  1:769088/71‑1 (MQ=255)                                                                             gACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACg                                                                                                                            >  2:612053/1‑71 (MQ=255)                                                                             gACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACg                                                                                                                            >  2:769088/1‑71 (MQ=255)                                                                                 gAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGc                                                                                                            >  1:495336/1‑83 (MQ=255)                                                                                 gAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGc                                                                                                            <  2:495336/83‑1 (MQ=255)                                                                                  aGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCt                                                   >  2:729349/1‑139 (MQ=255)                                                                                  aGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCt                                                   >  2:727414/1‑139 (MQ=255)                                                                                   gCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACGGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTa                                                  >  2:233095/1‑139 (MQ=255)                                                                                      aaaTGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGc                                               <  2:1256445/139‑1 (MQ=255)                                                                                        aTGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAt                                                                      <  1:390671/114‑1 (MQ=255)                                                                                        aTGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAt                                                                      >  2:390671/1‑114 (MQ=255)                                                                                           gTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGaaa                                          <  1:647447/139‑1 (MQ=255)                                                                                            ttttGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGaaaa                                         >  1:979837/1‑139 (MQ=255)                                                                                                     ccGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGaaaaa                                        <  2:274452/131‑1 (MQ=255)                                                                                                     ccGAGCAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGaaaaa                                        >  1:274452/1‑131 (MQ=255)                                                                                                          cAGATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGaa                           <  2:979837/139‑1 (MQ=255)                                                                                                            gATTGAATTGATAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAg                         <  2:61894/139‑1 (MQ=255)                                                                                                               tgaattgaTAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCgg                      <  2:925270/139‑1 (MQ=255)                                                                                                                     gaTAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAg                                                                                   >  1:915882/1‑72 (MQ=255)                                                                                                                     gaTAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAg                                                                                   <  2:915882/72‑1 (MQ=255)                                                                                                                       tAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGa                            >  2:220561/1‑125 (MQ=255)                                                                                                                       tAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGa                            <  1:220561/125‑1 (MQ=255)                                                                                                                       tAGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTg              >  2:585347/1‑139 (MQ=255)                                                                                                                        aGATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGt             <  2:194082/139‑1 (MQ=255)                                                                                                                         gATTCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTc            >  2:210603/1‑139 (MQ=255)                                                                                                                           ttCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTg                                             <  2:1246887/104‑1 (MQ=255)                                                                                                                           ttCTGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTg                                             >  1:1246887/1‑104 (MQ=255)                                                                                                                              tGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACtata                                                              >  1:1195170/1‑84 (MQ=255)                                                                                                                              tGGAGACTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACtata                                                              <  2:1195170/84‑1 (MQ=255)                                                                                                                                   aCTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTc                                                                             <  1:98258/64‑1 (MQ=255)                                                                                                                                   aCTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTc                                                                             >  2:98258/1‑64 (MQ=255)                                                                                                                                   aCTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGCCGTa  <  1:729349/139‑1 (MQ=255)                                                                                                                                   aCTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGCCGTa  >  1:333661/1‑139 (MQ=255)                                                                                                                                   aCTAACCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGCCGTa  <  1:233095/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                       GCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGCCGTA  >  NC_000913/4425700‑4425968

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