Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	671,325	A→T	9.3%	Q127L (CAG→CTG)	PP_0573 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	671,325	0	A	T	9.3%	111.1 / 3.4	43	Q127L (CAG→CTG)	PP_0573	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (19/20);  new base T (2/2);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

CCCGACTGGGCCCCGGCCAACACCCGGCTGGTGCTGCAAAAAGTGCGCTTCGAGGGCGGCACGGTGTTCCCCTTGGGCGACCTGCGCGAACACTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACTTGCTGTCGTATGCCTACCTGCCTCCTCAGGATGTGGCCGATGGCAGGGTCGACGTGGTGCTGGTCGAAGGTTACGTACGTGATTACCACGTGGAAGGCGA  >  NC_002947/671189‑671472                                                                                                                                         |                                                                                                                                                    cccGACTGGGCCCCGGCCAACACCCGGCTGGTGCTGCAAAAAGTGCGCTTCGAGGGCGGCACGGTGTTCCCCTTGGGCGACCTGCGCGAACACTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCa                                                                                                                                         >  4:45654/1‑149 (MQ=255)   cGACTGGGCCCCGGCCAACACCCGGCTGGTGCTGCAAAAAGTGCGCTTCGAGGGCGGCACGGTGTTCCCCTTGGGCGACCTGCGCGAACACTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAgc                                                                                                                                       <  5:307804/149‑1 (MQ=255)    gACTGGGCCCCGGCCAACACCCGGCTGGTGCTGCAAAAAGTGCGCTTCGAGGGCGGCACGGTGTTCCCCTTGGGCGACCTGCGCGAACACTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGTATTGCTGCAATTCACCCAacg                                                                                                                                      <  4:230924/149‑1 (MQ=255)         ggCCCCGGCCAACACCCGGCTGGTGCTGCAAAAAGTGCTCTTCGAGGGCGGCACGGTGTTCCCCGTGGGCGACCTGCGCGAACACTATCAGCCGCTGATCGGGCCCGAAACCACGTTGGGTGAATTGCAGCACTACACCCAGCGCCAGa                                                                                                                                 >  5:242594/1‑149 (MQ=255)                        ccGGCTGGTGCTGCAAATAGTGCGCTTCGAGGGCGGCACGGTGTTCCCCTTGGGCGACCTGCGCGAACACTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGATCCAACGCTACCGGc                                                                                                                  <  7:316564/149‑1 (MQ=255)                             tGGTGCTGCAAAAAGTGCGCTTCGAGGGCGGCACGGTGTTCCCCTTGGGCGACCTGCGCGAACACTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGTATTGCTGCAATTCACCCAacgtggt                                                                                                                                  >  4:229839/1‑124 (MQ=39)                             tGGTGCTGCAAAAAGTGCGCTTCGAGGGCGGCACGGTGTTCCCCTTGGGCGACCTGCGCGAACACTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGTATTGCTGCAATTCACCCAacgtggt                                                                                                                                  <  3:229839/128‑5 (MQ=39)                                  ctgcAAAAAGTGCGCTTCGAGGGCGGCACGGTGTTCCCCTTGGGCGACCTGCGCGAACACTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTa                                                                                                        >  5:178884/1‑149 (MQ=255)                                     cAAAAAGTGCGCTTCGAGGGCGGCACGGTGTTCCCCTTGGGCGACCTGCGCGAACACTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACt                                                                                                      <  3:349661/148‑1 (MQ=255)                                     cAAAAAGTGCGCTTCGAGGGCGGCACGGTGTTCCCCTTGGGCGACCTGCGCGAACACTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACt                                                                                                      <  3:344213/148‑1 (MQ=255)                                      aaaaaGTGCGCTTCGAGGGCGGCACGGTGTTCCCCTTGGGCGACCTGCGCGAACACTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACTTg                                                                                                    >  7:254778/1‑149 (MQ=255)                                         aaGTGCGCTTCGAGGGCGGCACGGTGTTCCCCTTGGGCGACCTGCGCGAACACTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACTTGCTg                                                                                                 <  4:169820/149‑1 (MQ=255)                                              cgcTTCGAGGGCGGCACGGTGTTCCCCTTGGGCGACCTGCGCGAACACTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACTTGCTGTCGTa                                                                                            <  4:226269/149‑1 (MQ=255)                                               gcTTCGAGGGCGGCACGGTGTTCCCCTTGGGCGACCTGCGCGAACACTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACTTGCTGTCGTAt                                                                                           <  7:301173/149‑1 (MQ=255)                                                       ggcggcACGGTGTTCCCCTTGGGCGACCTGCGCGAACACTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACTTGCTGTCGTATGCCTACCt                                                                                   <  4:40850/149‑1 (MQ=255)                                                                      cccTTGGGCGACCTGCGCGAACACTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACTTGCTGTCGTATGCCTACCTGCCTCCTCAGGAtgt                                                                    <  3:45654/149‑1 (MQ=255)                                                                       gcTTGGGCTAGCTGCGCGAACACTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGAGCCAACGCTACCGGCAGAAGGGCTACTTGCTGTCGTATGCCTACCTGCCTCCTCTGGAtgtg                                                                   <  6:242594/148‑1 (MQ=255)                                                                       ccTTGGGCGACCTGCGCGAACACTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACTTGCTGTCGTATGCCTACCTGc                                                                                 <  6:75211/135‑1 (MQ=255)                                                                       ccTTGGGCGAACTGCGCGAACACTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACTTGCTGTCGTATGCCTACCTGc                                                                                 >  5:75211/1‑135 (MQ=255)                                                                        cTTGGGCGACCTGCGCGAACACTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACTTGCTGTCGTATGCCTACCTGCCTCCTCAGGATGTgg                                                                  >  2:196049/1‑149 (MQ=255)                                                                            ggCGACCTGCGCGAACACTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACTTGCTGTCGTATGCCTACCTGc                                                                                 >  2:12522/1‑130 (MQ=255)                                                                            ggCGACCTGCGCGAACACTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACTTGCTGTCGTATGCCTACCTGc                                                                                 <  1:12522/130‑1 (MQ=255)                                                                            ggCGACCTGCGCGAACACTATAAGCAGATGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACCCCCAGCGCCTGACCCAACGCAACCGGCAGAAGGGCTACTTGCTGTCGTATGCCTACCTGCCTCCTGAGGATGTGGACGa                                                              >  6:111930/1‑149 (MQ=255)                                                                             gCGACCTGCGCGAACACTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACTTGCTGTCGTATGCCTACCTGCCTCCTCAGGATGTgg                                                                  <  8:278694/144‑1 (MQ=255)                                                                             gCGACCTGCGCGAACACTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACTTGCTGTCGTATGCCTACCTGCCTCCTCAGGATGTgg                                                                  >  7:278694/1‑144 (MQ=255)                                                                              cGACCTGCGCGAACACTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACTTGCTGTCGTATGCCTACCTGCCTCCTCAGGATGTGGCCGATg                                                            <  6:178884/149‑1 (MQ=255)                                                                                aCCTGCGCGAACACTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACTTGCTGTCGTATGCCTACCTGCCTCCTCAGGATGTGGCCGATGGc                                                          <  1:196049/149‑1 (MQ=255)                                                                                          acacTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACTTGCTGTCGTATGCCTACCTGCCTCCTCAGGATGTGGCCGATGGCAGGGTCGACg                                                >  4:123734/1‑149 (MQ=255)                                                                                            acTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACTTGCTGTCGTATGCCTACCTGCCTCCTCAGGATGTGGCCGATGGCAGGGTTGACgtg                                              <  8:254778/149‑1 (MQ=255)                                                                                             cTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACTTGCTGTCGTATGCCTACCTGCCTCCTCAGGATGTGGCCGATGGCAGGGTCGACgtgg                                             >  3:175793/1‑149 (MQ=255)                                                                                              tATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACTTGCTGTCGTATGCCTACCTGCCTCCTCAGGATGTGGCCGATGGCAGGGTCGACgtggt                                            >  1:234594/1‑149 (MQ=255)                                                                                                tCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACTTGCTGTCGTATGCCTACCTGCCTCCTCAGGATGTGGCCGATGGCAGGGTCGACGTGGTGc                                          <  4:175793/149‑1 (MQ=255)                                                                                                  aGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACTTGCTGTCGTATGCCTACCTGCCTCCTCAGGATGTGGCCGATGGCAGGTTCGACGTGGTGCTg                                        >  5:38091/1‑149 (MQ=255)                                                                                                   gCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGCCAGAAGGGCTACTTGCTGTCGTATGCCTACCTCCCTCCTCAGGATGTGGCCGATGGCAGGGTCGACGTGGTGCTgg                                       >  7:114840/1‑149 (MQ=255)                                                                                                    ccGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACTTGCTGTCGTATGCCTACCTGCCTCCTCAGGATGGGGCCGATGGCAGGGTCGACGTGCTGCTGGt                                      >  5:100592/1‑149 (MQ=255)                                                                                                       cTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACTTGCTGTCGTATGCCTACCTGCCTCCTCAGGATGTGGCCGATGGCAGGGTCGACGTGGTGCTGGTCg                                    >  2:206130/1‑148 (MQ=255)                                                                                                          aTCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACTTGCTGTCGTATGCCTACCTGCCTCCTCAGGATGTGGCCGATGGCAGGGTCGACGTGGTGCTGGTCGAAgg                                <  3:123734/149‑1 (MQ=255)                                                                                                                ggcgACACCTCGTTGGGTGTATTGGAGCACTACAGCCAGCGCATGACGCAACGCTACGGGCTGAAGGGCTACTTGCTGTGGTATGCCTACCTGCCTCCTCAGGATGTGGCCTTTGGCAGGGTCGACGTGGTGCTGGTCGAAGGTtacgt                          <  6:100592/148‑1 (MQ=255)                                                                                                                 gcgAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACTTGCTGTCGTATGCCTACCTGCCTCCTCAGGATGTGGCCGATGGCAGGGTCGACGTGGTGCTGGTCGAAGGTtacgta                         >  2:173866/1‑149 (MQ=255)                                                                                                                   gAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACTTGCTGTCGTATGCCTACCTGCCTCCTCAGGATGTGGCCGATGGCAGGGTCGACGTGGTGCTGGTCGAAGGTtacgtacg                       >  8:27120/1‑149 (MQ=255)                                                                                                                             tGGGTGTATTGCTGCAATTCACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACTTGCTGTCGTATGCCTACCTGCCTCCTCAGGATGTGGCCGATGGCAGGGTCGACGTGGTGCTGGTCGAAGGTTACGTACGTGATTACCAc             >  2:44345/1‑149 (MQ=255)                                                                                                                                 tGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACTTGCTGTCGTATGCCTACCTGCCTCCTCAGGATGTGGCCGATGGCAGGGTCGACGTGGTGCTGGTCGAAGGTTACGTACGTGATTACCAGGTg          >  8:171431/1‑148 (MQ=255)                                                                                                                                    aTTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACTTGCTGTCGTATGCCTACCTGCCTCCTCAGGATGTGGCCGATGGCAGGGTCGACGTGGTGCTGTTCGAAGGTTACGTACGTGATTACCACGTGGaa       <  6:38091/148‑1 (MQ=255)                                                                                                                                         agcaATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACTTGCTGTCGTATGCCTACCTGCCTCCTCAGGATGTGGCCGATGGCAGGGTCGACGTGGTGCTGGTCGAAGGTTACGTACGTGATTAc                <  7:129822/134‑1 (MQ=255)                                                                                                                                         agcaATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACTTGCTGTCGTATGCCTACCTGCCTCCTCAGGATGTGGCCGATGGCAGGGTCGACGTGGTGCTGGTCGAAGGTTACGTACGTGATTAc                >  8:129822/1‑134 (MQ=255)                                                                                                                                         agcaATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACTTGCTGTCGTATGCCTACCTGCCTCCTCAGGATGTGGCCGATGGCAGGGTCGACGTGGTGCTGGTCGAAGGTTACGTACGTGATTACAAAGTGGAAGGCGa  >  6:186214/1‑148 (MQ=255)                                                                                                                                         |                                                                                                                                                    CCCGACTGGGCCCCGGCCAACACCCGGCTGGTGCTGCAAAAAGTGCGCTTCGAGGGCGGCACGGTGTTCCCCTTGGGCGACCTGCGCGAACACTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACTTGCTGTCGTATGCCTACCTGCCTCCTCAGGATGTGGCCGATGGCAGGGTCGACGTGGTGCTGGTCGAAGGTTACGTACGTGATTACCACGTGGAAGGCGA  >  NC_002947/671189‑671472

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