Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,243,725	G→A	10.8%	intergenic (+14/+85)	tsaA → / ← rnt	peroxiredoxin/ribonuclease T

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,243,725	0	G	A	10.8%	120.9 / 4.9	46	intergenic (+14/+85)	tsaA/rnt	peroxiredoxin/ribonuclease T
Reads supporting (aligned to +/- strand):  ref base G (23/18);  new base A (3/2);  total (26/20)
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.63e-01

CTCTGCAGTTCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAA  >  NC_002947/1243587‑1243869                                                                                                                                           |                                                                                                                                                   ctctGCAGTTCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCgg                                                                                                                                          >  4:17802/1‑149 (MQ=255)        gTTCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATAAAAAAACCGGccacg                                                                                                                                     >  3:256357/1‑144 (MQ=255)            cACGAAGAGCACGGCTAAGTGTGCCCTGCCTACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGc                                                                                                                                <  3:98507/148‑1 (MQ=255)              cGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCgg                                                                                                                             >  8:304054/1‑149 (MQ=255)                      aCGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATAAAAAAACCGGCCACGAAGGCCGGttttttca                                                                                                                     <  4:256357/149‑3 (MQ=39)                         gCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGc                                                                                                                  >  8:197194/1‑149 (MQ=255)                                         cAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGcc                                                                                                 >  3:298294/1‑148 (MQ=255)                                               gAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGcc                                                                                                     <  8:35924/138‑1 (MQ=255)                                                 aaaaaGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTAt                                                                                                                    >  5:264822/1‑123 (MQ=255)                                                 aaaaaGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTAt                                                                                                                    <  6:264822/123‑1 (MQ=255)                                                 aaaaaGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCgg                                                                                        <  3:17802/149‑1 (MQ=255)                                                         gACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGccc                                                                                >  4:285251/1‑149 (MQ=255)                                                               ggCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAAt                                                                          >  1:158431/1‑149 (MQ=255)                                                                gCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTATTCGCGGGCGCGCCCGCGaa                                                                           >  6:255544/1‑147 (MQ=255)                                                                 cATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTATTCGCGGGCGCGCCCGCGaagag                                                                        <  5:255544/149‑4 (MQ=255)                                                                   tGAACGGTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTTTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCt                                                                      <  7:299941/149‑1 (MQ=255)                                                                    gAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCtt                                                                     >  2:50103/1‑149 (MQ=255)                                                                     aaCGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCtt                                                                     <  7:304054/148‑1 (MQ=255)                                                                      aCGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCtt                                                                     >  8:62287/1‑147 (MQ=255)                                                                       cGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAg                                                                  <  4:298294/149‑1 (MQ=255)                                                                        gCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAg                                                                  >  1:276112/1‑148 (MQ=255)                                                                          ttCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTTCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAg                                                                  >  5:246538/1‑146 (MQ=255)                                                                                tGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCa                                                                >  6:219384/1‑142 (MQ=255)                                                                                tGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCa                                                                <  5:219384/142‑1 (MQ=255)                                                                                  aaGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATc                                                        <  3:285251/148‑1 (MQ=255)                                                                                  aaGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCg                                                       <  7:197194/149‑1 (MQ=255)                                                                                         tgctgcCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACGGGCCTTCGTGGCCGTTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGaaa                                                 <  1:50103/148‑1 (MQ=255)                                                                                                   cTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAgcc                                      >  8:129113/1‑149 (MQ=255)                                                                                                    tGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGACCGCGAATGCCTTTAGCAGATCAATCGTTGAAATAGCGCCAgccg                                     >  7:285056/1‑149 (MQ=255)                                                                                                     gAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAgccgc                                    >  2:297274/1‑149 (MQ=255)                                                                                                        cGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATAAAAAAACCGGCCACGAAGGCCGGTTTTTTCAt                                                                                                                    <  2:60717/68‑1 (MQ=17)                                                                                                        cGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATAAAAAAACCGGCCACGAAGGCCGGTTTTTTCAt                                                                                                                    >  1:60717/1‑68 (MQ=17)                                                                                                         gagaACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAAt                                                         <  8:253533/124‑1 (MQ=255)                                                                                                         gagaACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAAt                                                         >  7:253533/1‑124 (MQ=255)                                                                                                         gagaACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCATAAGCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAgccgcc                                   >  1:108291/1‑148 (MQ=255)                                                                                                              cGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttctt                           >  8:75362/1‑149 (MQ=255)                                                                                                                cTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttctttc                         <  2:276112/149‑1 (MQ=255)                                                                                                                           tAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGa               <  7:62287/148‑1 (MQ=255)                                                                                                                            aaTTGCTGCGTGCATGAAAAAAGCGGCCTTCGTGGCGGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGGAGCTCAATGGTTGAAATCGCGGCCGCCGCCCATTTGTTTCCAGCGGTTGAc              >  7:53452/1‑148 (MQ=255)                                                                                                                              ttGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCGGATCAATCGtt                                                     >  1:27697/1‑107 (MQ=255)                                                                                                                              ttGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGtt                                                     <  2:27697/107‑1 (MQ=255)                                                                                                                                   gcGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCg       >  2:215914/1‑148 (MQ=255)                                                                                                                                    cGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGc      <  2:158431/148‑1 (MQ=255)                                                                                                                                    cGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCa     <  5:56570/149‑1 (MQ=255)                                                                                                                                    cGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCCAATGCCATTAGCAGATAAATAGTTTAAATACCGCCAGCCGCCCATTTCATTACAGCGGTTGACGATGccccc     >  6:56570/1‑146 (MQ=255)                                                                                                                                       gCATAAAAAAACCGGCCACGAAGGCCGGTTTTTTCATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGaa  >  3:157922/1‑149 (MQ=38)                                                                                                                                           gAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttctttc                         >  4:195099/1‑122 (MQ=255)                                                                                                                                           |                                                                                                                                                   CTCTGCAGTTCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAA  >  NC_002947/1243587‑1243869

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