Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,491,397	A→G	T166A (ACA→GCA)	cybB →	cytochrome b561

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,491,397	0	A	G	97.8%	150.8 / NA	45	T166A (ACA→GCA)	cybB	cytochrome b561
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (18/26);  minor base C (0/1);  total (18/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCTATTATATTTATTTGTATGAT  >  NC_000913/1491267‑1491531                                                                                                                                   |                                                                                                                                       ttCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTTGCACACCCCTATTTCTGGAAGAACAACGCACTTCTAc                                                                                                                               >  1:306271/1‑140 (MQ=255) ttCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACGCACTTCTAc                                                                                                                               >  2:460134/1‑140 (MQ=255)   cAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACGCACTTCTACGc                                                                                                                             >  1:73690/1‑140 (MQ=255)    agagGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACGCACTTCTACGCa                                                                                                                            <  2:177356/140‑1 (MQ=255)        gCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACGCACTTCTACGCatgat                                                                                                                        <  2:779204/140‑1 (MQ=255)             tttCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACGCACTTCTACGCATGATGCcgc                                                                                                                   >  2:666574/1‑140 (MQ=255)             tttCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGGTGCACGCTGCGGCTGCACTGGCACACCACCATTACAGCCAGGACAACGCACAACTACGCAGGATGCcgc                                                                                                                   >  2:397464/1‑140 (MQ=255)                  aaCGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACGCACTTCTACGCATGATGCCGCGTaaa                                                                                                              <  2:135531/140‑1 (MQ=255)                       gTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACGCACTTCTACGCATGATGCCGCGTAAACGTTc                                                                                                         <  1:731528/140‑1 (MQ=255)                              gCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACGCACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAgg                                                                                                  <  1:193630/140‑1 (MQ=255)                              gCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACGCACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAgg                                                                                                  <  1:759977/140‑1 (MQ=255)                                 tAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACGCACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGata                                                                                               >  2:459833/1‑140 (MQ=255)                                    aGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACGCACTTCGACGCATGTTGCCGCGTAAACGCTCCTGAAGGATAttt                                                                                            >  1:352096/1‑140 (MQ=255)                                        gTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACGCACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAg                                                                                        >  1:75519/1‑140 (MQ=255)                                         tGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACGCACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGa                                                                                       >  2:411193/1‑140 (MQ=255)                                         tGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACGCACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGa                                                                                       >  2:16360/1‑140 (MQ=255)                                            cATGAAACGCTGGCGAATCTGGAATCTTTGGTCATCGTGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACGCACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGaaaa                                                                                    <  1:368872/140‑1 (MQ=255)                                            cATGAAACGCTGGCGAAACTGGTATATTTTGTCATTGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACGCACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGaaaa                                                                                    <  2:306271/140‑1 (MQ=255)                                               gAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACGCACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGc                                                                                 <  1:459833/140‑1 (MQ=255)                                                aaaCGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACGCACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGcc                                                                                <  2:73690/140‑1 (MQ=255)                                                aaaCGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACGCACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGcc                                                                                <  1:16360/140‑1 (MQ=255)                                                     cTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACGCACTTCTACGCATGATGCCGCGTAAACGTTCGTGACGGATATTTAGAGAAAACGCCTGTAc                                                                           >  1:362297/1‑140 (MQ=255)                                                        gCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACGCACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTaa                                                                        >  2:635937/1‑140 (MQ=255)                                                            aTCTGGGATATTTTGTTTTCGGCTTGCACGCTGCGCCTGCACTGGCACGCCACTATTTCTGGAAGGACAACGCACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAAcc                                                                    <  2:352096/140‑1 (MQ=255)                                                             tCTGGGATATTTTGTTATCGGGGTGCACGTTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACGCACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCg                                                                   <  1:397464/140‑1 (MQ=255)                                                             tCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACGCACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCg                                                                   <  1:411193/140‑1 (MQ=255)                                                             tCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACGCACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCg                                                                   <  2:538083/140‑1 (MQ=255)                                                                gggATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACGCACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGAcc                                                                <  2:615157/140‑1 (MQ=255)                                                                gggATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACGCACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGAcc                                                                <  1:635937/140‑1 (MQ=255)                                                                       tttGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACGCACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTa                                                         <  1:419013/140‑1 (MQ=255)                                                                       tttGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACGCACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTa                                                         <  2:189621/140‑1 (MQ=255)                                                                           tCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACGCACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAg                                                      >  2:482857/1‑139 (MQ=255)                                                                           tCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACGCACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAg                                                      <  1:482857/139‑1 (MQ=255)                                                                           tCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACGCACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTGAAAGAAAGCGCCTGTACTAAAGCCGACCCGTGGTACAgg                                                     >  2:244897/1‑140 (MQ=255)                                                                             aTCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACGCACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCg                                                   <  1:666574/140‑1 (MQ=255)                                                                               cGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACGCACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCgaa                                                 <  1:84281/140‑1 (MQ=255)                                                                                    tGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACGCACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATa                                            <  1:738732/140‑1 (MQ=255)                                                                                    tGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACGCACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATa                                            <  1:460134/140‑1 (MQ=255)                                                                                     gCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACGCACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATAc                                           >  2:526758/1‑140 (MQ=255)                                                                                           tGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACGCACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCt                                     >  1:418023/1‑140 (MQ=255)                                                                                                      tGGCACACCACTATTTCTGGAAGGACAACGCACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAgc                          <  2:75519/140‑1 (MQ=255)                                                                                                          acacCACTATTTCTGGAAGGACAACGCACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCt                      >  1:621563/1‑140 (MQ=255)                                                                                                           caccacTATTTCTGGAAGGACAACGCACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCta                     <  1:693015/140‑1 (MQ=255)                                                                                                             ccacTATTTCTGGAAGGGCAACCCACTTCTACGAACGATTCCGCGTAATCGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCtatt                   <  2:362297/140‑1 (MQ=255)                                                                                                                        ggAAGGACAACGCACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCTATTATATTTATTTg        >  2:694539/1‑140 (MQ=255)                                                                                                                              acaacgCACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCTATTATATTTATTTGTATGat  <  2:621563/140‑1 (MQ=255)                                                                                                                                   |                                                                                                                                       TTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCTATTATATTTATTTGTATGAT  >  NC_000913/1491267‑1491531

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