Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	772,414	T→G	R319R (CGT→CGG)	cydA →	cytochrome d terminal oxidase, subunit I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	772,414	0	T	G	97.8%	158.0 / ‑4.0	45	R319R (CGT→CGG)	cydA	cytochrome d terminal oxidase, subunit I
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (26/18);  total (27/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-01

CTGTTCGGCATTCCTGATCAGGAAGAGGAGACGAACAAATTTGCGATTCAGATCCCTTACGCACTGGGCATCATTGCAACGCGTTCCGTGGATACCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGTAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGC  >  NC_000913/772268‑772557                                                                                                                                                   |                                                                                                                                                ctGTTCGGCATTCCTGATCAGGAAGAGGAGACGAACAAATTTGCGATTCAGATCCCTTACGCACTGGGCATCATTGCAACGCGTTCCGTGGATACCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGGAACg                                                                                                                                             >  2:652768/1‑151 (MQ=255)  tGTTCGGCATTCCTGATCAGGAAGAGGAGACGAACAAATTTGCGATTCAGATCCCTTACGCACTGGGCATCATTGCAACGCGTTCCGTGGATACCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGGAACgg                                                                                                                                            >  1:524355/1‑151 (MQ=255)   gTTCGGCATTCCTGATCAGGAAGAGGAGACGAACAAATTTGCGATTCAGATCCCTTACGCACTGGGCATCATTGCAACGCGTTCCGTGGATACCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGGAAc                                                                                                                                              <  2:164912/148‑1 (MQ=255)     tCGGCATTCCTGATCAGGAAGAGGAGACGAACAAATTTTCGATTCAGTTCCCTTACTCTCTGGGCATCATTGCATCGCGTTCCTTGGATACCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGGaa                                                                                                                                               <  1:411965/145‑1 (MQ=255)          aTTCCTGATCAGGAAGAGGAGACGAACAAATTTGCGATTCAGATCCCTTACGCACTGGGCATCATTGCAACGCGTTCCGTGGATACCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGGAACGGGATGAAgg                                                                                                                                    >  1:196179/1‑151 (MQ=255)          aTTCCTGATCAGGAAGAGGAGACGAACAAATTTGCGATTCAGATCCCTTACGCACTGGGCATCATTGCAACGCGTTCCGTGGATACCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGGAACGGGATGAAgg                                                                                                                                    >  2:702867/1‑151 (MQ=255)               tGATCAGGAAGAGGAGACGAACAAATTTGCGATTCAGATCCCTTACGCACTGGGCATCATTGCAACGCGTTCCGTGGATACCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGGAACGGGATGAAGGCGTAc                                                                                                                               >  2:417873/1‑151 (MQ=255)                    aGGAAGAGGAGACGAACAAATTTGCGATTCAGATCCCTTACGCACTGGGCATCATTGCAACGCGTTCCGTGGATACCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGGAACGGGATGAAGGCGTActctct                                                                                                                          >  2:803195/1‑151 (MQ=255)                        agagGAGACGAAAAAATTTGCGATTCAGATCCCTTACGCACTGGGCATCATTGCAACGCGTTCCGTGGATACCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGGAACGGGATGAAGGCGTACTCTCTGCTc                                                                                                                      <  2:80378/151‑1 (MQ=255)                            gagaCGAACAAATTTGCGATTCAGATCCCTTACGCACTGGGCATCATTGCAACGCGTTCCGTGGATACCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGGAACGGGATGAAGGCGTACTCTCTGCTCGaac                                                                                                                  >  1:267196/1‑151 (MQ=255)                              gaCGAACAAATTTGCGATTCAGATCCCTTACGCACTGGGCATCATTGCAACGCGTTCCGTGGATACCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGTAACGGGATGAAGGCGTACTCTCTGCTCGaaca                                                                                                                 >  2:867628/1‑150 (MQ=255)                                 gAACAAATTTGCGATTCAGATCCCTTACGCACTGGGCATCATTGCAACGCGTTCCGTGGATACCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGGAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTg                                                                                                             >  1:471566/1‑151 (MQ=255)                                   aCAAATTTGCGATTCAGATCCCTTACGCACTGGGCATCATTGCAACGCGTTCCGTGGATACCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGGAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCg                                                                                                           <  2:11709/151‑1 (MQ=255)                                     aaaTTTGCGATTCAGATCCCTTACGCACTGGGCATCATTGCAACGCGTTCCGTGGATACCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGGAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGtt                                                                                                         <  2:312519/151‑1 (MQ=255)                                            cGATTCAGATCCCTTACGCACTGGGCATCATTGCAACGCGTTCCGTGGATACCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGGAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTc                                                                                                  <  2:194223/151‑1 (MQ=255)                                                tCAGATCCCTTACGCACTGGGCATCATTGCAACGCGTTCCGTGGATACCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGGAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTAcc                                                                                              >  1:779658/1‑151 (MQ=255)                                                        cTTACGCACTGGGCATCATTGCAACGCGTTCCGTGGATACCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGGAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGc                                                                                      >  2:499857/1‑151 (MQ=255)                                                          tACGCACTGGGCATCATTGCAACGCGTTCCGTGGATACCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGGAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCg                                                                                     <  2:624322/150‑1 (MQ=255)                                                            cGCTCTGGGCATCATTGCAACGCGTTCCGTGGATACCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGGAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCg                                                                                     <  2:533091/148‑1 (MQ=255)                                                            cGCACTGGGCATCATTGCAACGCGTTCCGTGGATACCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGGAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGt                                                                                   >  2:407145/1‑150 (MQ=255)                                                             gCACTGGGCATCATTGCAACGCGTTCCGTGGATACCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGGAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTc                                                                                 >  1:781803/1‑151 (MQ=255)                                                             gCACTGGGCATCATTGCAACGCGTTCCGTGGATACCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGGAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTc                                                                                 <  1:261767/151‑1 (MQ=255)                                                                      atcatTGCAACGCGTTCCGTGGATACCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGGAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGt                                                                        <  2:798119/151‑1 (MQ=255)                                                                          ttGCAACGCGTTCCGTGGATACCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGGAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGCTCGTGACCAGTTCa                                                                     >  2:263856/1‑150 (MQ=255)                                                                            gCAACGCGTTCCGTGGATACCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGGAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAAt                                                                   <  2:413629/150‑1 (MQ=255)                                                                                cgcgTTCCGTGGATACCCCGGTTATCGACCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGGAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCa                                                               >  1:768581/1‑150 (MQ=255)                                                                                   gTTCCGTGGATACCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGGAACGGGATGAAGGCGTACTCTATGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATgaa                                                           <  2:398545/151‑1 (MQ=255)                                                                                        gTGGATACCCCGGCTATCGGCCTGGAAGAGCTGATGGTGCAGCATGAAGCACCCATTCGGAACGGGATGAAGGCGGACCCTCTATTCGAACAACTGCGGTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATTCCATgtagaa                                                        >  1:2418/1‑149 (MQ=255)                                                                                           gATACCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGGAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGAcc                                                   >  2:710532/1‑151 (MQ=255)                                                                                              aCCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGGAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTc                                                 <  1:297112/150‑1 (MQ=255)                                                                                                cccGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGGAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGt                                              >  2:475354/1‑151 (MQ=255)                                                                                                          ggCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGGAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTct                                      >  1:155040/1‑149 (MQ=255)                                                                                                           gCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGGAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTctgct                                   >  1:141582/1‑151 (MQ=255)                                                                                                           gCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGGAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTctgct                                   <  1:320706/151‑1 (MQ=255)                                                                                                               gAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGGAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTctgctgct                                <  2:882762/150‑1 (MQ=255)                                                                                                               gAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGGAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTctgctgct                                <  2:460068/150‑1 (MQ=255)                                                                                                                aaaGAGCTGATGGTGCAGCATGAAGAACGCATTCGGAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTctgctgctg                               >  1:857796/1‑150 (MQ=255)                                                                                                                      cTGATGGTGCAGCATGAAGAACGCATTCGGAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGc                         >  1:110013/1‑150 (MQ=255)                                                                                                                           ggTGCAGCATGAAGAACGCATTCGGAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACg                   <  1:148840/151‑1 (MQ=255)                                                                                                                            gTGCAGCATGAAGAACGCATTCGGAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGCTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGc                  >  1:298181/1‑151 (MQ=255)                                                                                                                               cagcaTGAAGAACGCATTCGGAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCa                >  2:794052/1‑150 (MQ=255)                                                                                                                                      aagaaCGCATTCGGAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGc        <  1:29542/151‑1 (MQ=255)                                                                                                                                       agaaCGCATTCGGAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCt       >  2:79057/1‑151 (MQ=255)                                                                                                                                       agaaCGCATTCGGAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCt       >  1:27358/1‑151 (MQ=255)                                                                                                                                             cATTCGGAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCCGCTGAAACGCTATACCCCAAACGTGGCTGATGc  >  2:16928/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                CTGTTCGGCATTCCTGATCAGGAAGAGGAGACGAACAAATTTGCGATTCAGATCCCTTACGCACTGGGCATCATTGCAACGCGTTCCGTGGATACCCCGGTTATCGGCCTGAAAGAGCTGATGGTGCAGCATGAAGAACGCATTCGTAACGGGATGAAGGCGTACTCTCTGCTCGAACAACTGCGTTCTGGTTCTACCGACCAGGCGGTTCGTGACCAGTTCAATAGCATGAAGAAAGACCTCGGTTACGGTCTGCTGCTGAAACGCTATACGCCAAACGTGGCTGATGC  >  NC_000913/772268‑772557

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