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	96.7%	84.7 / NA	30	T68T (ACC→ACT)	cybB	superoxide oxidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (18/11);  minor base G (0/1);  total (18/12)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GTTATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAAATATTCAAGGTTACAAATCAGCATTCACTGGCTGGTCTTTTTACTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCAT  >  NC_000913/1490863‑1491344                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                gTTATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAAATATTCAAGGTTACAAATCAGCATTCACTGGCTGGTCTTTTTACTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGAt                                                                                                                                                                                                                                         >  2:140843/1‑251 (MQ=255)                aaTTATGTTAAAGCAGGAAAGGTTATGGAAAATAAATATTCAAGGTTACAAATCAGCATTCACTGGCTGGTCTTTTTACTGTTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTGCCCACGTAGCGATCGGCCACTCATCAACATGAGTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCATCTCCGTCGATTATACCTAAGCCaaa                                                                                                                                                                                                                          <  2:26328/251‑1 (MQ=255)                             cAGGAAATGTTATGGAAAATAAATATTCAAGGTTACAAATCAGCATTCACTGGCTGGTCTTTTTACTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACg                                                                                                                                                                                                             <  1:405319/251‑1 (MQ=255)                                   aTGTTATGGAAAATAAATATTCAAGGTTCCAAATCAGCATTCACTGGCTGGTCTTTTTACTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGGCGTTCGTCTTTTGTTAAGGCTGAAATACCCAAATCCCCCGATTATACCTAAGCCAAAACCGATGATGACgggcc                                                                                                                                                                                                         >  2:375515/1‑247 (MQ=255)                                   aTGTTATGGAAAATAAATATTCAAGGTTACAAATCAGCATTCACTGGCTGGTCTTTTTACTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTg                                                                                                                                                                                                       <  2:202401/251‑1 (MQ=255)                                    tGTTATGGAAAATAAATATTCAAGGTTACAAATCAGCATTCACTGGCTGGTCTTTTTACTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTgg                                                                                                                                                                                                      >  1:229214/1‑251 (MQ=255)                                           gAAAATAAATATTCAAGGTTACAAATCAGCATTCACTGGCTGGTCTTTTTACTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCCttt                                                                                                                                                                                               >  2:18484/1‑251 (MQ=255)                                                    tatTCAAGGTTACAAATCAGCATTCACTGGCTGGTCTTTTTACTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACAttt                                                                                                                                                                                      >  1:397425/1‑251 (MQ=255)                                                                     cAGCAATCACTGGCTGGTCTTTGTACTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGTTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTGGTTAAGGCTGAAATACCCAACTCCGCCGATTATGCCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGt                                                                                                                                                                     <  2:471445/251‑1 (MQ=255)                                                                            cACTGGCTGGTCTTTTTACTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTgc                                                                                                                                                              >  1:32133/1‑251 (MQ=255)                                                                                  ctggTCTTTTTACTGGTTATCGCCGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGcc                                                                                                                                                        >  2:118878/1‑251 (MQ=255)                                                                                     gTCTTTTTACTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGt                                                                                                                                                     >  2:571475/1‑251 (MQ=255)                                                                                     gTCTTTTTACTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGt                                                                                                                                                     >  2:474378/1‑251 (MQ=255)                                                                                           ttACTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAttg                                                                                                                                                <  1:420218/250‑1 (MQ=255)                                                                                              cTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGATTATACCTAAGCCACAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATATCCTGTTTATTGCGCTGCCAGTGAttggttt                                                                                                                                            >  1:364233/1‑251 (MQ=255)                                                                                                ggTTATCGCAGCGTATTGCGCAATTGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACGCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTAACTCCTGTTTATTGCGCTGCCAGTGAttggttt                                                                                                                                            <  2:397425/249‑1 (MQ=255)                                                                                                                  cgcAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCACCAAAATGATTCATGTTTCCTGTGCCATCTCAATTCTCGTGCTGATGGTCGATCGTCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGATTATACCTAAGCCAGAACCGCTGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTACTGTTTATTGCACTGCCAGTGATTGGTTTGGTGATGATGTATAACCggg                                                                                                                        >  1:528349/1‑251 (MQ=255)                                                                                                                    cAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATGCCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGc                                                                                                                      <  1:406438/251‑1 (MQ=255)                                                                                                                                            cACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCgtttggtttg                                                                                              <  1:18484/251‑1 (MQ=255)                                                                                                                                                        ggCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTCTTGTTAAGGTTGAAATGCCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCAATTGGGACATTTGGCGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACAGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTAc                                                                                  >  1:346907/1‑251 (MQ=255)                                                                                                                                                                   aaCATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGCTGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGc                                                                       >  1:265008/1‑251 (MQ=255)                                                                                                                                                                    acatgatgtatttttCCTGTGGCATCTGAAATCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTATGGTTTGACGATGCCTTACGCTTCATAGGcc                                                                      <  1:421033/239‑1 (MQ=255)                                                                                                                                                                        taTTCATGTTTCCTGGGTCATCGCAATTCGCGTGCTGATGGTCGTTCTTCTTTTGTTAAGGCTGAAAAACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTTCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGAGGCCGTACGCTTCAGAGGCCAAt                                                                   <  2:364233/249‑1 (MQ=255)                                                                                                                                                                                      gtgGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGAt                                                    <  2:32133/251‑1 (MQ=255)                                                                                                                                                                                            tctcAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCGAACTCCGCCGATTATACATCAGCCAAAACCGATGATGCAGGGACAAGCGGATTAGGTACATTTGGTGATTTATCTCCAGTTTATTGCGCTGCCAGTGATTGGTTCGGTGATGATGCATAACAGGGGCAACCCGGGGTTTGCGTTTTGTTTGACGATGCCTTAAGCTTCAGAGGCAAATTTCGACCGGGGAGATAGCtt                                               >  1:578815/1‑250 (MQ=255)                                                                                                                                                                                               cAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTTTTCCCCTGCCAGGAATTGGTTTGGGGTTGTTGTATAACCGGGCCACCCCGGGGTTTGCGTTGGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAAggggaggttggtttaa                                             >  1:427365/1‑237 (MQ=255)                                                                                                                                                                                                                          ttttGTTAAGGCTGAAATACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTg                >  1:258544/1‑251 (MQ=255)                                                                                                                                                                                                                               ttAAGGCTGAAATACCCAACTCCGCCGCTTATACCTAAGCCAAAACAGATGATGACGGGACTGGCGCCTTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGCTGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGACACGCTGGCGCATCTggggt            >  2:499980/1‑248 (MQ=255)                                                                                                                                                                                                                                       gAAATACCCAACTCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTc    <  1:571475/250‑1 (MQ=255)                                                                                                                                                                                                                                        aaaTACCCAACTACGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCCCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGATTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGATGGCGACTCTGGGATATTTTGTCAt  >  1:285766/1‑251 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                GTTATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAAATATTCAAGGTTACAAATCAGCATTCACTGGCTGGTCTTTTTACTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCAT  >  NC_000913/1490863‑1491344

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