Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,491,105	C→T	T68T (ACC→ACT)	cybB →	superoxide oxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,491,105	0	C	T	100.0%	81.8 / NA	29	T68T (ACC→ACT)	cybB	superoxide oxidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/17);  total (12/17)

CCCGTTATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAAATATTCAAGGTTACAAATCAGCATTCACTGGCTGGTCTTTTTACTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCG  >  NC_000913/1490860‑1491324                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                            cccGTTATGTTAACGGAAACTTATGTTAAAGCAGGAAATGGTATGGACAATAAATATTCAAGGTTACAAATCAGCATTCACTGGCTGGTCTTTTTACTGGTTATCACAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAAAATGATTCATGTTTCCTGTGGCATCTCCATTCTCGTGCTGATGGTCATTAGTCTTTTGTTTAAGCTGAAAATCCCAACTccgcc                                                                                                                                                                                                                        >  2:518195/1‑251 (MQ=255) cccGTTATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAAATATTCAAGGTTACAAATCAGCATTCACTGGCTGGTCTTTTTACTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACTccgcc                                                                                                                                                                                                                        <  1:569828/251‑1 (MQ=255)            aaCGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAAATATTCAAGGTTACAAATCAGCATTCACTGGCTGGTCTTTTTACTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTAGTTCGTCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGTTTCTACCt                                                                                                                                                                                                              >  2:550302/1‑250 (MQ=255)              cGGAAAATTAAGTTAAAGCAAGAAATGGTATTGAAAATAAATATTCAAGGTTACAAATCAGCATTCACTGGCTGGTCTTTTTACTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCTATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGATTATACCTAAg                                                                                                                                                                                                           <  1:79527/251‑1 (MQ=255)                       aTGTTAAAGCAGGAAATGTTATGGAAAATAAATATTCAAGGTTACAAATCAGCATTCACTGGCTGGTCTTTTTACTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGATTATACCTAAGCCAAAACCg                                                                                                                                                                                                  >  1:539135/1‑251 (MQ=255)                                                    aaataTTCAAGGTTACAAATCAGCATTCACTGGCTGGTCTTTTTACTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGAc                                                                                                                                                                      <  1:461555/250‑1 (MQ=255)                                                                    aaaTCAGCATTCACTGGCTGGTCTTTTTACTGGTTATCGCAGAGTATTGCGCAATGGAGGTTCGTGGTTTCTTCCCACGTAGCGATCTGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGGCGTTCGTCTTTTGTTAAGGATGAAATAACCAACTCCCCCGTTTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTTGTGATTTAtatc                                                                                                                                                     >  2:153738/1‑248 (MQ=255)                                                                      aTCAGCATTCACTGGCTGGTCTTTTTACTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCt                                                                                                                                                   <  1:550302/251‑1 (MQ=255)                                                                            aTTCACTGGCTGGTCTTTTTACTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACAt                                                                                                                                                                    <  1:242335/228‑1 (MQ=255)                                                                            aTTCACTGGCTGGTCTTTTTACTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCAATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGTTGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACAtagatcggtagcgcggcgtgtagg                                                                                                                                             >  2:242335/1‑228 (MQ=255)                                                                               cACTGGCTGGTCTTTTTACTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGt                                                                                                                                                               <  1:17655/230‑1 (MQ=255)                                                                                      tggTATTTTTACTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCACACGTAGCGCTCGGCCAATCATCAACATGCTTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGCAATACCCAACTCCGCCGATTATACCTACTCCAGAACCGATGATGACGGGACTGGCGCATGTGGTACGTTTGGTTGTTTATCTCCTGTTTATTGCGCTGcc                                                                                                                                    >  2:99900/1‑250 (MQ=255)                                                                                                       aTCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGtga                                                                                                                   <  1:324082/250‑1 (MQ=255)                                                                                                              gttATTGCGCAATGGAGTTTCGTTGTTTCTTCCCACTTAGCTATCGGCCCCTCATCAACATGAGTCATTTTTCCTTTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATGTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGtat                                                                                                           <  2:230701/249‑1 (MQ=255)                                                                                                                   tGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCCGTGATTGGTTTGGTGATGATGTATAACCg                                                                                                      >  1:45136/1‑251 (MQ=255)                                                                                                                     cgcAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCggg                                                                                                    <  1:112427/251‑1 (MQ=255)                                                                                                                            gAGTTTCGTGGTTGCTTCCCACGTAGCGATCGGCCACTCACCAACATGATTCATGGTTCCTTTGGCATCTCAATTCTCTTGCTGAGGGTCGTTCGTCTTGTGTTAAGGCGGAAATACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAAccc                                                                                             <  2:212385/251‑1 (MQ=255)                                                                                                                                                  gTAGCGATCGGCCACTCATCAACAGGATTCATGTTTCCTGGGGCATCACAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAGCTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGGTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACg                                                                       <  1:90390/251‑1 (MQ=255)                                                                                                                                                    aGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCCGTGATTGGTTTGGTGATGATGTATAACCGGGGCACCCCGTGGTTTGCGTTTGGTTTGACGAt                                                                     >  2:50943/1‑251 (MQ=255)                                                                                                                                                                             ttCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAAttt                                             <  1:50943/250‑1 (MQ=255)                                                                                                                                                                                             cATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCt                            <  1:99900/251‑1 (MQ=255)                                                                                                                                                                                                  cAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAg                       <  1:518195/251‑1 (MQ=255)                                                                                                                                                                                                         cGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTgg                 <  1:399374/250‑1 (MQ=255)                                                                                                                                                                                                          ttGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGATATAACCAACTCCGCCGATTATACCTCAGCCAAAACCGATGATGACGGGCCTGGCGCATTTGGGCCATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGGTTGGTTTGGGGATGATGTATAACCGGGACAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTACGAGGCCAATTTCGAACGGGTAGATGGCTTAAAGTCGTGGc                >  2:343171/2‑250 (MQ=255)                                                                                                                                                                                                               gATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGaaa          <  1:486056/251‑1 (MQ=255)                                                                                                                                                                                                               gATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGaa           >  1:164220/1‑250 (MQ=255)                                                                                                                                                                                                                        tcgtCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCg  >  1:477600/1‑250 (MQ=255)                                                                                                                                                                                                                                              tACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACg        >  1:569667/1‑222 (MQ=255)                                                                                                                                                                                                                                              tACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACg        <  2:569667/222‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                            CCCGTTATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAAATATTCAAGGTTACAAATCAGCATTCACTGGCTGGTCTTTTTACTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCG  >  NC_000913/1490860‑1491324

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