Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,843,877	T→G	9.3%	I242S (ATC→AGC)	PP_4259 →	iron‑sulfur cluster‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,843,877	0	T	G	9.3%	126.7 / 3.9	43	I242S (ATC→AGC)	PP_4259	iron‑sulfur cluster‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (22/17);  new base G (2/2);  total (24/19)
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.81e-01

AGGCCGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTT‑CGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTGACGGCCTGCA  >  NC_002947/4843733‑4844024                                                                                                                                                  |                                                                                                                                                    aGGCCGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTT‑CGACAAGGACACGCTGATCgtg                                                                                                                                                 >  4:67185/1‑149 (MQ=255)     cGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTT‑CGACAAGGACACGCTGATCGTGTCCt                                                                                                                                             <  8:282893/149‑1 (MQ=255)            tGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTT‑CGACAAGGACACGCTGATCgtg                                                                                                                                                 >  4:10959/1‑138 (MQ=255)            tGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTT‑CGACAAGGACACGCTGATCgtg                                                                                                                                                 <  3:10959/138‑1 (MQ=255)                ccTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTT‑CGACAAGGACACGCTGATCGTGTCCTACGACCCACGc                                                                                                                                  >  1:138087/1‑149 (MQ=255)                ccTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTT‑CGACAAGGACACGCTGATCGTGTCCTACGACCCACGc                                                                                                                                  >  6:11589/1‑149 (MQ=255)                 cTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGACGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCACCTATGCCCGCTTCCAGAGCGTGATGTT‑CGACAAGGACACGCTGATc                                                                                                                                                    >  7:538/1‑130 (MQ=255)                    cTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTT‑CGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTg                                                                                                                              <  4:26651/149‑1 (MQ=255)                                       aCCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTT‑CGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCcgc                                                                                                            <  6:181569/148‑1 (MQ=255)                                                         aaCGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTT‑CGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTa                                                                                         >  4:1414/1‑149 (MQ=255)                                                          aCGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTT‑CGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTa                                                                                         >  6:89784/1‑148 (MQ=255)                                                            gCCGGCTGGCTACACGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTT‑CGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCcgc                                                                                                            <  3:216664/127‑1 (MQ=255)                                                            gCCGGCTGGCTACACGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTT‑CGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCcgc                                                                                                            >  4:216664/1‑127 (MQ=255)                                                             ccGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCCGACAAGGACACGATCAGCGTGTCCttgt                                                                                                                                          <  1:21230/97‑4 (MQ=25)                                                             ccGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTT‑CGACAAGGACACGATCAGCGTGTCCttgt                                                                                                                                          >  2:21230/1‑93 (MQ=25)                                                                         gcgAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTT‑CGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTg                                                                         <  3:153216/149‑1 (MQ=255)                                                                          cgAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTT‑CGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTgg                                                                        <  5:89784/149‑1 (MQ=255)                                                                           gAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTT‑CGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTggg                                                                       >  4:219687/1‑149 (MQ=255)                                                                                        tgtACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTT‑CGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTg                                                                         >  5:261942/1‑134 (MQ=255)                                                                                        tgtACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTT‑CGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTg                                                                         <  6:261942/134‑1 (MQ=255)                                                                                          tACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTT‑CGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTg                                                        <  3:1414/149‑1 (MQ=255)                                                                                             aTGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTT‑CGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCAc                                                     >  5:52729/1‑149 (MQ=255)                                                                                                  ccccTATGCCCGCTTCCAGAGCGTGATGTT‑CGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCGTCGACTGCACCAtgt                                                >  8:142627/1‑149 (MQ=255)                                                                                                  ccccTATGCCCGCTTCCAGAGCGTGATGTT‑CGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCAtgt                                                <  3:219687/149‑1 (MQ=255)                                                                                                     cTATGCCCGCTTCCAGAGCGTGATGTT‑CGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATAc                                                                                               <  5:6693/99‑1 (MQ=255)                                                                                                     cTATGCCCGCTTCCAGAGCGTGATGTT‑CGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATAc                                                                                               >  6:6693/1‑99 (MQ=255)                                                                                                               ttCCAGAGCGTGATGTT‑CGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACt                                                         >  8:121218/1‑127 (MQ=255)                                                                                                               ttCCAGAGCGTGATGTT‑CGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACt                                                         <  7:121218/127‑1 (MQ=255)                                                                                                                            tGTT‑CGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATc                      <  6:52729/149‑1 (MQ=255)                                                                                                                                  gACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGt                                      <  7:98074/128‑1 (MQ=255)                                                                                                                                  gACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGt                                      >  8:98074/1‑128 (MQ=255)                                                                                                                                  gACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACAt                 >  5:221642/1‑149 (MQ=255)                                                                                                                                    cAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATcc               <  3:67185/149‑1 (MQ=255)                                                                                                                                     taGGACACGATCAGCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAg                                        >  7:272946/2‑123 (MQ=255)                                                                                                                                     taGGACACGATCAGCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAg                                        <  8:272946/122‑1 (MQ=255)                                                                                                                                       ggACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCAGCATGTGCGTGCAGGTCTGCCCGACCGGCATCAACATCCGt             >  5:308480/1‑148 (MQ=255)                                                                                                                                       ggACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGt             >  7:135184/1‑148 (MQ=255)                                                                                                                                       ggACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTg            >  6:187923/1‑149 (MQ=255)                                                                                                                                           acGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTGACgg        >  1:170374/1‑149 (MQ=255)                                                                                                                                             gCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGAcc                            >  6:259557/1‑127 (MQ=255)                                                                                                                                             gCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGAcc                            <  5:259557/127‑1 (MQ=255)                                                                                                                                             gCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATc                      <  6:315376/133‑1 (MQ=255)                                                                                                                                             gCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATc                      >  5:315376/1‑133 (MQ=255)                                                                                                                                                  tCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCTTCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTGACGGCCTGCa  <  7:142627/148‑1 (MQ=255)                                                                                                                                                  tCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTGACGGCCTg    >  7:241135/1‑146 (MQ=255)                                                                                                                                                  |                                                                                                                                                    AGGCCGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTT‑CGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTGACGGCCTGCA  >  NC_002947/4843733‑4844024

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