Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	672,098	G→A	9.0%	A385T (GCG→ACG)	PP_0573 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	672,098	0	G	A	9.0%	122.3 / 2.4	44	A385T (GCG→ACG)	PP_0573	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (19/21);  new base A (2/2);  total (21/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 = 1.00e+00

TGCGATCGGCCTGCGTCAGCCGGTGGTTGTCAGGCCGAATGAGTGGCTGGCGGTGCAGGGCCACCTGTACTCGGTCAGTGAGCATGATGAAGACCGGATGGCCGAGCAGCCAACGTCTCGCAAATACGACACCTATGTGCGGGCGCTGTCATTCGAAGGTGACTGGCGCAAGGTGGAGCAGGGGCGTTTGCGTATCGTCAGTGCCGGGGTCCATCAGGGTCTGGACTACCTGGGGGCACGCAGCGATGCTGACTACGATCTCGACTTCCTGCGCTTGCGGCTTTCCGG  >  NC_002947/671956‑672243                                                                                                                                               |                                                                                                                                                  tGCGATCGGCCTGCGTCAGCCGGTGGTTGTCAGGCCGAATGAGTGGCTGGCGGTGCAGGGCCACCTGTACTCGGTCAGTGAGCATGATGAAGACCGGATGGCCGAGCAGCCAACGTCTCGCAAATACGACACCTATGTGCGGGCGCt                                                                                                                                               <  6:87811/147‑1 (MQ=255) tGCGATCGGCCTGCGTCAGCCGGTGGTTGTCAGGCCGAATGAGTGGCTGGCGGTGCAGGGCCACCTGTACTCGGTCAGTGAGCATGATGAAGACCGGATGGCCGAGCAGCCAACGTCTCGCAAATACGACACCTATGTGCGGGCGCTg                                                                                                                                              <  5:110157/148‑1 (MQ=255) tGCGATCGGCCTGCGTCAGCCGGTGGTTGTCAGGCCGAATGAGTGGCTGGCGGTGCAGGGCCACCTGTACTCGGTCAGTGAGCATGATGAAGACCGGATGGCCGAGCAGCCAACGTCTCGCAAATACGACACCTATGTGCGGGCGCTGt                                                                                                                                             >  5:6807/1‑149 (MQ=255)                    ccGGTGGTTGTCAGGCCGAATGAGTGGCTGGCGGTGCAGGGCCACCTGTACTCGGTCAGTGAGCATGATGAAGACCGGATGGCCGAGCAGCCAACGTCTCGCAAATACGACACCTATGTGCGGGCGCTGTCATTCGAAGGTGACTGgcg                                                                                                                          >  5:162286/1‑149 (MQ=255)                          gtTGTCAGGCCGAATGAGTGGCTGGCGGTGCAGGGCCACCTGTACTCGGTCAGTGAGCATGATGAAGACCGGATGGCCGAGCAGCCAACGTCTCGCAAATACGACACCTATGTGCGGGCGCTGTCATTCGAAGGTGACTGGCGCAAgg                                                                                                                     <  6:6807/148‑1 (MQ=255)                                 ggCCGAATGAGTGGCTGGCGGTGCAGGGCCACCTGTACTCGGTCAGTGAGCATGATGAAGACCGGATGGCCGAGCAGCCAACGTCTCGCAAATACGACACCTATGTGCGGGCGCTGTCATTCGAAGGTGACTGGCGCAAGGTGGAGCAg                                                                                                             <  7:68572/149‑1 (MQ=255)                                        tGAGTGGCTGGCGGTGCAGGGCCACCTGTACTCGGTCAGTGAGCATGATGAAGACCGGATGGCCGAGCAGCCAACGTCTCGCAAATACGACACCTATGTGCGGGCGCTGTCATTCGAAGGTGACTGGCGCAAGGTGGAGCAGGGGCGtt                                                                                                      >  2:155907/1‑149 (MQ=255)                                          aGTGGCTGGCGGTGCAGGGCCACCTGTACTCGGTCAGTGAGCATGATGAAGACCGGATGGCCGAGCAGCCAACGTCTCGCAAATACGACACCTATGTGCGGGCGCTGTCATTCGAAGGTGACTGGCGCAAGGTGGAGCAGGGGCGTTTg                                                                                                    <  5:235741/149‑1 (MQ=255)                                           gTGGCTGGCGGTGCAGGGCCACCTGTACTCGGTCAGTGAGCATGATGAAGACCGGATGGCCGAGCAGCCAACGTCTCGCAAATACGACACCTATGTGCGGGCGCTGTCATTCGAAGGTGACTGGCGCAAGGTGGAGCAGGGGCGTTTGc                                                                                                   >  8:94881/1‑149 (MQ=255)                                                tggcGGTGCAGGGCCACCTGTACTCGGTCAGTGAGCATGATGAAGACCGGATGGCCGAGCAGCCAACGTCTCGCAAATACGACACCTATGTGCGGGCGCTGTCATTCGAAGGTGACTGGCGCAAGGTGGAGCAGGGGCGTTTGCGTATc                                                                                              >  8:122101/1‑149 (MQ=255)                                                  gcGGTGCAGGGCCACCTGTACTCGGTCAGTGAGCATGATGAAGACCGGATGGCCGAGCAGCCAACGTCTGGCAAATACGACACCTATGTGCGGGCGCTGTCATTCGAATGTGACTGGCGCAAGGTGGAGCAGGGGCGTTTGCGTATc                                                                                              <  1:85170/147‑1 (MQ=255)                                                     gTGCAGGGCCACCTGTACTCGGTCAGTGAGCATGATGAAGACCGGATGGCCGAGCAGCCAACGTCTCGCAAATATGACACCTATGTGCGGGCGCTGTCATTCGAAGGTGACTGGCGCAAGGTGGAGCAGGGGCGTTTGc                                                                                                   <  4:142888/139‑1 (MQ=255)                                                     gTGCAGGGCCACCTGTACTCGGTCAGTGAGCATGATGAAGACCGGATGGCCGAGCAGCCAACGTCTCGCAAATATGACACCTATGTGCGGGCGCTGTCATTCGAAGGTGACTGGCGCAAGGTGGAGCAGGGGCGTTTGc                                                                                                   >  3:142888/1‑139 (MQ=255)                                                       gCAGGGCCACCTGTACTCGGTCAGTGAGCATGATGAAGACCGGATGGCCGAGCAGCCAACGTCTCGCAAATACGACACCTATGTGCGGGCGCTGTCATTCGAAGGTGACTGGCGCAAGGTGGAGCAGGGGCGTTTGCGTATCGTCAGTg                                                                                       >  3:54075/1‑149 (MQ=255)                                                                      cTCGGTCAGTGAGCATGATGAAGACCGGATGGCCGAGCAGCCAACGTCTCGCAAATACGACACCTATGTGCGGGCGCTGTCATTCGAAGGTGACTGGCGCAAGGTGGAGCAGGGGCGTTTGCGTATCGTCAGTGCCGGGGTCCATCAgg                                                                        <  1:89028/149‑1 (MQ=255)                                                                       tCGGTCAGTGAGCATGATGAAGACCGGATGGCCGAGCAGCCAACGTCTCGCAAATACGACACCTATGTGCGGGCGCTGTCATTCGAAGGTGACTGGCGCAAGGTGGAGCAGGGGCGTTTGCGTATCGTCAGTGCCGGGGTCCATCAggg                                                                       >  4:126036/1‑149 (MQ=255)                                                                         ggTCAGTGAGCATGATGAAGACCGGATGGCCGAGCAGCCAACGTCTCGCAAATACGACACCTATGTGCGGGCGCTGTCATTCGAAGGTGACTGGCGCAAGGTGGAGCAGGGGCGTTTGCGTATCGTCAGTGCCGGGGTCCATCAGGGTc                                                                     >  8:249069/1‑149 (MQ=255)                                                                         ggTCAGTGAGCATGATGAAGACCGGATGGCCGAGCAGCCAACGTCTCGCAAATACGACACCTATGTGCGGGCGCTGTCATTCGAAGGTGACTGGCGCAAGGTGGAGCAGGGGCGTTTGCGTATCGTCAGTGCCGGGGTCCATCAGGGTc                                                                     >  1:75678/1‑149 (MQ=255)                                                                            cAGTGAGCATGATGAAGACCGGATGGCCGAGCAGCCAACGTCTCGCAAATACGACACCTATGTGCGGGCGCTGTCATTCGAAGGTGACTGGCGCAAGGTGGAGCAGGGGCGTTTGCGTATCGTCAGTGCCGGGGTCCATCAGGGTCTg                                                                   >  6:201197/1‑148 (MQ=255)                                                                             agagAGCATGATGAAGAACGGATGGCCGAGCAGCCAACGTCTAGCAAATACGACACCTATGTGCGGACCCTGTCATTCGAAGGAGACTGCCGCAAGGTGGAGCAGGGGCGATTGCGTATCGTCAGTTCCGGCGTCCATCAGGGTCTga                                                                  >  7:191859/4‑147 (MQ=255)                                                                             aGTGAGCATGATGAAGACCGGATGGCCGAGCAGCCAACGTCTCGCAAATACGACACCTATGTGCGGACGCTGTCATTCGAAGGTGACTGGCGCAAGGTGGAGCAGGGGCGTTTGCGTATCGTCAGTGCCGGGGTCCATCAGGGTCTGGa                                                                 >  7:191855/1‑149 (MQ=255)                                                                                 aGCATGATGAAGACCGGATGGCCGAGCAGCCAACGTCTCGCAAATACGACACCTATGTGCGGGCGCTGTCATTCGAAGGTGACTGGCGCAAGGTGGAGCAGGGGCGTTTGCGTATCGTCAGTGCCGGGGTCCATCAGGGTCTGGACTAc                                                             <  1:155907/149‑1 (MQ=255)                                                                                     tgatgaAGACCGGATGGCCGAGCAGCCAACGTCTCGCAAATACGACACCTATGTGCGGGCGCTGTCATTCGAAGGTGACTGGCGCAAGGTGGAGCAGGGGCGTTTGCGTATCGTCAGTGCCGGGGTCCATCAGGGTCTGGACTACCTgg                                                         <  2:75678/149‑1 (MQ=255)                                                                                        tgaAGACCGGATGGCCGAGCAGCCAACGTCTCGCAAATACGACACCTATGTGCGGGCGCTGTCATTCGAAGGTGACTGGCGCAAGGTGGAGCAGGGGCGTTTGCGTATCGTCAGTGCCGGGGTCCATCAGGGTCTGGACTACCTggggg                                                      <  3:129430/149‑1 (MQ=255)                                                                                             aCCGGATGGCCGAGCAGCCAACGTCTCGCAAATACGACACCTATGTGCGGGCGCTGTCATTCGAAGGTGACTGGCGCAAGGTGGAGCAGGGGCGTTTGCGTATCGTCAGTGCCGGGGTCCATCAGGGTCTGGACTACCTGGGGGCACGc                                                 <  7:249069/149‑1 (MQ=255)                                                                                               cGGATGGCCGAGCAGCCAACGTCTCGCAAATACGACACCTATGTGGGGGCGCTGTCATTCGaa                                                                                                                                     >  1:189851/1‑63 (MQ=255)                                                                                               cGGATGGCCGAGCAGCCAACGTCTCGCAAATACGACACCTATGTGCGGGCGCTGTCATTCGaa                                                                                                                                     <  2:189851/63‑1 (MQ=255)                                                                                                 gATGGCCGAGCAGCCAACGTCTCGCAAATACGACACCTATGTGCGGGCGCTGTCATTCGAAGGTGACTGGCGCAAGGTGGAGCAGGGGCGTTTGCGTATCGTCAGTGCCGGGGTCCATCAGGGTCTGGACTACCTGGGGGCACGCAGc                                              >  1:46103/1‑148 (MQ=255)                                                                                                  aTGGCCGAGCAGCCAACGTCTCGCAAATACGACACCTATGTGCGGGCGCTGTCATTCGAAGGTGACTGGCGCAAGGTGGAGCAGGGGCGTTTGCGTATCGTCAGTGCCGGGGTCCATCAGGGTCTGGACTACCTGGGGGCACGCAGCGa                                            <  2:228540/149‑1 (MQ=255)                                                                                                        gAGCAGCCAACGTCTCGCAAATACGACACCTATGTGCGGGCGCTGTCATTCGAAGGTGACTGGCGCAAGGTGGAGCAGGGGCGTTTGCGTATCGTCAGTGCCGGGGTCCATCAGGGTCTGGACTACCTGGGGGCACGCAGCGATGCTGa                                      <  6:162286/149‑1 (MQ=255)                                                                                                        gAGCAGCCAACGTCTCGCAAATACGACACCTATGTGCGGGCGCTGTCATTCGAAGGTGACTGGCGCAAGGTGGAGCAGGGGCGTTTGCGTATCGTCAGTGCCGGGGTCCATCAGGGTCTGGACTACCTGGGGGCACGCAGCGATGCTGa                                      <  7:188648/149‑1 (MQ=255)                                                                                                         agcagcCAACGTCTCGCAAATACGACACCTATGTGCGGGCGCTGTCATTCGAAGGTGACTGGCGCAAGGTGGAGCAGGGGCGTTTGCGTATCGTCAGTGCCGGGGTCCATCAGGGTCTGGACTACCTGGGGGCACGCAGCGATGCTGAc                                     >  2:24150/1‑149 (MQ=255)                                                                                                               cAACGTCTCGCAAATACGACACCTATGTGCGGGCGCTGTCATTCGAAGGTGACTGGCGCAAGGTGGAGCAGGGGCGTTTGCGTATCGTCAGTGCCGGGGTCCATCAGGGTCTGGACTACCTGGGGGCACGCAGCGATGCTGACTACGAt                               <  4:54075/149‑1 (MQ=255)                                                                                                                   gTCTCGCAAATACGACACCTATGTGCGGGCGCTGTCATTCGAAGGTGACTGGCGCAAGGTGGAGCAGGGGCGTTTGCGTATCGTCAGTGCCGGGGTCCATCAGCGTCTGGACTACCTGGGGGCACGCAGCGATGCTGACTACGAtctc                            >  3:123167/1‑148 (MQ=255)                                                                                                                           aaTACGACACCTATGTGCGGGCGCTGTCATTCGAAGGTGACTGGCGCAAGGTGGAGCAGGGGCGTTTGCGTATCGTCAGTGCCGGGGTCCATCAGGGTCTGGACTACCTGGGGGCACGCAGCGATGCTGACTACGATCTCGACTTCCTg                   <  2:46103/149‑1 (MQ=255)                                                                                                                             tACGACACCTATGTGCGGGCGCTGTCATTCGAAGGTGACTGGCGCAAGGTGGAGCAGGGGCGTTTGCGTATCGTCAGTGCCGGGGTCCATCAGGGTCTGGACTACCTGGGGGCACGCAGCGATGCTGACTACGATCTCGACTTCCTgcg                 <  5:201197/149‑1 (MQ=255)                                                                                                                                     cTATGTGCGGGCGCTGTCATTCGAAGGTGACTGGCGCAAGGTGGAGCAGGGGCGTTTGCGTATCGTCAGTGCCGGGGTCCATCAGGGTCTGGACTACCTgg                                                         <  2:141763/101‑1 (MQ=255)                                                                                                                                     cTATGTGCGGGCGCTGTCATTCGAAGGTGACTGGCGCAAGGTGGAGCAGGGGCGTTTGCGTATCGTCAGTGCCGGGGTCCATCAGGGTCTGGACTACCTgg                                                         >  1:141763/1‑101 (MQ=255)                                                                                                                                     cTATGTGCGGGCGCTGTCATTCGAAGGTGACTGGCGCAAGGTGGAGCAGGGGCGTTTGCGTATCGTCAGTGCCGGGGTCCATCAGGGTCTGGACTACCTGGGGGCACGCAGCGATGCTGACTACGATCTCGACTTCCTGCGCGTGCGGc         >  2:83887/1‑149 (MQ=255)                                                                                                                                      tATGTGCGGGCGCTGTCATTCGAAGGTGACTGGCGCAAGGTGGAGCAGGGGCGTTTGCGTATCGTCAGTGCCGGGGTCCATCAGGGTCTGGACTACCTGGGGGCACGCAGCGATGCTGACTACGATCTCGACTTc                      >  2:260474/1‑135 (MQ=255)                                                                                                                                      tATGTGCGGGCGCTGTCATTCGAAGGTGACTGGCGCAAGGTGGAGCAGGGGCGTTTGCGTATCGTCAGTGCCGGGGTCCATCAGGGTCTGGACTACCTGGGGGCACGCAGCGATGCTGACTACGATCTCGACTTc                      <  1:260474/135‑1 (MQ=255)                                                                                                                                           gCGGGCGCTGTCATTCGAAGGTGACTGGCGCAAGGTGGAGCAGGGGCGTTTGCGTATCGTCAGTGCCGGGGTCCATCAGGGTCTGGACTACCTGGGGGCACGCAGCGATGCTGACTACGATCTCGACTTCCTGCGCTTGCGGCTTTCCg   >  2:190378/1‑149 (MQ=255)                                                                                                                                            cGGACGCTGTCATTCGACGGTGACTGGCGCAAGGTGGAGCAGGGGCGTTTGCGTATCGTCAGTGCCGGTGTCCATCAGGGTCTGGACTACCTGGGGTCACGCAGCGATGCTGACTACGATCTCGACTTCCTGCGCTTGCGGCTTTCCgg  <  8:191859/149‑1 (MQ=255)                                                                                                                                            cGGACGCTGTCATTCGAAGGTGACTGGCGCAAGGTGGAGCAGGGGCGTTTGCGTATCGTCAGTGCCGGGGTCCATCAGGGTCTGGACTACCTGGGGGCACGCAGCGATGCTGACTACGATCTCGACTTCCTGCGCTTGCGGCTTTCCgg  <  8:191855/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                  TGCGATCGGCCTGCGTCAGCCGGTGGTTGTCAGGCCGAATGAGTGGCTGGCGGTGCAGGGCCACCTGTACTCGGTCAGTGAGCATGATGAAGACCGGATGGCCGAGCAGCCAACGTCTCGCAAATACGACACCTATGTGCGGGCGCTGTCATTCGAAGGTGACTGGCGCAAGGTGGAGCAGGGGCGTTTGCGTATCGTCAGTGCCGGGGTCCATCAGGGTCTGGACTACCTGGGGGCACGCAGCGATGCTGACTACGATCTCGACTTCCTGCGCTTGCGGCTTTCCGG  >  NC_002947/671956‑672243

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