Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,667,919	C→T	14.5%	A307A (GCC→GCT)	ffh →	signal recognition particle GTPase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,667,919	0	C	T	14.5%	114.5 / 8.9	48	A307A (GCC→GCT)	ffh	signal recognition particle GTPase subunit
Reads supporting (aligned to +/- strand):  ref base C (19/22);  new base T (3/4);  total (22/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

GCGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATGGAC  >  NC_002947/1667773‑1668054                                                                                                                                                   |                                                                                                                                         gcgcTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGcc                                                                                                                                          >  6:327437/1‑147 (MQ=255)          gTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCAt                                                                                                                               <  5:327437/149‑1 (MQ=255)           tACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATc                                                                                                                              <  3:243469/149‑1 (MQ=255)              gTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAg                                                                                                                                       >  6:85102/1‑137 (MQ=255)              gTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAg                                                                                                                                       <  5:85102/137‑1 (MQ=255)               tGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGAc                                                                                                                           >  2:144984/1‑148 (MQ=255)                   aTCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGAcaaggc                                                                                                                      <  7:116026/149‑1 (MQ=255)                    tCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGAcaaggc                                                                                                                      >  1:181153/1‑148 (MQ=255)                          gCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGAcaag                                                                                                                        <  1:284855/140‑1 (MQ=255)                          gCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGAcaag                                                                                                                        >  2:284855/1‑140 (MQ=255)                             taCCGATGAACTTGATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGAACGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGa                                                                                                             >  2:79029/2‑148 (MQ=255)                                     aaGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTgg                                                                                                     <  2:181153/148‑1 (MQ=255)                                       gTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCGGCCTGCTCGATCAGGCTGAGc                                                                                                                                      <  4:82474/113‑1 (MQ=255)                                           aTCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGaa                                                                                              >  4:230119/1‑149 (MQ=255)                                                  tGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTgaag                                                                                       <  3:160106/149‑1 (MQ=255)                                                  tGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTg                                                                                          >  5:286566/1‑146 (MQ=255)                                                    ggTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTgaagaa                                                                                     <  1:144984/149‑1 (MQ=255)                                                        agaagaCCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGaag‑ggc                                                                                <  2:201302/149‑1 (MQ=255)                                                         gaagaCCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCGGCCTGCTCGATCAGGCTGAGCacg                                                                                                                                   >  3:250016/1‑96 (MQ=38)                                                         gaagaCCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCGGCCTGCTCGATCAGGCTGAGCacg                                                                                                                                   <  4:250016/98‑3 (MQ=38)                                                               cGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGaag‑ggcaag                                                                             <  6:355184/145‑1 (MQ=255)                                                               cGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGaag‑ggcaag                                                                             >  5:355184/1‑145 (MQ=255)                                                                      cTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCt                                                                  >  2:277829/1‑149 (MQ=255)                                                                       tCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGaaa                                                                                             <  8:156243/122‑1 (MQ=255)                                                                       tCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGaaa                                                                                             >  7:156243/1‑122 (MQ=255)                                                                                 ccACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTc                                                        <  1:32247/148‑1 (MQ=255)                                                                                 ccACCCCGACCGCCTCGCTTCGCGCATCCGCGGCATGGGCGACGTGCTCCGCCTGATCGAGCAGGACGAGCAGACCATCGACAAGTCCAAGGCCGACAACCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTc                                                        <  2:43855/148‑1 (MQ=255)                                                                                   aCCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCcgc                                                     <  6:286566/149‑1 (MQ=255)                                                                                    ccccGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGc                                                                                                    <  2:332225/102‑1 (MQ=255)                                                                                       cGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCCGCGAc                                                  >  7:278977/1‑148 (MQ=255)                                                                                       cGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCCGCGAc                                                  <  8:278977/148‑1 (MQ=255)                                                                                                     cgcgCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGTGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCaaa                                     <  1:79029/148‑1 (MQ=255)                                                                                                      gcgcATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGTCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATga                                  <  1:277829/149‑1 (MQ=255)                                                                                                      gcgcATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATg                                   >  2:179769/1‑148 (MQ=255)                                                                                                              tCGGCGTGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGaa                                                              >  8:267520/1‑113 (MQ=255)                                                                                                              tCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCTAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGAc                                                           >  1:109592/1‑116 (MQ=255)                                                                                                              tCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGAc                                                           <  2:109592/116‑1 (MQ=255)                                                                                                                       gCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTc                 >  8:78018/1‑149 (MQ=255)                                                                                                                           tctGCTCGGCCTGCTCGATCAGGCTGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTcggcg             >  7:106493/3‑149 (MQ=255)                                                                                                                               cTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACTAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCt         >  6:301176/1‑149 (MQ=255)                                                                                                                                tcggCCTGCTCGATCAGGCTGAGCAGTCCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCt                  <  2:143332/136‑1 (MQ=255)                                                                                                                                tcggCCTGCTCGATCAGGCTGAGCAGGCCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCt                  >  1:143332/4‑139 (MQ=255)                                                                                                                                tCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCggcggc           >  6:64244/1‑146 (MQ=255)                                                                                                                                tCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCggcggc           <  5:64244/146‑1 (MQ=255)                                                                                                                                   gCCTGCTCGATCAGGCTGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATg     <  8:106493/149‑1 (MQ=255)                                                                                                                                   gCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATg     <  7:78018/149‑1 (MQ=255)                                                                                                                                    ccTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAg                                               <  2:4756/106‑1 (MQ=255)                                                                                                                                    ccTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAg                                               >  1:4756/1‑106 (MQ=255)                                                                                                                                       gATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATGGAc  >  6:41810/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                         GCGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATGGAC  >  NC_002947/1667773‑1668054

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