Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,505,544	A→C	11.5%	D219A (GAC→GCC)	petB →	cytochrome b

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,505,544	0	A	C	11.5%	112.9 / 3.9	41	D219A (GAC→GCC)	petB	cytochrome b
Reads supporting (aligned to +/- strand):  ref base A (21/15);  new base C (2/3);  total (23/18)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGTGGTGACTACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGTTCTTCTTCCCTGAAATGGGCGGT  >  NC_002947/1505397‑1505689                                                                                                                                                    |                                                                                                                                                  ttCGTGGTGACTACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGAc                                                                                                                                                  <  1:164166/149‑1 (MQ=255)        tgACTACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGTCGACGCCATCaggg                                                                                                                                           <  6:312258/149‑4 (MQ=255)               cTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGTCGACGCCATCag                                                                                                                                             >  1:216261/1‑139 (MQ=255)               cTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGTCGACGCCATCag                                                                                                                                             <  2:216261/140‑2 (MQ=255)                  aTCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGTCGACGCCATCagg                                                                                                                                            >  7:53933/1‑136 (MQ=255)                  aTCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGTCGACGCCATCagg                                                                                                                                            <  8:53933/138‑3 (MQ=255)                                  tGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCg                                                                                                                 >  1:38940/1‑149 (MQ=255)                                    aaCCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCt                                                                                                               >  4:217191/1‑149 (MQ=255)                                      ccGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTg                                                                                                              >  4:172924/1‑148 (MQ=255)                                         cttcttCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATACTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCATGATGGCGTCGACATCAAGAAGAAAAAGGAc                                                                                                                                >  3:146529/1‑127 (MQ=255)                                            cttcGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTa                                                                                                       >  2:180761/1‑149 (MQ=255)                                                gCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCaag                                                                                                                                            >  8:209597/1‑108 (MQ=255)                                                gCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCaag                                                                                                                                            <  7:209597/108‑1 (MQ=255)                                                 cccTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCg                                                                                                  >  8:52824/1‑149 (MQ=255)                                                      cATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTcc                                                                                              <  1:180761/148‑1 (MQ=255)                                                                cccTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGtt                                                                                                >  7:324762/1‑136 (MQ=255)                                                                cccTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGtt                                                                                                <  8:324762/136‑1 (MQ=255)                                                                   tGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACAcc                                                                                >  1:317984/1‑149 (MQ=255)                                                                      cAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTc                                                                                                                   <  3:290761/111‑1 (MQ=255)                                                                      cAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTc                                                                                                                   >  4:290761/1‑111 (MQ=255)                                                                      cAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGtt                                                                             <  3:217191/149‑1 (MQ=255)                                                                                  tcGGCCTGCTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCATGATGACGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTc                                                                 >  6:349992/1‑149 (MQ=255)                                                                                     gCCTGGTTGTATTCCACATCCTCGCCTTGCACGAAGTGGGATCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGc                                                              <  8:274148/149‑1 (MQ=255)                                                                                                              cTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTCCCTGTTCGTGTTCt                                     >  1:321426/1‑149 (MQ=255)                                                                                                                 gCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTg                                    <  2:38940/147‑1 (MQ=255)                                                                                                                   aCGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGc                                                                                                                >  3:330630/1‑69 (MQ=255)                                                                                                                   aCGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGc                                                                                                                <  4:330630/69‑1 (MQ=255)                                                                                                                         tGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCgtggtg                          >  3:239394/1‑149 (MQ=255)                                                                                                                           ggTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGATATTGTCGGCGTGGTGGTGTTCCTGt                                              <  2:63477/127‑1 (MQ=255)                                                                                                                           ggTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGATATTGTCGGCGTGGTGGTGTTCCTGt                                              >  1:63477/1‑127 (MQ=255)                                                                                                                             ttCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGACGGGGTGttct                      >  8:333848/1‑149 (MQ=255)                                                                                                                                 aacaacCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGttcttctt                  >  4:264940/1‑149 (MQ=255)                                                                                                                                       ccTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGTTCTTCTTCCCTg             >  3:293146/1‑148 (MQ=255)                                                                                                                                         tGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCgtg                             <  4:245046/130‑1 (MQ=255)                                                                                                                                         tGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCgtg                             >  3:245046/1‑130 (MQ=255)                                                                                                                                             ggCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCgtggtgg                                                       <  8:103642/100‑1 (MQ=255)                                                                                                                                             ggCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCgtggtgg                                                       >  7:103642/1‑100 (MQ=255)                                                                                                                                             ggCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGTTCTTCTTCCCTGAAATgg       <  4:293146/148‑1 (MQ=255)                                                                                                                                                gTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGGTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATGGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGttg                       >  7:54511/1‑128 (MQ=255)                                                                                                                                                ctCGACATCAAGAAGAAAAAGCACGAAAACGGTATTCCGCTGTATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGttc                       <  8:54511/128‑1 (MQ=255)                                                                                                                                                  cGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTCCCTGTTCGTGTTCTGCGCCGTGGTGTTCTTCTTCCCTGAAATGGGCGGt  <  2:321426/148‑1 (MQ=255)                                                                                                                                                    aCATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGTTCTTCTTCCCTGAAATGGGc     <  8:183267/143‑1 (MQ=255)                                                                                                                                                    aCATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGTTCTTCTTCCCTGAAATGGGc     >  7:183267/1‑143 (MQ=255)                                                                                                                                                    |                                                                                                                                                  TTCGTGGTGACTACCTGATCTCGGGCATCACCCTGAACCGCTTCTTCGCCCTGCATGTGGTTGCCCTGCCAATCGTGATTCTCGGCCTGGTGGTATTGCACATCCTCGCCTTGCACGAAGTGGGTTCGAACAACCCTGATGGCGTCGACATCAAGAAGAAAAAGGACGAAAACGGTATTCCGCTGGATGGTATTCCGTTCCACCCGTACTACACCGTTAAAGACATTGTCGGCGTGGTGGTGTTCCTGTTCGTGTTCTGCGCCGTGGTGTTCTTCTTCCCTGAAATGGGCGGT  >  NC_002947/1505397‑1505689

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