Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	671,319	A→T	10.6%	E125V (GAA→GTA)	PP_0573 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	671,319	0	A	T	10.6%	121.6 / 4.0	45	E125V (GAA→GTA)	PP_0573	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (20/20);  new base T (2/3);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

CGGCTATCGACAACCCCGACTGGGCCCCGGCCAACACCCGGCTGGTGCTGCAAAAAGTGCGCTTCGAGGGCGGCACGGTGTTCCCCTTGGGCGACCTGCGCGAACACTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACTTGCTGTCGTATGCCTACCTGCCTCCTCAGGATGTGGCCGATGGCAGGGTCGACGTGGTGCTGGTCGAAGGTTACGTACGTGATTACCACGTG  >  NC_002947/671175‑671464                                                                                                                                                 |                                                                                                                                                  cggcTATCGACAACCCCGACTGGGCCCCGGCCAACACCCGGCTGGTGCTGCAAAAAGTGCGCTTCGAGGGCGGCACGGTGTTCCCCTTGGGCGACCTGCGCGAACACTATCAGCCGCTGATCGGGCGCTAAACCACGTTGGGTGAATTg                                                                                                                                               <  5:219706/149‑1 (MQ=255) cggcTATCGACAACCCCGACTGGGCCCCGGCCAACACCCGGCTGGTGCTGCAAAAAGTGCGCTTCGAGGGCGGCACGGTGTTCCCCTTGGGCGACCTGCGCGAACACTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTg                                                                                                                                               <  1:186361/149‑1 (MQ=255)  ggcTATCGACAACCCCGACTGGGCCCCGGCCAACACCCGGCTGGTGCTGCAAAAAGTGCGCTTCGAGGGCGGCACGGTGTTCCCCTTGGGCGACCTGCGCGAACACTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTg                                                                                                                                               >  3:40850/1‑148 (MQ=255)               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                                                                                                                          <  3:229839/128‑5 (MQ=39)                                           tGGTGCTGCAAAAAGTGCGCTTCGAGGGCGGCACGGTGTTCCCCTTGGGCGACCTGCGCGAACACTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGTATTGCTGCAATTCACCCAacgtggt                                                                                                                          >  4:229839/1‑124 (MQ=39)                                                ctgcAAAAAGTGCGCTTCGAGGGCGGCACGGTGTTCCCCTTGGGCGACCTGCGCGAACACTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTa                                                                                                >  5:178884/1‑149 (MQ=255)                                                   cAAAAAGTGCGCTTCGAGGGCGGCACGGTGTTCCCCTTGGGCGACCTGCGCGAACACTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACt                                                                                              <  3:344213/148‑1 (MQ=255)                                                   cAAAAAGTGCGCTTCGAGGGCGGCACGGTGTTCCCCTTGGGCGACCTGCGCGAACACTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACt                                                                                              <  3:349661/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                                                                         <  1:12522/130‑1 (MQ=255)                                                                                          ggCGACCTGCGCGAACACTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACTTGCTGTCGTATGCCTACCTGc                                                                         >  2:12522/1‑130 (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)                                                                                                                                                 |                                                                                                                                                  CGGCTATCGACAACCCCGACTGGGCCCCGGCCAACACCCGGCTGGTGCTGCAAAAAGTGCGCTTCGAGGGCGGCACGGTGTTCCCCTTGGGCGACCTGCGCGAACACTATCAGCCGCTGATCGGGCGCGAAACCACGTTGGGTGAATTGCAGCAATACACCCAGCGCCTGACCCAACGCTACCGGCAGAAGGGCTACTTGCTGTCGTATGCCTACCTGCCTCCTCAGGATGTGGCCGATGGCAGGGTCGACGTGGTGCTGGTCGAAGGTTACGTACGTGATTACCACGTG  >  NC_002947/671175‑671464

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