Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP006881	2,181,354	A→T	11.8%	L124L (CTT→CTA)	U712_10920 ←	Menaquinol‑cytochrome c reductase cytochrome b/c subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	2,181,354	0	A	T	11.8%	84.4 / 4.8	34	L124L (CTT→CTA)	U712_10920	Menaquinol‑cytochrome c reductase cytochrome b/c subunit
Reads supporting (aligned to +/- strand):  ref base A (15/15);  new base T (2/2);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.81e-01

TGTTCTTCTGCCTTGGCCCAGTCGTGAGTAGCAACTGACTGCCACGTCAAAAATACAGCAGCGGAAATCGCCAAAAGCATCATACCAACAGCAACAGGACGCTTCCACGGCCTTCTCTCTGTTCCTCTGTCAAGGAATGGCGCCAAAAGCAGGGCGCCAAAAGCGAGTCCCGGCATAATCATGGCGCCGACCACTGTAAAGCTTCCGGCGGCGTATTCATACTTCAACAGCTGGTATAAGAAAAGAAAATACCAATCCGGCAGC  >  CP006881/2181223‑2181486                                                                                                                                    |                                                                                                                                     tGTTCTTCTGCCTTGGCCCAGTCGTGAGTAGCAACTGACTGCCACGTCAAAATTACAGCAGCGGAAATCGCCAAAAGCATCATACCAACAGCAACAGGACTCTTCCACGGCCTTCTCTCTGTTCCTCTGTCAAGGAATg                                                                                                                               <  6:174970/139‑1 (MQ=255)   ttcttcTGCCTTGGCCCAGTCGTGAGTAGCAACTGACTGCCAAGTCAAAAAAAAAGAAGCGGAAATCGCCAAAAGCAACATACCAACACCAACAGGACGCTACCACGGCCTACTCTCTCTTCCACTGACAATGAATGgc                                                                                                                             >  8:189051/1‑139 (MQ=255)               ggCCCAGTCGTGAGTAGCATCTGTCTGCCACGTCAAAAATACAGCAGCGGAAATCGCCAAAAGCATCTTACCAACAGCAACAGGACGCTTCCACGGCCTTCTCTCTGTTCCTCTGTCAAGGAATGGCGCCAAAAGCAg                                                                                                                  <  7:189051/138‑1 (MQ=255)               ggCCCAGTCGTGAGTAGCAACTGACTGCCACGTCAAAAATACAGCAGCGGAAATCGCCAAAAGCATCATACCAACAGCAACAGGACGCTTCCACGGCCTTCTCTCTGTTCCTCTGTCAAGGAATGGCGCCAAAAGCAg                                                                                                                  <  5:15009/138‑1 (MQ=255)               ggCCCAGTCGTGAGTAGCAACTGACTGCCACGTCAAAAATACAGCAGCGGAAATCGCCAAAAGCATCATACCAACAGCAACAGGACGCTTCCACGGCCTTCTCTCTGTTCCTCTGTCAAGGAATGGCGCCAAAAGCAg                                                                                                                  <  7:185329/138‑1 (MQ=255)                     gTCGTGAGTAGCAACTGACTGCCCCGTCAAAAATACAGCAGCGGAAATCGCCAAAAGCATCATACCAACACCAACAGGACGCTTCCACGGCCTTCTCGCTGTTGAACTGGAAAGGAAATGCGCCAAAACCACGGCGccc                                                                                                           >  6:15009/1‑138 (MQ=255)                     gTCGTGAGTAGCAACTGACTGCCACGTCAAAAATACAGCAGCGGAAATCGCCAAAAGCATCATACCAACAGCAACAGGACGCTTCCACGGCCTTCTCTCTGTTCCTCTGTCAAGGAATGGCGCCAAAAGCAGGGCGcc                                                                                                            >  5:71035/1‑138 (MQ=255)                                cAACTGACTGCCACGTCAAAAATACAGCAGCGGAAATCGCCAAAAGCATCATACCAACAGCAACAGGACGCTTCCACGGCCTTCTCTCTGTTCCTCTGTCAAGGAATGGCGCCAAAAGCAGGGCGCCAAAAGCGAGTcc                                                                                                >  2:235974/1‑139 (MQ=255)                                cAACTGACTGCCACGTCAAAAATACAGCAGCGGAAATCGCCAAAAGCATCATACCAACAGCAACAGGACGCTTCCACGGCCTTCTCTCTGTTCCTCTGTCAAGGAATGGCGCCAAAAGCAGGGCGCCAAAAGCGAGTcc                                                                                                <  6:71035/139‑1 (MQ=255)                                                cAAAAATACAGCAGCGGAAATCGCCAAAAGCATCATACCAACAGCAACAGGACGCTTCCACGGCCTTCTCTCTGTTCCTCTGTCAAGGAATGGCGCCAAAAGCAGGGCGCCAAAAGCGAGTCCCGGCATAATCATGgcg                                                                                <  1:235974/139‑1 (MQ=255)                                                           cagcGGAAATCGCCAAAAGCATCATACCAACAGCAACAGGACGCTTCCACGGCCTTCTCTCTGTTCCTCTGTCAAGGAATGGCGCCAAAAGCAGGGCGCCAACAGCGAGTCCCGGCCTAATCATGGCGCCGCCCCCTg                                                                      >  4:256569/1‑138 (MQ=255)                                                              cGGAAATCGCCAAAAGCATCATACCAACAGCAACAGGACGCTTCCACGGCCTTCGCTCGGTTCGTCAGGCAAGGAATGGCGCCAAAAGAAGGGCGCCAAAAGCGACTACCGGCATAATCATGGCGCCGACCACTGTa                                                                    >  4:78982/1‑137 (MQ=255)                                                                        cAAAAGCATCATACCAACAGCAACAGGACGCTTCCACGGCCTTCTCTCTGTTCCTCTGTCAAGGAATGGCGCCAAAAGCAGGGCGCCAAAAGCGAGTCCCGGCATAATCATGGCGCCGCCCACTGTAAAGCTTCcggcg                                                        >  2:195629/1‑139 (MQ=255)                                                                           aaGCATCATACCGACAGCGACCGGACGCTTCCACGGCCTTCTTTCTGTCCCTCTGTCTAGGAATGGAGCCAAAAGCAGGGCCCCAAAAGCGAGACCTGGCATAATCATGGCGCCGAc                                                                           <  2:41862/117‑1 (MQ=255)                                                                           aaGCATCATACCGACAGCGACCGGACGCTTCCACGGCCTTCTTTCTGTCCCTCTGTCTAGGAATGGAGCCAAAAGCAGGGCCCCAAAAGCGAGACCTGGCATAATCATGGCGCCGAc                                                                           >  1:41862/1‑117 (MQ=255)                                                                               atcatACCAACAGCAACAGGACGCTTCCACGGCCTTCTCTCTGTTCCTCTGTCAAGGAATGGCGCCAAAACCAGGGCGCCAAAAGCGAGTCCCGGCATAATCATGGCGCCGACCACTGTAAAGCTTCCGGCGGCGTAtt                                                 <  5:89452/139‑1 (MQ=255)                                                                                 catACCAACAGCAACAGGACGCTTCCACGGCCTTCTCTCTGTTCCTCTGTCAAGGAATGGCGCCAAAAGCAGGGCGCCAAAAGCGAGTCCCGGCATAATCATGGCGCCGACCACTGTAAAGCTTCCGGCGGCGTATTCa                                               <  1:195629/139‑1 (MQ=255)                                                                                                    cGCTTCCACGGCCTTCTCTCTGTTCCTCTGTCAAGGAATGGCGCCAAAAGCAGGGCGCCAAAAGCGAGTCCCGGCATAATCATGGCGCCGACCACTGTAAAGCTTCCGGCGGCGTATTCATACTTCAACAGCTGGtat                             >  5:110516/1‑138 (MQ=255)                                                                                                       ttCCACGGCCTTCTCTCTGTTCCTCTGTCAAGGAATGGCGCCAAAAGCAGGGCGCCAAAAGCGAGTCCCGGCATAATCATGGCGCCGACCCCTGTAAAGCTTCCGGTGGCGTATTCATACTTGAACAGCTGGttttag                          >  8:164222/1‑133 (MQ=255)                                                                                                         ccACGGCCTTCTCTCTGTTCCTCTGTCAAGGAATGGCGCCAAAAGCAGGGCGCCAAAAGCGAGTCCCGGCATAATCATGGCGCCGACCACTGTAAAGCTTCCGGCGGCGTATTCATACTTCAACAGCTGGTATaagaa                        <  5:81071/138‑1 (MQ=255)                                                                                                               ccTTCTCTCTGTTCCTCTGTCAAGGAATGGCGCCAAAAGCAGGGCGCCAAAAGCGAGTCCCGGCATAATCATGGCGCCGAc                                                                           >  7:87914/1‑81 (MQ=255)                                                                                                               ccTTCTCTCTGTTCCTCTGTCAAGGAATGGCGCCAAAAGCAGGGCGCCAAAAGCGAGTCCCGGCATAATCATGGCGCCGAc                                                                           <  8:87914/81‑1 (MQ=255)                                                                                                                 ttCACTCTGTTCCTCTGTCAAGGAATGGCGCCAAAAGCAGGGCGCCAAAAGCGAGTCCCGGCATAATCATGGCGCCGACCACTGTAAAGCTTCCCGCGGCGTATTCATAATTAAACACCGTGTTTAAgaaaagaaaaa                >  6:81071/1‑137 (MQ=255)                                                                                                                   ctctctGTTCCTATGTCAAGGAATGGCGCCAAAAGCAGGGCGCCAAAAGCGAGTCCCGGCATAATCATGGCGCCGACCACTGTAAAGCATCCGGCGGCGTATTCATACATCAACACCTGGTATAAGAAAAGAAAAAAc              >  4:176204/1‑138 (MQ=255)                                                                                                                         gtcccTCTGTCTAGGAATGGAGCCAAAAGCAGGGCCCCAAAAGCGAGACCTGGCATAATCATGGCGCCGACCACTGTAAAGCTTCCGGCGGCGTATTCATACTTCAACAGCTGGTATAAGAAAAGAAAATAc              <  4:286317/129‑1 (MQ=255)                                                                                                                         gtcccTCTGTCTAGGAATGGAGCCAAAAGCAGGGCCCCAAAAGCGAGACCTGGCATAATCATGGCGCCGACCACTGTAAAGCTTCCGGCGGCGTATTCATACTTCAACAGCTGGTATAAGAAAAGAAAATAc              >  3:286317/4‑132 (MQ=255)                                                                                                                         gTTCCTCTGTCAAGGAATGGCGCCAAAAGCAGGGCGCCAAAAGCGAGTCCCGGCATAATCATGGCGCCGACCACTGTa                                                                    >  5:97654/1‑78 (MQ=255)                                                                                                                         gTTCCTCTGTCAAGGAATGGCGCCAAAAGCAGGGCGCCAAAAGCGAGTCCCGGCATAATCATGGCGCCGACCACTGTa                                                                    <  6:97654/78‑1 (MQ=255)                                                                                                                         gTTCCTCTGTCAAGGAATGGCGCCAAAAGCAGGGCGCCAAAAGCGAGTCCCGGCATAATCATGGCGCCGACCACTGTAAAGCTTCCGGCGGCGTATTCATACTTCAACAGCTGGTATAAGAAAAGAAAATACCAATcct       >  5:147294/1‑138 (MQ=255)                                                                                                                          ttCCTATGTAAAGGAATGGCGCCAAAAGCTGGGCGCCAAAAGCGAGTCCCGGCATAATCATGGCGCCGACCACTGTAAAGCTTCCGGCGGCGTATTCATACTTCAACAGCTGGTATAAGAAAAGAAAATACCAATCCgg      <  1:112293/139‑1 (MQ=255)                                                                                                                             ctctGTCAAGGAATGGCGCCAAAAGCAGGGCGCCAAAAGCGAGTCCCGGCATAATCATGGCGCCGACCACTGTAAAGCTTCCGGCGGCGTATTCATACTTCAACa                                     >  7:269033/1‑105 (MQ=255)                                                                                                                             ctctGTCAAGGAATGGCGCCAAAAGCAGGGCGCCAAAAGCGAGTCCCGGCATAATCATGGCGCCGACCACTGTAAAGCTTCCGGCGGCGTATTCATACTTCAACa                                     <  8:269033/105‑1 (MQ=255)                                                                                                                                tGTCAAGGAATGGCGCCAAAAGCAGGGCGCCAAAAGCGAGTCCCGGCATAATCATGGCGCCGACCACTGTAAAGCTTCCGGTGGCGTATTCATACTTCAACAGCTGGTATAAGAAAAGAAAATACCAATCCGGCAGc  <  7:164222/137‑1 (MQ=255)                                                                                                                                tGTCAAGGAATGGCGCCAAAAGCAGGGCGCCAAAAGCGAGTCCCGGCATAATCATGGCGCCGACCACTGTAAAGCTTCCGGCGGCGTATTCATACTTCAACAGCTGGTATAAGAAAAGAAAATACCAATCCGGCAGc  <  5:216154/137‑1 (MQ=255)                                                                                                                                    |                                                                                                                                     TGTTCTTCTGCCTTGGCCCAGTCGTGAGTAGCAACTGACTGCCACGTCAAAAATACAGCAGCGGAAATCGCCAAAAGCATCATACCAACAGCAACAGGACGCTTCCACGGCCTTCTCTCTGTTCCTCTGTCAAGGAATGGCGCCAAAAGCAGGGCGCCAAAAGCGAGTCCCGGCATAATCATGGCGCCGACCACTGTAAAGCTTCCGGCGGCGTATTCATACTTCAACAGCTGGTATAAGAAAAGAAAATACCAATCCGGCAGC  >  CP006881/2181223‑2181486

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