Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	847,680	T→C	9.5%	N297S (AAC→AGC)	PP_0730 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	847,680	0	T	C	9.5%	111.7 / 2.3	42	N297S (AAC→AGC)	PP_0730	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (20/18);  new base C (2/2);  total (22/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 = 9.95e-01

TGCGCGCTGAAGCACTGAGTCGTCAATTAGAATGGGCTCGCCCAGATGCAGACGGACTAGCAGTTCGAACGGGTCGATGCCGGTCAGCATCCGGTACATCTGCGGCGCTACGCTACCCATCATGCGTGCGTTGATCTCCACTAGACAGGGCCCTTTGGCATCGCGCAACAGTTCAACATGATAAAGGCCAACATCCAGCCCCAGAGCGCTGAACACTGCCTGCACATATGCTTTGATCTCGTCGCGCTGTGCTGCGTCAAGGCCTGACGGCATGGTCG  >  NC_002947/847550‑847827                                                                                                                                   |                                                                                                                                                    tGCGCGCTGAAGCACTGAGTCGTCAATTAGAATGGGCTCGCCCAGATGCAGACGGACTAGCAGTTCGAACGGGTCGATGCCGGTCAGCATCCGGTACATCTGCGGCGCTACGCTACCCATCATGCGTGCGTTGATCTCCACTAGACAgg                                                                                                                                   >  1:314388/1‑149 (MQ=255)   cgcgcTGAATCACTGAGTCGTCAATTAGACTGGGCTCGCCCAGATGCAGACGGACTAGCAGTTCGAACGGGTCGATGCCGGTCAGCATCCGGTACATCTGCGGCGCTACGCTACCCATCATGCGTGCGTTGATCTCCACTAGACAGGGc                                                                                                                                 <  4:38914/149‑1 (MQ=255)        tGAAGCACTGAGTCGTCAATTAGAATGGGCTCGCCCAGATGCAGACGGACTAGCAGTTCGAACGGGTCGATGCCGGTCAACATCCGGTACATCTGCGGCGCTACGCTACACATCAAGCGTGCGTTGATCTCCACTAGACAGGGCCCttt                                                                                                                            >  4:113733/1‑149 (MQ=255)              aCTGAGTCGTCAATTAGAATGGGCTCGCCCAGATGCAGACGGACTAGCAGTTCGAACGGGTCGATGCCGGTCAGCATCCGGTACATCTGCGGCGCTACGCTACCCATCATGCGTGCGTTGATCTCCACTAGACAGGGCCCTTTGGCa                                                                                                                        >  4:164094/1‑147 (MQ=255)                      gtcAATTAGAATGGGCTCGCCCAGATGCAGACGGACTAGCAGTTCGAACGGGTCGATGCCGGTCAGCATCCGGTACATCTGCGGCGCTACGCTACCCATCATGCGTGCGTTGATCTCCACTAGACAGGGCCCTTTGGCATCGCGCAAc                                                                                                               <  3:164094/148‑1 (MQ=255)                       tcAATTAGAATGGGCTCGCCCAGATGCAGACGGACTAGCAGTTCGAACGGGTCGATGCCGGTCAGCATCCGGTACATCTGCGGCGCTACGCTACCCATCATGCGTGCGTTGATCTCCACTAGACAGGGCCCTTTGGCATCGCGCAACAg                                                                                                             >  3:289829/1‑149 (MQ=255)                            tAGAATGGGCTCGCCCAGATGCAGACGGACTAGCAGTTCGAACGGGTCGATGCCGGTCAGCATCCGGTACATCTGCGGCGCTACGCTACCCATCATGCGTGCGTTGATCTCCACTAGACAGGGCCCTTTGGCATCGCGCAACAGTTCaa                                                                                                        >  8:55910/1‑149 (MQ=255)                              gAATGGGCTCGCCCAGATGCAGACGGACTAGCAGTTCGAACGGGTCGATGCCGGTCAGCATCCGGTACATCTGCGGCGCTACGCTACCCATCATGCGTGCGTTGATCTCCACTAGACAGGGCCCTTTGGCATCGCGCAACAGTTCAAc                                                                                                       <  2:314388/148‑1 (MQ=255)                               aaTGGGCTCGCCCAGATGCAGACGGACTAGCAGTTCGAACGGGTCGATGCCGGTCAGCATCCGGTACATCTGCGGCGCTACGCTACCCATCATGCGTGCGTTGATCTCCACTAGACAGGGCCCTTTGGCATc                                                                                                                      <  1:214287/132‑1 (MQ=255)                               aaTGGGCTCGCCCAGATGCAGACGGACTAGCAGTTCGAACGGGTCGATGCCGGTCAGCATCCGGTACATCTGCGGCGCTACGCTACCCATCATGCGTGCGTTGATCTCCACTAGACAGGGCCCTTTGGCATc                                                                                                                      >  2:214287/1‑132 (MQ=255)                                                  aGACGGACTAGCAGTTCGAACGGGTCGATGCCGGTCAGCATCCGGTACATCTGCGGCGCTACGCTACCCATCATGCGTGCGTTGAtctc                                                                                                                                              <  7:136584/89‑1 (MQ=255)                                                  aGACGGACTAGCAGTTCGAACGGGTCGATGCCGGTCAGCATCCGGTACATCTGCGGCGCTACGCTACCCATCATGCGTGCGTTGAtctc                                                                                                                                              >  8:136584/1‑89 (MQ=255)                                                   gACGGACTAGCAGTTCGAACGGGTCGATGCCGGTCAGCATCCGGTACATCTGCGGCGCTACGCTACCCATCATGCGTGCGTTGATCTCCACTAGACAGGGCCCTTTGGCATCGCGCAACAGTTCAACATGATAAAGGCCAACATCCAGc                                                                                 <  4:289829/149‑1 (MQ=255)                                                    aCGGACTAGCAGTTCGAACGGGTCGATGCCGGTCAGCATCCGGTACATCTGCGGCGCTACGCTACCCATCATGCGTGCGTTGATCTCCACTAGACAGGGCCCTTTGGCATCGCGCAACAGTTCAACATGATAAAGGCCAACATCCAGcc                                                                                >  3:278778/1‑149 (MQ=255)                                                     cGGACTAGCAGTTCGAACGGGTCGATGCCGGTCAGCATCCGGTACATCTGCGGCGCTACGCTACCCATCATGCGTGCGTTGATCTCCACTAGACAGGGCCCTTTGGCATCGCGCAACAGTTCAACATGATAAAGGCCAACATCCAGccc                                                                               <  6:188629/149‑1 (MQ=255)                                                         cTAGCAGTTCGAACGGGTCGATGCCGGTCAGCATCCGGTACATCTGCGGCGCTACGCTACCCATCATGCGTGCGTTGATCTCCACTAGACAGGGCCCTTTGGCATCGCGCAACAGTTCAACATGATAAAGGCCAACATCCAGCCCCaga                                                                           <  7:55910/149‑1 (MQ=255)                                                         cTAGCAGTTCGAACGGGTCGATGCCGGTCAGCATCCGGTACATCTGCGGCGCTACGCAACCCATCATGCGTGCGTTGATCTCCACTAGACAGGGCCCTTTGGCATCGGGCAACAGTTCAACATGATAAAGGCCAACATCCAGCCCCaga                                                                           >  4:108750/1‑149 (MQ=255)                                                          tAGCAGTTCGAACGGGTCGATGCCGGTCAGCATCCGGTACATCTGCGGCGCTACGCTACCCATCATGCGTGCGTTGATCTCCACTAGACAGGGCCCTTTGGCATCGCGCAACAGTTCAACATGATAAAGGCCAACATCCAGCCCCagag                                                                          >  3:277130/1‑149 (MQ=255)                                                            gCAGTTCGAACGGGTCGATGCCGGTCAGCATCCGGTACATCTGCGGCGCTACGCTACCCATCATGCGTGCGTTGATCTCCACTAGACAGGGCCCTTTGGCATCGCGCAACAGTTCAACATGATAAAGGCCAACATCCAGCCCCAGAgcg                                                                        >  5:154867/1‑149 (MQ=255)                                                            gCAGTTCGAACGGGTCGATGCCGGTCAGCATCCGGTACATCTGCGGCGCTACGCTACCCATCATGCGTGCGTTGATCTCCACTAGACAGGGCCCTTTGGCATCGCGCAACAGTTCAACATGATAAAGGCCAACATCCAGCCCCAGAgcg                                                                        >  6:26684/1‑149 (MQ=255)                                                              aGTTCGAACGGGTCGATGCCGGTCAACATCCGGTACATCTGCGGCGCTACGCTACCCATCTTGCGTGCGTTGATCTCCACTAGACAGGGCCCTTTGGCATCGCGCAACATTTCAACATGATAAAGGCCAACATCCAGCCCCAGAGCGCt                                                                      <  3:113733/149‑1 (MQ=255)                                                               gTTCGAACGGGTCGATGCCGGTCAGCATCCGGTACATCTGTGGCGCTACGCTACCCATCATGCGTGCGCTGATCTCCACTAGACAGGGCCCTTTGGCATCGCGCAACAGTTCAACATGATAAAGGCCAACATCCAGCCCCAGAGCGCt                                                                      >  2:233380/1‑148 (MQ=255)                                                               gTTCGAACGGGTCGATGCCGGTCAGCATCCGGTACATCTGTGGCGCTACGCTACCCATCATGCGTGCGCTGATCTCCACTAGACAGGGCCCTTTGGCATCGCGCAACAGTTCAACATGATAAAGGCCAACATCCAGCCCCAGAGCGCt                                                                      >  2:237935/1‑148 (MQ=255)                                                                ttCGAACGGGTCGATGCCGGTCAGCATCCGGTACATCTGCGGCGCTACGCTACCCATCATGCGTGCGTTGATCTCCACTAGACAGGGCCCTTTGGCATCGCGCAACAGTTCAACATGATAAAGGCCAACATCCAGCCCCAGAGCGCTg                                                                     >  3:82824/1‑148 (MQ=255)                                                                    aaCGGGTCGATGCCGGTCAGCATCCGGTACATCTGCGGCGCTACGCTACCCATCATGCGTGCGTTGATCTCCACTAGACAGGGCCCTTTGGCATCGCGCAACAGTTCAACATGATAAAGGCCAACATCCAGCCCCAGAGCGCTGAacac                                                                >  8:270361/1‑149 (MQ=255)                                                                         gTCGATGCCGGTCAGCATCCGGTACATCTGCGGCGCTACGCTACCCATCATGCGTGCGTTGATCTCCACTAGACAGGGCCCTTTGGCATCGCGCAACAGTTCAACATGATAAAGGCCAACATACAGCCCCAGAGCGCTGAACActgcct                                                           >  5:124074/1‑149 (MQ=255)                                                                         gTCGATGCCGGTCAGCATCCGGTACATCGGCGGCGCTACGCTACCCATCATGCGTGCGTTGATCTCCACTAGACAGGGCACTTTGGCATCGCGCAACAGTTCAACATGATAAAGGCCAACATCCAGCCCCAGAGCGCTGAACActgcct                                                           >  1:20884/1‑149 (MQ=255)                                                                                ccGGTCAGCATCCGGTACATCTGCGGCGCTACGCTACCCATCATGCGTGCGTTGATCTCCACTAGACAGGGCCCTTTGGCATCGCGCAACAGTTCAACATGATAAAGGCCAACATCCAGCCCCAGAGCGCTGAACACTGCCTGcac                                                       <  5:325141/146‑1 (MQ=255)                                                                                ccGGTCAGCATCCGGTACATCTGCGGCGCTACGCTACCCATCATGCGTGCGTTGATCTCCACTAGACAGGGCCCTTTGGCATCGCGCAACAGTTCAACATGATAAAGGCCAACATCCAGCCCCAGAGCGCTGAACACTGCCTGcac                                                       >  6:325141/1‑146 (MQ=255)                                                                                   gTCAGCATCCGGTACATCTGCGGCGCTACGCTACCCATCATGCGTGCGTTGATCTCCACTAGACAGGGCCCTTTGGCATCGCGCAACAGTTCAACATGATAAAGGCCAACATCCAg                                                                                  >  8:239440/1‑116 (MQ=255)                                                                                   gTCAGCATCCGGTACATCTGCGGCGCTACGCTACCCATCATGCGTGCGTTGATCTCCACTAGACAGGGCCCTTTGGCATCGCGCAACAGTTCAACATGATAAAGGCCAACATCCAg                                                                                  <  7:239440/116‑1 (MQ=255)                                                                                       gCATCCGGTACATCTGCGGCGCTACGCTACCCATCATGCGTGCGTTGATCTCCACTAGACAGGGCCCTTTGGCATCGCGCAACAGTTCAACATGATAAAGGCCAACATCCAGCCCCAGAGCGCTGAACACTGCCTGCACATATGCTTTg                                             <  8:177098/149‑1 (MQ=255)                                                                                        cATCCGGTACATCTGCGGCGCTACGCTACCCATCATGCGTGCGTTGATCTCCACTAGACAGGGCCCTTTGGCATCGCGCAACAGTTCAACATGATAAAGGCCAACATCCAGCCCCAGAGCGCTGAACACTGCCTGCACATATGCTTTGa                                            >  6:317462/1‑149 (MQ=255)                                                                                               tACATCTGCGGCGCTACGCTACCCATCATGCGTGCGTTGATCTCCACTAGACAGGGCCCTTTGGCATCGCGCAACAGTTCAACATGATAAAGGCCAACATCCAGCCCCAGAGCGCTGAACACTGCCTGCACATATGCTTTGATCtcgtc                                     <  6:154867/149‑1 (MQ=255)                                                                                               tACATCTGCGGCGCTACGCTACCCATCATGCGTGCGTTGATCTCCACTAGACAGGGCCCTTTGGCATCGCGCAACAGTTCAACATGATAAAGGCCAACATCCAGCCCCAGAGCGCTGAACACTGCCTGCACATATGCTTTGATCtcgtc                                     <  6:124074/149‑1 (MQ=255)                                                                                                       cggcgCTACGCTACCCATCATGCGTGCGTTGATCTCCACTAGACAGGGCCCTTTGGCATCGCGCAACAGTTCAACATGATAAAGGCCAACATCCAGCCCCAGAGCGCTGAACACTGCCTGCACATATGCTTTGATCTCGTCGCGCtgt                              <  5:317462/148‑1 (MQ=255)                                                                                                           gctacgctacCCATCATGCGTGCGTTGATCTCCACTAGACAGGGCCCTTTGGCATCGCGCAACAGTTCAACATGATAAAGGCCAACATCCAGCCCCAGAGCGCTGAACACTGCCTGCACATATGCTTTGATCTCGTCGCGCTGtgctgc                         <  5:26684/149‑1 (MQ=255)                                                                                                                  tacCCATCATGCGTGCGTTGATCTCCACTAGACAGGGCCCTTTGGCATCGCGCAACAGTTCAACATGATAAAGGCCAACATCCAGCCCCAGAGCGCTGAACACTGCCTGCACATATGCTTTGAtctc                                        >  8:67366/1‑127 (MQ=255)                                                                                                                  tacCCATCATGCGTGCGTTGATCTCCACTAGACAGGGCCCTTTGGCATCGCGCAACAGTTCAACATGATAAAGGCCAACATCCAGCCCCAGAGCGCTGAACACTGCCTGCACATATGCTTTGAtctc                                        <  7:67366/127‑1 (MQ=255)                                                                                                                    cccATCATGCGTGCGTTGATCTCCACTAGACAGGGCCCTTTGGCATCGCGCAACAGTTCAACATGATAAAGGCCAACATCCAGCCCCAGAGCGCTGAACACTGCCTGCACATATGCTTTGATCTCGTCGCGCTGTGCTGCGTCAAgg                  <  3:108750/147‑1 (MQ=255)                                                                                                                    cccATCATGCGTGCGCTGATCTCCACTAGACAGGGCCCTTTGGCATCGCGCAACAGTTCAACATGATAAAGGCCAACATCCAGCCCCAGAGCGCTGAACACTGCCTGCACATATGCTTTGATCTCGTCGCGCTGTGCTGCGTCAAgg                  <  1:237935/147‑1 (MQ=255)                                                                                                                    cccATCATGCGTGCGCTGATCTCCACTAGACAGGGCCCTTTGGCATCGCGCAACAGTTCAACATGATAAAGGCCAACATCCAGCCCCAGAGCGCTGAACACTGCCTGCACATATGCTTTGATCTCGTCGCGCTGTGCTGCGTCAAgg                  <  1:233380/147‑1 (MQ=255)                                                                                                                             cgtgcgtTGATCTCCACTAGACAGGGCCCTTTGGCATCGCGCAACAGTTCAACATGATAAAGGCCAACATCCAGCCCCAGAGCGCTGAACACTGCCTGCACATATGCTTTGATCTCGTCGCGCTGTGCTGCGTCAAGGCCTGACGGCa        <  4:82824/148‑1 (MQ=255)                                                                                                                              gtgcgtTGATCTCCACTAGACAGGGCCCTTTGGCATCGCGCAACAGTTCAACATGATAAAGGCCAACATCCAGCCCCAGAGCGCTGAACACTGCCTGCACATATGCTTTGATCTCGTCGCGCTGTGCTGCGTCAAGGCCTGACGGCATg      <  4:277130/149‑1 (MQ=255)                                                                                                                                  gtTGATCTCCACTAGACAGGGCCCTTTGGCATCGCGCAACAGTTCAACATGATAAAGGCCAACATCCAGCCCCAGAGCGCTGAACACTGCCTGCACATATGCTTTGATCTCGTCGCGCTGTGCTGCGTCAAGGCCTGACGGCATGGTcg  >  6:217086/1‑149 (MQ=255)                                                                                                                                   |                                                                                                                                                    TGCGCGCTGAAGCACTGAGTCGTCAATTAGAATGGGCTCGCCCAGATGCAGACGGACTAGCAGTTCGAACGGGTCGATGCCGGTCAGCATCCGGTACATCTGCGGCGCTACGCTACCCATCATGCGTGCGTTGATCTCCACTAGACAGGGCCCTTTGGCATCGCGCAACAGTTCAACATGATAAAGGCCAACATCCAGCCCCAGAGCGCTGAACACTGCCTGCACATATGCTTTGATCTCGTCGCGCTGTGCTGCGTCAAGGCCTGACGGCATGGTCG  >  NC_002947/847550‑847827

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