Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,275,084	G→C	14.3%	intergenic (‑281/+31)	PP_5462 ← / ← PP_5463	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,275,084	0	G	C	14.3%	66.2 / 5.7	28	intergenic (‑281/+31)	PP_5462/PP_5463	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (14/10);  new base C (2/2);  total (16/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

CGAATGGGCATTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGATAAGTTGTAT  >  NC_002947/1274942‑1275220                                                                                                                                               |                                                                                                                                         cGAATGGGCATTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCGAAGCCCGCCCATcggc                                                                                                                                    >  7:3987/1‑149 (MQ=255) cGAATGGGCATTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATgggc                                                                                                                                    <  3:12052/149‑1 (MQ=255)     tGGGCATTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATgggcgggc                                                                                                                                <  1:51616/149‑1 (MQ=255)          aTTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCGAAGCCCGCCCATCGGCGGGCtttgc                                                                                                                           <  8:3987/149‑4 (MQ=255)           ttttGGTCTATGCGGCGGAGAGTTAATGTAGATGCCTGATTCCCGGTTACGTCGGCTCCAATAATGTGCGGGAGAGCGTTAATGATGGGGGTGAGGGTGGTAAAGGTGATAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAgg                                                                                                                          <  4:105915/149‑1 (MQ=255)                                            gCCTGCTTCCCGGTTCCGTCGGCTCCAAAAAGGGGCGGGCGAGCGGTAATGATGGGGGGGATGGTGGTAAAGGTGAACGGAGAATGACCAAACCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACg                                                                                               <  6:10699/143‑1 (MQ=255)                                            gCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACg                                                                                               >  5:10699/1‑143 (MQ=255)                                                             gTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAgg                                                                        >  2:83105/1‑149 (MQ=255)                                                                      aaTAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCa                                                               >  7:106001/1‑149 (MQ=255)                                                                      aaTAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCACAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAACGTAGGAAGCa                                                               >  7:96171/1‑149 (MQ=255)                                                                                          tAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGACCGCACTCCtcat                                             >  3:139947/1‑147 (MQ=255)                                                                                              gatgGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTAt                                                                                                        >  5:89096/1‑84 (MQ=255)                                                                                                  ggggTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCg                                    >  8:173260/1‑148 (MQ=255)                                                                                                   gggTGATGTGGCTAAAGGTCAAAGCAGACTGACCACCCCCCGCCGCTGTGCGGGCTTCAGCGGGGCGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCATGAATCGCa                                  <  8:106001/149‑1 (MQ=255)                                                                                                        aTGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGAAGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACttt                              >  3:196125/1‑148 (MQ=255)                                                                                                                aaaGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAg                           <  5:26832/143‑1 (MQ=255)                                                                                                                aaaGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAg                           >  6:26832/1‑143 (MQ=255)                                                                                                                           gagaATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAgg                          <  5:215152/133‑1 (MQ=255)                                                                                                                           gagaATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAgg                          >  6:215152/1‑133 (MQ=255)                                                                                                                            agaATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGAt           >  7:95263/1‑147 (MQ=255)                                                                                                                             gaATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCa                                                     >  5:17499/1‑104 (MQ=255)                                                                                                                             gaATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCa                                                     <  6:17499/104‑1 (MQ=255)                                                                                                                              aaTGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGATAAGt       <  4:117235/149‑1 (MQ=255)                                                                                                                              aaTGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGATAAGt       >  3:117235/1‑149 (MQ=255)                                                                                                                               atccgCAAAGCCCGCCGATGGGCCGGCCTCAGGGGGGCGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGATAAGtt      <  1:60296/145‑1 (MQ=255)                                                                                                                                 ccctgaaGCCCGCCCATCGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGATAAGtt      >  4:280438/6‑147 (MQ=255)                                                                                                                                 ccctgaaGCCCGCCCATCGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGATAAGtt      <  3:280438/142‑1 (MQ=255)                                                                                                                                   gCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATGAGGAATCGGACTTTAAGGTGACTGGCTTTGGATAAGTTGTAt  >  2:270178/1‑149 (MQ=255)                                                                                                                                           cgccgATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGATAAGtt      <  2:294575/137‑1 (MQ=255)                                                                                                                                           cgccgATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGATAAGtt      >  1:294575/1‑137 (MQ=255)                                                                                                                                               |                                                                                                                                         CGAATGGGCATTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGATAAGTTGTAT  >  NC_002947/1274942‑1275220

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: ‑