Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,433,164	C→T	*287* (TAG→TAA)	qorB ←	NAD(P)H:quinone oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,433,164	0	C	T	100.0%	115.0 / NA	37	*287* (TAG→TAA)	qorB	NAD(P)H:quinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/17);  total (20/17)

GGCGTACAGCACGCCGCTAATCATTACTGGCATCCATTGATCTGTCAGAAGAAAGTTCAGGAAGCAGACATGATTGCTCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAACTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTAACGTTGTCGTTGGGTGGCCAATCAATTTGCTAAGTGTTTTGCTGTCATCAAACAGACCGCCTTTCGATGCGCCAACGTCAG  >  NC_000913/4433033‑4433285                                                                                                                                    |                                                                                                                          ggCGTACAGCACGCCGCTAATCATTACTGGCATCCATTGATCTGTCAGAAGAAAGTTCAGGAAGCAGACATGATTGCTCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGttat                                                                                                                    <  2:711036/139‑1 (MQ=255)  gCGTACAGCACGCCGCTAATCATTACTGGCATCCATTGATCTGTCAGAAGAAAGTTCAGGTAGCAGACATGATTGCTCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGttatt                                                                                                                   <  2:1015331/139‑1 (MQ=255)   cGTACAGCACGCCGCTAATCATTACTGGCATCCATTGATCTGTCAGAAGAAAGTTCAGGAAGCAGACATGATTGCTCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGttatta                                                                                                                  <  1:1272660/139‑1 (MQ=255)            cgccgcTAATCATTACTGGCATCCATTGATCTGTCAGAAGAAAGTTCAGGAAGCAGACATGATTGCTCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTaaaa                                                                                                         >  2:1294117/1‑139 (MQ=255)              ccgcTAATCATTACTGGCATCCATTGATCTGTCAGAAGAAAGTTCAGGAAGCAGACATGATTGCTCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAg                                                                                                       <  2:991516/139‑1 (MQ=255)               cgcTAATCATTACTGGCATCCATTGATCTGTCAGAAGAAAGTTCAGGAAGCAGACATGATTGCTCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGa                                                                                                      >  1:476062/1‑139 (MQ=255)                    aTCATTACTGGCATCCATTGATCTGTCAGAAGAAAGTTCAGGAAGCAGACATGATTGCTCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGttat                                                                                                                    >  1:197927/1‑120 (MQ=255)                    aTCATTACTGGCATCCATTGATCTGTCAGAAGAAAGTTCAGGAAGCAGACATGATTGCTCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGttat                                                                                                                    <  2:197927/120‑1 (MQ=255)                                aTCCATTGATCTGTCAGAAGAAAGTTCAGGAAGCAGACATGATTGCTCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGc                                                                                     <  2:395266/139‑1 (MQ=255)                                      tGATCTGTCAGAAGAAAGTTCAGGAAGCAGACATGATTGCTCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTAACGt                                                                               >  2:559689/1‑139 (MQ=255)                                        aTCTGTCAGAAGAAAGTTCAGGAAGCAGACATGATTGCTCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTAACGTTg                                                                             >  2:831370/1‑139 (MQ=255)                                             tCAGAAGAAAGTTCAGGAAGCAGACATGATTGCTCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTAACGTTGTCGtt                                                                        >  1:603328/1‑139 (MQ=255)                                                       gTTCAGGAAGCAGACATGATTGCTCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTAACGTTGTCGTTGGGTGGCcaa                                                              >  1:115441/1‑139 (MQ=255)                                                            ggAAGCAGACATGATTGCTCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTAACGTTGTCGTTGGGTGGCcaatcaat                                                         <  1:831370/139‑1 (MQ=255)                                                                         aTTGCTCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGttat                                                                                                                    >  1:689977/1‑67 (MQ=255)                                                                         aTTGCTCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGttat                                                                                                                    <  2:689977/67‑1 (MQ=255)                                                                           tGCTCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTAACGTTGTCGTTGGGTGGCCAATCAATTTGCTAAGTGTTTTg                                          >  2:1134926/1‑139 (MQ=255)                                                                             cTCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCAGCTAACGTTGTCGTTGGGTGGCCAATCAATTTGCTAAGTGTTTTGCt                                        >  1:787006/1‑139 (MQ=255)                                                                              tCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTAACGTTGTCGTTGGGTGGcc                                                                <  2:964365/114‑1 (MQ=255)                                                                              tCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTAACGTTGTCGTTGGGTGGcc                                                                >  1:964365/1‑114 (MQ=255)                                                                              tCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTAACGTTGTCGTTGGGTGGCCAATCAATTTGCTAAGTGTTTTGCTg                                       >  1:903548/1‑139 (MQ=255)                                                                              tCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTAACGTTGTCGTTGGGTGGCCAATCAATTTGCTAAGTGTTTTGCTg                                       >  1:1138976/1‑139 (MQ=255)                                                                              tCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTAACGTTGTCGTTGGGTGGCCAATCAATTTGCTAAGTGTTTTGCTg                                       <  1:1159752/139‑1 (MQ=255)                                                                                   gCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTAACGTTGTCGTTGGGTGGCCAATCAATTTGCTAAGTGTTTTGCTGtcatc                                  <  1:1294117/139‑1 (MQ=255)                                                                                       cccGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTAACGTTGTCGTTGGGTGGcc                                                                <  2:1521929/105‑1 (MQ=255)                                                                                       cccGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTAACGTTGTCGTTGGGTGGcc                                                                >  1:1521929/1‑105 (MQ=255)                                                                                           tCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTAACGTTGTCGTTGGGTGGCCAATCAATTTGCTAAGTGTTTTGCTGTCATCAAACAGAc                          <  2:903548/139‑1 (MQ=255)                                                                                                     agagGATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCt                                                                                    <  2:1273423/71‑1 (MQ=255)                                                                                                     agagGATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCt                                                                                    >  1:1273423/1‑71 (MQ=255)                                                                                                         gATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTAACGTTGTCGTTGGGTGGCCAATCAATTTGCt                                                    >  2:1097930/1‑99 (MQ=255)                                                                                                         gATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTAACGTTGTCGTTGGGTGGCCAATCAATTTGCt                                                    <  1:1097930/99‑1 (MQ=255)                                                                                                         gATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTAACGTTGTCGTTGGGTGGCCAATCAATTTGCTAAGTGTTTTGCTGTCATCAAACAGACCGCCTTTCGAt               <  2:632175/136‑1 (MQ=255)                                                                                                         gATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTAACGTTGTCGTTGGGTGGCCAATCAATTTGCTAAGTGTTTTGCTGTCATCAAACAGACCGCCTTTCGAt               >  1:632175/1‑136 (MQ=255)                                                                                                           tGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTAACGTTGTCGTTGGGTGGCCAATCAATTTGCTAAGTGTTTTGCTGTCATCAAACAGACCGCCTTTCGATGCGcc          >  2:184773/1‑139 (MQ=255)                                                                                                                   gaTGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTAACGTTGTCGTTGGGTGGCCAATCAATTTGCTAAGTGTTTTGCTGTCATCAAACAGACCGCCTTTCGATGCGCCAACGTCAg  >  2:608783/1‑139 (MQ=255)                                                                                                                      gCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTAACGTTGTCGTTGGGTGGCcaatcaat                                                         >  2:22155/1‑81 (MQ=255)                                                                                                                      gCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTAACGTTGTCGTTGGGTGGCcaatcaat                                                         <  1:22155/81‑1 (MQ=255)                                                                                                                                    |                                                                                                                          GGCGTACAGCACGCCGCTAATCATTACTGGCATCCATTGATCTGTCAGAAGAAAGTTCAGGAAGCAGACATGATTGCTCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAACTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTAACGTTGTCGTTGGGTGGCCAATCAATTTGCTAAGTGTTTTGCTGTCATCAAACAGACCGCCTTTCGATGCGCCAACGTCAG  >  NC_000913/4433033‑4433285

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