Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,313,733	A→G	intergenic (‑70/‑287)	yciC ← / → ompW	UPF0259 family inner membrane protein/outer membrane protein W

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,313,733	0	A	G	100.0%	150.0 / NA	42	intergenic (‑70/‑287)	yciC/ompW	UPF0259 family inner membrane protein/outer membrane protein W
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/22);  total (20/22)

ACAGAATGGTCATAAATTGATTACGGAAGAAATTTCCGGTGTCACGGTATACGGACTGCGCCGTGATAGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCAGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGCGTAGCGAAAATATTTAAAAATTGATTTAAATCACATTAACCAG  >  NC_000913/1313591‑1313880                                                                                                                                               |                                                                                                                                                    aCAGAATGGTCATAAATTGATTACGGAAGAAATTTCCGGTGTCACGGTATACGGACTGCGCCGTGATAGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAgg                                                                                                                                             >  1:322102/1‑151 (MQ=255)   aGAATGGTCATAAATTGATTACGGAAGAAATTTCCGGTGTCACGGTATACGGACTGCGCCGTGATAGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTc                                                                                                                                           <  1:700236/151‑1 (MQ=255)    gAATGGTCATAAATTGATTACGGAAGAAATTTCCGGTGTCACGGTATACGGACTGCGCCGTGATAGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCt                                                                                                                                          <  2:218114/151‑1 (MQ=255)        ggTCATAAATTGATTACGGAAGAAATTTCCGGTGTCACGGTATACGGACTGCGCCGTGATAGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTa                                                                                                                                      <  1:452554/151‑1 (MQ=255)            aTAAATTGATTACGGAAGAAATTTCCGGTGTCACGGTATACGGACTGCGCCGTGATAGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATgg                                                                                                                                   >  2:686298/1‑150 (MQ=255)            aTAAATTGATTACGGAAGAAATTTCCGGTGTCACGGTATACGGACTGCGCCGTGATAGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGa                                                                                                                                  <  1:224474/151‑1 (MQ=255)             tAAATTGATTACGGAAGAAATTTCCGGTGTCACGGTATACGGACTGCGCCGTGATAGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGa                                                                                                                                  >  2:258777/1‑150 (MQ=255)              aaaTTGATTACGGAAGAAATTTCCGGTGTCACGGTATACGGACTGCGCCGTGATAGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGaa                                                                                                                                 >  2:562561/1‑150 (MQ=255)               aaTTGATTACGGAAGAAATTTCCGGTGTCACGGTATACGGACTGCGCCGTGATAGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGaaa                                                                                                                                >  2:665710/1‑150 (MQ=255)                aTTGATTACGGAAGAAATTTCCGGTGTCACGGTATACGGACTGCGCCGTGATAGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGaaa                                                                                                                                >  2:906555/1‑149 (MQ=255)                  tGATTACGGAAGAAATTTCCGGTGTCACGGTATACGGACTGCGCCGTGATAGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGc                                                                                                                             >  2:721124/1‑150 (MQ=255)                   gATTACGGAAGAAATTTCCGGTGTCACGGTATACGGACTGCGCCGTGATAGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCa                                                                                                                            <  1:360126/150‑1 (MQ=255)                     ttACGGAAGAAATTTCTGGTGTCACGGGATACGGACTGCGCCGTGATAGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCatat                                                                                                                         >  2:703858/1‑151 (MQ=255)                             gaaATTTCCGGTGTCACGGTATACGGACTGCGCCGTGATAGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTaa                                                                                                                 <  2:748281/151‑1 (MQ=255)                             gaaATTTCCGGTGTCACGGTATACGGACTGCGCCGTGATAGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTaa                                                                                                                 <  1:293274/151‑1 (MQ=255)                             gaaATTTCCGGTGTCACGGTATACGGACTGCGCCGTGATAGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTa                                                                                                                  >  1:52554/1‑150 (MQ=255)                                  ttCCGGTGTCACGGTATACGGACTGCGCCGTGATAGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCg                                                                                                             <  2:694115/150‑1 (MQ=255)                                       gtgtCACGGTATACGGACTGCGCCGTGATAGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTa                                                                                                       <  2:150018/151‑1 (MQ=255)                                          tCACGGTATACGGACTGCGCCGTGATAGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAAtt                                                                                                    >  2:687478/1‑151 (MQ=255)                                           cACGGTATACGGACTGCGCCGTGATAGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAAtt                                                                                                    <  1:334368/150‑1 (MQ=255)                                            aCGGTATACGGACTGCGCCGTGATAGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAAtt                                                                                                    <  2:741179/149‑1 (MQ=255)                                                               gTGATAGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGc                                                                                >  1:779832/1‑150 (MQ=255)                                                                    aGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATc                                                                           <  2:558251/150‑1 (MQ=255)                                                                     gACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATcgc                                                                         <  1:907882/151‑1 (MQ=255)                                                                            aCTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGt                                                                  <  1:148931/151‑1 (MQ=255)                                                                                         tttATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGtt                                                      >  1:834783/1‑150 (MQ=255)                                                                                          ttATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTa                                                     >  1:204075/1‑150 (MQ=255)                                                                                               cAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATAttt                                               <  1:627120/151‑1 (MQ=255)                                                                                                aGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTg                                              >  2:757824/1‑151 (MQ=255)                                                                                                 ggTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGc                                             <  1:492194/151‑1 (MQ=255)                                                                                                        aTTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGCGTAGCg                                       >  2:899929/1‑150 (MQ=255)                                                                                                             cGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGCGTAGCGAAAata                                 <  1:719740/151‑1 (MQ=255)                                                                                                              gggCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGCGTAGCGAAAatat                                >  2:303245/1‑151 (MQ=255)                                                                                                                 cAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGCGTAGCGAAAATAttt                              >  2:549511/1‑150 (MQ=255)                                                                                                                     tGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGCGTAGCGAAAATATTTaaaa                          <  2:50911/150‑1 (MQ=255)                                                                                                                         cccTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGCGTAGCGAAAATATTTAAAAATTg                      >  1:174395/1‑150 (MQ=255)                                                                                                                             ggTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGCGTAGCGAAAATATTTAAAAATTGAttt                  <  1:651723/150‑1 (MQ=255)                                                                                                                             ggTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGCGTAGCGAAAATATTTAAAAATTGATTTa                 <  1:428944/151‑1 (MQ=255)                                                                                                                              gTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGCGTAGCGAAAATATTTAAAAATTGATTTa                 >  1:417632/1‑150 (MQ=255)                                                                                                                                    ccTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGCGTAGCGAAAATATTTAAAAATTGATTTAAATcaca          <  1:427825/151‑1 (MQ=255)                                                                                                                                       cGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGCGTAGCGAAAATATTTAAAAATTGATTTAAATCACATTa       <  1:902137/151‑1 (MQ=255)                                                                                                                                        gggTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGCGTAGCGAAAATATTTAAAAATTGATTTAAATCACATTa       >  2:471710/1‑150 (MQ=255)                                                                                                                                            ggcggcATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGCGTAGCGAAAATATTTAAAAATTGATTTAAATCACATTAACCAg  <  1:731640/151‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                    ACAGAATGGTCATAAATTGATTACGGAAGAAATTTCCGGTGTCACGGTATACGGACTGCGCCGTGATAGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCAGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGCGTAGCGAAAATATTTAAAAATTGATTTAAATCACATTAACCAG  >  NC_000913/1313591‑1313880

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