Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	874,539:1	+T	9.3%	coding (115/840 nt)	PP_0759 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	874,539	1	.	T	9.3%	152.5 / 11.2	43	coding (115/840 nt)	PP_0759	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (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 = 4.47e-01

GACTCGGCCTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCG‑TCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGG  >  NC_002947/874403‑874683                                                                                                                                          |                                                                                                                                                 gACTCGGCCTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCG‑TCGAACGCCAACa                                                                                                                                     >  1:3656/1‑150 (MQ=255)             gATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCG‑TCGAACGCCAACAGTACATGTTCAg                                                                                                                         <  1:28433/150‑1 (MQ=255)                  ggCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCG‑TCGAACGCCAACAGTAc                                                                                                                                 >  2:13202/1‑137 (MQ=255)                  ggCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCG‑TCGAACGCCAACAGTAc                                                                                                                                 <  1:13202/137‑1 (MQ=255)                  ggCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCG‑TCGAACGCCAACAGTACATGTTCAGCATGc                                                                                                                    >  1:469304/1‑150 (MQ=255)                      cATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCCGTGAATACCACCAGCGGCGGTGCTGTTGGCG‑TCGAACGCCAACAGTACATGTCcatcatcggcgc                                                                                                                <  2:612795/150‑7 (MQ=255)                       aTGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCG‑TCGAACGCCAACAGTACATGTTCAGCATGCTCTCg                                                                                                               >  2:206059/1‑150 (MQ=255)                        tgtgGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTTGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCG‑TCGAACGCCAACAGTACATGTTCAGCATGCTCTCGa                                                                                                              <  1:513660/150‑1 (MQ=255)                            gATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGAGGTGAATACCACCAGCGGCGGTGCTGTTGGCG‑TCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGa                                                                                                          >  1:272758/1‑150 (MQ=255)                             aTAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCG‑TCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGAt                                                                                                         <  2:656535/150‑1 (MQ=255)                             aTAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCG‑TCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGAt                                                                                                         <  1:58317/150‑1 (MQ=255)                                gggAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCG‑TCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAg                                                                                                      <  2:3656/150‑1 (MQ=255)                                gggAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCG‑TCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAg                                                                                                      <  2:469304/150‑1 (MQ=255)                                 ggAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTNCTGAGTGCTNGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCG‑TCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAgg                                                                                                     >  1:395786/1‑150 (MQ=255)                                 ggAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCG‑TCGAACGCCaa                                                                                                                                       <  2:355583/116‑1 (MQ=255)                                 ggAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCG‑TCGAACGCCaa                                                                                                                                       >  1:355583/1‑116 (MQ=255)                                 ggAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCG‑TCGAACGCCAACAGTACa                                                                                                                                <  2:218806/123‑1 (MQ=255)                                 ggAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCG‑TCGAACGCCAACAGTACa                                                                                                                                >  1:218806/1‑123 (MQ=255)                                        tttCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCG‑TCGAACGCCAACAGTACATGTCCAGCATGCCGACTACCCATGAGGGCAAccg                                                                                              <  2:503343/150‑2 (MQ=255)                                          tCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCG‑TCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGa                                                                                            <  1:127377/150‑1 (MQ=255)                                            aTGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTACTGTTGGCGTTCG‑ACGCCAAc                                                                                                                                      >  2:208919/1‑106 (MQ=255)                                            aTGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTACTGTTGGCGTTCG‑ACGCCAAc                                                                                                                                      <  1:208919/106‑1 (MQ=255)                                                 tAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGAGGTGAATACCACCAGCGGCGGTGCTGTTGGCG‑TCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTAcgc                                                                                     <  2:272758/150‑1 (MQ=255)                                                         ttGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCG‑TCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGt                                                                                         >  1:729414/1‑138 (MQ=255)                                                         ttGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCG‑TCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGt                                                                                         <  2:729414/138‑1 (MQ=255)                                                               tcaGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCG‑TCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATcagc                                                                       >  1:303377/1‑150 (MQ=255)                                                                   ccTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCG‑TCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGAc                                                                   >  2:100374/1‑150 (MQ=255)                                                                      gctgAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCG‑TCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCCGTGGTATCAGCAAACCCt                                                                >  2:581599/1‑150 (MQ=255)                                                                               tGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCG‑TCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGc                                                       >  2:237170/1‑150 (MQ=255)                                                                                cgCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCG‑TCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCa                                                      <  2:395786/149‑1 (MQ=255)                                                                                    tCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCG‑TCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGa                                                  >  2:623660/1‑150 (MQ=255)                                                                                                         cGGTGAATACCACCAGCGGCGGTGCTGTTGGCG‑TCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGt                             <  1:40619/150‑1 (MQ=255)                                                                                                             gAATACCACCAGCGGCGGTGCTGTTGGCG‑TCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAg                         >  2:210626/1‑150 (MQ=255)                                                                                                                  ccaccaGCGGCGGTGCTGTTGGCG‑TCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAg                                                                                 >  1:112594/1‑89 (MQ=255)                                                                                                                  ccaccaGCGGCGGTGCTGTTGGCG‑TCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAg                                                                                 <  2:112594/89‑1 (MQ=255)                                                                                                                               tGCTGTTGGCG‑TCGAACGCCAACAGTACATGTTCAGCATGCTGTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGAACCTGGAGGAAGTCTCGACCAAGGGGGTTCGCGACAAGGCCAGCGCCGATGCCAAGCgtgt       >  2:611531/1‑150 (MQ=255)                                                                                                                                gCTGTTGGCG‑TCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGtgtg      <  2:818690/150‑1 (MQ=255)                                                                                                                                  tGTTGGCGTTCG‑ACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGaa                                                         <  2:160254/97‑1 (MQ=255)                                                                                                                                  tGTTGGCGTTCG‑ACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGaa                                                         >  1:160254/1‑97 (MQ=255)                                                                                                                                  tGTTGGCG‑TTGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCa    >  2:583472/1‑150 (MQ=255)                                                                                                                                  tGTTGGCG‑TTCGACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGCGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCACGGGTGTGCTCGACAACTCCAGCGCCGATGCCAAGCGTGTGCa    >  2:536279/1‑150 (MQ=255)                                                                                                                                   gTTGGCG‑TCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAg   <  1:100374/150‑1 (MQ=255)                                                                                                                                    ttGGCG‑TCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAgg  <  1:210626/150‑1 (MQ=255)                                                                                                                                          |                                                                                                                                                 GACTCGGCCTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCG‑TCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGG  >  NC_002947/874403‑874683

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