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%	155.2 / NA	43	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 (22/21);  total (22/21)

ACCAACAGAATGGTCATAAATTGATTACGGAAGAAATTTCCGGTGTCACGGTATACGGACTGCGCCGTGATAGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCAGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGCGTAGCGAAAATATTTAAAAATTGATTTAAATCACATTA  >  NC_000913/1313587‑1313875                                                                                                                                                   |                                                                                                                                               aCCAACAGAATGGTCATAAATTGATTACGGAAGAAATTTCCGGTGTCACGGTATACGGACTGCGCCGTGATAGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGgcggca                                                                                                                                             <  1:423678/150‑1 (MQ=255)  ccAACAGAATGGTCATAAATTGATTACGGAAGAAATTTCCGGTGTCACGGTATACGGACTGCGCCGTGATAGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATc                                                                                                                                           <  1:446093/151‑1 (MQ=255)      cAGAATGGTCATAAATTGATTACGGAAGAAATTTCCGGTGTCACGGTATACGGACTGCGCCGTGATAGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGt                                                                                                                                       <  2:202146/151‑1 (MQ=255)       aGAATGGTCATAAATTGATTACGGAAGAAATTTCCGGTGTCACGGTATACGGACTGCGCCGTGATAGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTc                                                                                                                                      <  1:409527/151‑1 (MQ=255)               cATAAATTGATTACGGAAGAAATTTCCGGTGTCACGGTATACGGACTGCGCCGTGATAGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATgg                                                                                                                              <  1:488287/151‑1 (MQ=255)                aTAAATTGATTACGGAAGAAATTTCCGGTGTCACGGTATACGGACTGCGCCGTGATAGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGa                                                                                                                             <  2:63531/151‑1 (MQ=255)                  aaaTTGATTACGGAAGAAATTTCCGGTGTCACGGTATACGGACTGCGCCGTGATAGACAAGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGaa                                                                                                                            >  1:164003/1‑150 (MQ=255)                         ttACGGAAGAAATTTCCGGTGTCACGGTATACGGACTGCGCCGTGATAGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCatat                                                                                                                    >  1:114837/1‑151 (MQ=255)                          tACGGAAGAAATTTCCGGTGTCACGGTATACGGACTGCGCCGTGATAGACATGCACTCTCCTTGAGTTTTATGCATGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCatat                                                                                                                    <  1:386082/150‑1 (MQ=255)                               aagaaATTTCCGGTGTCACGGTATACGGACTGCGCCGTGATAGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGt                                                                                                              <  2:512012/151‑1 (MQ=255)                                               cACGGTATACGGACTGCGCCGTGATAGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTg                                                                                              >  1:518874/1‑151 (MQ=255)                                                    tataCGGACTGCGCCGTGATAGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTgg                                                                                          >  1:371960/1‑150 (MQ=255)                                                    tataCGGACTGCGCCGTGATAGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGc                                                                                         >  1:338491/1‑151 (MQ=255)                                                                          aCATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATcgcg                                                                   >  1:763473/1‑151 (MQ=255)                                                                            aTGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGc                                                                 <  2:407340/151‑1 (MQ=255)                                                                             tGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGc                                                                 >  1:368789/1‑150 (MQ=255)                                                                             tGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGc                                                                 >  1:51340/1‑150 (MQ=255)                                                                                aCTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAg                                                              >  1:58025/1‑150 (MQ=255)                                                                                     ccTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTgg                                                           <  1:611228/148‑1 (MQ=255)                                                                                     ccTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGc                                                          <  1:126338/149‑1 (MQ=255)                                                                                     ccTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCg                                                         >  1:512383/1‑150 (MQ=255)                                                                                        tGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGAc                                                       <  2:409632/149‑1 (MQ=255)                                                                                             tttATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGtt                                                 <  1:253467/150‑1 (MQ=255)                                                                                                  gCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATAtt                                           >  1:165568/1‑151 (MQ=255)                                                                                                     ggTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGc                                        <  1:296399/151‑1 (MQ=255)                                                                                                          tAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGCGTAg                                    <  2:673288/150‑1 (MQ=255)                                                                                                          tAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGCGTAGc                                   >  1:121003/1‑151 (MQ=255)                                                                                                                gCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGCGTAGCGAAAat                             <  1:525988/151‑1 (MQ=255)                                                                                                                gCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGCGTAGCGAAAat                             <  1:137407/151‑1 (MQ=255)                                                                                                                 cGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGCGTAGCGAAAat                             >  1:167061/1‑150 (MQ=255)                                                                                                                 cGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGCGTAGCGAAAat                             >  2:286208/1‑150 (MQ=255)                                                                                                                   ggCAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGCGTAGCGAAAATAtt                          <  2:455310/151‑1 (MQ=255)                                                                                                                   ggCAATTGTACCCTGGTTAAGCCTCGGGCGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGCGTAGCGAAAATAtt                          >  2:592574/1‑151 (MQ=255)                                                                                                                     cAATTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGCGTAGCGAAAATATTTa                        <  1:194810/151‑1 (MQ=255)                                                                                                                       aTTGTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGCGTAGCGAAAATATTTaa                       >  1:121919/1‑150 (MQ=255)                                                                                                                          gTACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGCGTAGCGAAAATATTTaaaaa                    <  2:690446/150‑1 (MQ=255)                                                                                                                           tACCCTGGTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGCGTAGCGAAAATATTTAAAAAtt                  >  2:615638/1‑151 (MQ=255)                                                                                                                                 ggTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGCGTAGCGAAAATATTTAAAAATTGAttt             >  2:241362/1‑150 (MQ=255)                                                                                                                                 ggTTAAGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGCGTAGCGAAAATATTTAAAAATTGAtt              <  1:437240/149‑1 (MQ=255)                                                                                                                                    tAAGCCTAGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGCGTAGCGAAAATATTTAAAAATTGATTTAAAt         >  1:571016/1‑151 (MQ=255)                                                                                                                                     aaGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGCGTAGCGAAAATATTTAAAAATTGATTTAAAt         >  2:754054/1‑150 (MQ=255)                                                                                                                                     aaGCCTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGCGTAGCGAAAATATTTAAAAATTGATTTAAATc        >  2:122516/1‑151 (MQ=255)                                                                                                                                        ccTCGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGCGTAGCGAAAATATTTAAAAATTGATTTAAATcac      >  1:405231/1‑150 (MQ=255)                                                                                                                                           cGGGTGGCGGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGCGTAGCGAAAATATTTAAAAATTGATTTAAATCACATTa  <  2:109115/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                               ACCAACAGAATGGTCATAAATTGATTACGGAAGAAATTTCCGGTGTCACGGTATACGGACTGCGCCGTGATAGACATGCACTCTCCTTGAGTTTTATGCAGGTGTTAATTAGCGGGCAATTGTACCCTGGTTAAGCCTCGGGTGGCAGCATCAGGTCTTGTATGGAAAAGCATATCTTTGTAAAGCGGAGGTAATTGCTGGCCTTCAGGCATGGGCAAATCGCGGCAAGTGGCGACCGAGTTAATATTTGCGTAGCGAAAATATTTAAAAATTGATTTAAATCACATTA  >  NC_000913/1313587‑1313875

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