Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,645,945	A→T	15.5%	L32Q (CTG→CAG)	PP_3212 ←	Rieske 2Fe‑2S family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,645,945	0	A	T	15.5%	61.3 / 3.6	25	L32Q (CTG→CAG)	PP_3212	Rieske 2Fe‑2S family protein
Reads supporting (aligned to +/- strand):  ref base A (9/12);  new base T (2/2);  total (11/14)
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

AGTCGTATTCCCAGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCA  >  NC_002947/3645796‑3646080                                                                                                                                                      |                                                                                                                                        aGTCGTATTCCCAGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCca                                                                                                                                         >  2:11841/1‑150 (MQ=255)   tCGTATTCCCAGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCcagc                                                                                                                                       >  1:3933/1‑150 (MQ=255)           ccAGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCTGCTGGGTGaa                                                                                                                               >  1:166314/1‑150 (MQ=255)            cAGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAAt                                                                                                                              >  1:386107/1‑150 (MQ=255)               ccGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCTGCTGGGTGaagcgc                                                                                                                           <  2:166314/150‑5 (MQ=255)               ccGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGc                                                                                                                           >  2:132823/1‑150 (MQ=255)                     tACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTgggg                                                                                                                     <  1:11841/150‑1 (MQ=255)                        ccGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCTCGCATTATCACCTGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCg                                                                                                                  >  1:296998/1‑150 (MQ=255)                                      gCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGt                                                                                                    <  2:386107/150‑1 (MQ=255)                                        ccTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCa                                                                                                  >  2:628477/1‑150 (MQ=255)                                                 tGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCa                                                                                         >  2:732518/1‑150 (MQ=255)                                                    cGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAggg                                                                                      <  1:132823/150‑1 (MQ=255)                                                          cTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATAcc                                                                                <  2:3933/150‑1 (MQ=255)                                                               gACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGaaa                                                                           <  2:296998/150‑1 (MQ=255)                                                                   gTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGgc                                                                       <  1:732518/150‑1 (MQ=255)                                                                    tttGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGgcg                                                                      >  1:116046/1‑150 (MQ=255)                                                                             gcggtgGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCGTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAgg                                                             <  2:761491/150‑1 (MQ=255)                                                                                  ggCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGgttgtt                                                        <  2:736758/150‑1 (MQ=255)                                                                                   gcagcaGCGGTCCCGCATTGCCACCGGGGGACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGttt                                                       >  1:415947/1‑150 (MQ=255)                                                                                                     tGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGTCGCAATACAGGTTGTTTGGTAACCAACATAATCAg                                     <  2:347564/150‑1 (MQ=255)                                                                                                                  ccATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCg                                                                                                            >  2:480612/1‑66 (MQ=255)                                                                                                                  ccATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCg                                                                                                            <  1:480612/66‑1 (MQ=255)                                                                                                                   cATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCg                       >  2:547379/1‑150 (MQ=255)                                                                                                                                  acaacaCCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATgg        <  1:628477/150‑1 (MQ=255)                                                                                                                                        ccAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCa  <  2:116046/150‑1 (MQ=255)                                                                                                                                             cactCACCCTGCTGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTg                                                      >  2:599449/4‑93 (MQ=255)                                                                                                                                             cactCACCCTGCTGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTg                                                      <  1:599449/90‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                        AGTCGTATTCCCAGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCA  >  NC_002947/3645796‑3646080

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