Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,275,084	G→C	21.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	21.3%	84.5 / 14.8	46	intergenic (‑281/+31)	PP_5462/PP_5463	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (21/15);  new base C (3/7);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 1.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTGCGAATGGGCATTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCA‑GGGGGGGGCAGAAGGGCTAT‑CGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGATAAGTTGTATCCAGC  >  NC_002947/1274938‑1275225                                                                                                                                                   |                                                                                                                                                aGTGCGAATGGGCATTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGAt                                                                                                                                               <  2:77254/149‑1 (MQ=255)  gTGCGAATGGGCATTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATg                                                                                                                                              <  6:86277/149‑1 (MQ=255)     cGAATGGGCATTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATgggc                                                                                                                                           <  3:89926/149‑1 (MQ=255)          gggCATTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCt                                                                                                                                      >  3:335055/1‑149 (MQ=255)           ggCATTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCCATCGGCGGGCtt                                                                                                                                     <  8:254627/149‑1 (MQ=255)            gaaTTTTGGTCTATACTGCGGAGAGTTAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCTGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCCATCGGCGGGCttt                                                                                                                                    <  3:243539/147‑2 (MQ=255)                  ttgTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCA‑gggg                                                                                                                              <  4:98399/146‑1 (MQ=255)                  tGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCA‑gggg                                                                                                                              <  5:81849/148‑1 (MQ=255)                     ccTATCCGGCGGAGAGGGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCA‑gggggggg                                                                                                                          <  2:312788/148‑1 (MQ=255)                          ccGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCA‑GGGGGGGGCAg                                                                                                                       >  1:157749/1‑147 (MQ=255)                           cggcggAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGGGGGAGAGCGGGAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGTGCAAAGCCCCCCCATTGGCGGGCTTCA‑GGGGGGGGCAGAAg                                                                                                                    >  6:297260/1‑149 (MQ=255)                           cggcggAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTTGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCA‑GGGGGGGGCAGAAg                                                                                                                    <  3:257222/149‑1 (MQ=255)                             gcggAGAGGGAAGGTCGAAGCCTGCCTCCCGGTTACGTCGGCACCAATAAGGTGCGGGAGCGCGGTTCTGCTGGGGGCTAGGGTGGCACAGGTTAAAGGAGACTGAGCAACGCCCGCCGATGGGCGGGCTTCA‑GGGGGGGGCAGAAggg                                                                                                                  <  2:293663/149‑1 (MQ=255)                                      gAATGTCGAAGCCTGCTTCCCGGTTACGTCGGGTCCAATAAGGTGCGGGGGAGCGGTTATGGTGGGGGTGGTGGGGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCGTCAGGGGGGGGGCGGAGGGGCGAT‑CGta                                                                                                         >  2:220621/1‑149 (MQ=255)                                             gAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCA‑ggggg                                                                                                                             <  2:16049/122‑1 (MQ=255)                                             gAAGCCTGCTTCCCGGTTACGTCGGCTACAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCA‑ggggg                                                                                                                             >  1:16049/1‑122 (MQ=255)                                                gCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCA‑GGGGGGGGCAGAAGGGCTAT‑CGTACTACGGTCGCt                                                                                              >  1:177337/1‑149 (MQ=255)                                                 ccTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCA‑GGGGGGGGCAGCAGGGCTAT‑CGTACTACGGTCGCt                                                                                              >  6:327318/1‑148 (MQ=255)                                                  cTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCA‑GGGGGGGGCAGAAGGGCTAT‑CGTACTACGGTCGCTtc                                                                                            <  4:335055/149‑1 (MQ=255)                                                     cTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCA‑GGGGGGGGCAGAAGGGCTAT‑CGTACTACGGTCGCTTCTCg                                                                                         >  5:147908/1‑149 (MQ=255)                                                        cccGGTTACGTCGGCTCCAATTAGGTGGGGGAGAGCGGGAATGATGGGGGTGATGGTGGTTAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCA‑GGGGGGGGCAGAAGGGGTTA‑CCTACTACGGTCGCTTCTGGTc                                                                                       >  2:146917/1‑148 (MQ=255)                                                               ccGTCGGCTCCCCTAAGGTGCGGGAGCGCGGTATTCACGGGGGCGCTGGTGCGAACGGTGAAAGGAGAATGAGCACACCCCCCCCATGGGCGGGCTTCA‑GGGGGGCGCAGAAGGGCTAT‑CGTACTACGGTCGCTTCTCGTCCACCGaaa                                                                               <  8:298707/148‑1 (MQ=255)                                                                 gTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCCATcggcgg                                                                                                                                         >  7:153035/1‑91 (MQ=255)                                                                 gTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCCATcggcgg                                                                                                                                         <  8:153035/91‑1 (MQ=255)                                                                  tCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCA‑GGGGGGGGCAGAAGGGCTAT‑CGTACTACGGTCGCTTCTCGTCCACCGAAAgg                                                                             >  6:139918/1‑148 (MQ=255)                                                                  tCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTACA‑GGGGGGGGCAGAAGGGCTAT‑CGTACTACGGTCGCGTCTAGTCCAACGAAAgg                                                                             >  7:330525/1‑148 (MQ=255)                                                                       tCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCA‑GGGGGGGGCAGAAGGGCTAT‑CGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAgg                                                                         >  7:20954/1‑147 (MQ=255)                                                                                gTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCA‑GGGGGGGGCAGAAGGGCTAT‑CGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCAcccgcc                                                              >  7:17292/1‑149 (MQ=255)                                                                                     ggAGCGCGGTACGGATGGGGGGCATGGTGGTATAGGTGACAGGAGACTGAGCAATGCCCGGCGATGGGCGGGCTTCA‑GGGGGGGGCCGAAGGGCTAT‑CGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCAc                                                         <  4:70223/149‑1 (MQ=255)                                                                                     ggAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCA‑GGGGGGGGCAGAAGGGCTAT‑CGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCAc                                                         >  4:162121/1‑149 (MQ=255)                                                                                             gTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCA‑GGGGGGGGCAGAAGGGCTAT‑CGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCtcat                                                  >  7:279034/1‑148 (MQ=255)                                                                                             gTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCA‑GGGGGGGGCAGAAGGGCTAT‑CGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCtcat                                                  >  1:47801/1‑148 (MQ=255)                                                                                              tAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCA‑GGGGGGGGCAGAAGGGCTAT‑CGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCtcatc                                                 >  3:33108/1‑148 (MQ=255)                                                                                                   atgGGGGTCATGGAGCTAAAGGGCAAAGGCGAACCCGCACAGCCCCCCGACCGGCGGGCTTCC‑GGGCCGCGCAGAAGGGCTAT‑CGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAAt                                           <  5:273281/149‑1 (MQ=255)                                                                                                     gggggTGCAGGTGGGAAAGGCGACCGGACCACGAGCAACCCCCGCCGCTGGGCGGGCTTCA‑GGGGGCGGCAGAAGGGCTAT‑CGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCg                                         <  8:17292/149‑1 (MQ=255)                                                                                                        ggTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCA‑GGGGGGGGCAGAAGGGCTAT‑CGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCAc                                      >  2:167924/1‑149 (MQ=255)                                                                                                        ggTGATGGTGGTAAAGGCGAACGGAGAATGAGCAACGCCCGCCGAGGGGCGGGCTTCA‑GGGGGGGGCAGAAGGGCTAT‑CGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCAc                                      <  6:147908/149‑1 (MQ=255)                                                                                                                    aaaGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCA‑GGGGGGGGCAGAAGGGCGGTACGTAATACGGTCGCTTCTCGTCCACCGAAAGGAAGGAACCACCCGCCCGCACTCCTCATAAGGAAACGCACATTCAGGTGa                            >  1:113990/1‑148 (MQ=255)                                                                                                                          gacccgagaATGCGCACAGCCCGCCCATGGGCGGGCGTCA‑GGGGGGGGCAGAAGGGCTAT‑CGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCttt                    <  6:173850/144‑1 (MQ=255)                                                                                                                              ggAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCA‑GGGGGGGGCAGAAGGGCTAT‑CGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCAACAGGAATCGCACTTTCATGTGACTGCCGTAGGa                 >  5:27809/1‑148 (MQ=255)                                                                                                                                 gaATGAGCAAAGCCCGCCGATGGGCGGGCTTCA‑GGGGGGGGCAGAAGGGCTAT‑CGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGATAAg             >  4:203167/1‑149 (MQ=255)                                                                                                                                 caatgcccAAAGCCCGCCCAGCGGCGGGCTTCA‑GGGGGCGGCAGAAGGGCTAT‑CGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGATAAg             <  1:146917/142‑1 (MQ=255)                                                                                                                                   ccccctgaaGCCCGCCCATCGGCGGGCTTTA‑GGGGGGGGCAGAAGGGCTAT‑CGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGATAAGtt           >  3:2528/8‑149 (MQ=255)                                                                                                                                   aggAGCCAACCCCGCCCCTGGGCGGGAGTCA‑GCGGTGGTCAGAAGGGCTAT‑CGTACTAAGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCTCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGATAAGtt           <  8:330525/147‑1 (MQ=255)                                                                                                                                   aTGACCAAAGCCCGCCGACCGGCGGCCATAA‑GGGGCGCGCAGAAGGGCTAT‑CGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGATAAGtt           <  8:279034/149‑1 (MQ=255)                                                                                                                                            gCCCGCCGATGGGCGGGCTTCA‑GGGGGGGGCAGAAGGGCTAT‑CGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGATAAGTTGTATCCAGc  >  2:123021/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                AGTGCGAATGGGCATTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCA‑GGGGGGGGCAGAAGGGCTAT‑CGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGATAAGTTGTATCCAGC  >  NC_002947/1274938‑1275225

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