Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,918,174	C→G	10.9%	L191L (CTC→CTG)	PP_4328 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,918,174	0	C	G	10.9%	160.6 / 9.7	55	L191L (CTC→CTG)	PP_4328	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (21/28);  new base G (3/3);  total (24/31)
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.20e-01

CGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCC  >  NC_002947/4918029‑4918315                                                                                                                                                   |                                                                                                                                              cGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGc                                                                                                                                            <  6:112848/149‑1 (MQ=255) cGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGc                                                                                                                                            >  5:112848/1‑149 (MQ=255)  ggCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCg                                                                                                                                           >  1:270888/1‑149 (MQ=255)      gtgtACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCggcag                                                                                                                                       >  6:162422/1‑149 (MQ=255)       tgtACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCggcagg                                                                                                                                      <  8:271017/149‑1 (MQ=255)         tACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCCCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTGAGCggcacga                                                                                                                                     <  1:278065/148‑4 (MQ=255)         tACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGGCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATAGCGCTTCGTGCCGCTGAGCggcacgaa                                                                                                                                    <  5:228958/149‑5 (MQ=255)           ccGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGcaggcagg                                                                                                                                  <  1:300304/149‑1 (MQ=255)            cGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGa                                                                                                                                 <  5:114793/149‑1 (MQ=255)             gggCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGTCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGAc                                                                                                                                <  1:46661/149‑1 (MQ=255)              ggCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGa                                                                                                                                 <  6:219394/147‑1 (MQ=255)              ggCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGa                                                                                                                                 >  5:219394/1‑147 (MQ=255)                     gTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTc                                                                                                                          <  7:299004/147‑1 (MQ=255)                     gTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAg                                                                                                                        >  7:110454/1‑149 (MQ=255)                          ccTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCACCTCAGCATTg                                                                                                                   >  5:68267/1‑149 (MQ=255)                            tGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGcc                                                                                                                 >  4:346796/1‑149 (MQ=255)                               tCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCAATGCCCaa                                                                                                              >  5:158956/1‑149 (MQ=255)                                       aaCCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCGCTCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAgctgct                                                                                                      <  6:155130/149‑1 (MQ=255)                                         ccAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGa                                                                                                    <  5:102018/149‑1 (MQ=255)                                             gCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTc                                                                                                                          <  3:317118/123‑1 (MQ=255)                                             gCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTc                                                                                                                          >  4:317118/1‑123 (MQ=255)                                             gCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGAcc                                                                                                  >  2:173573/1‑147 (MQ=255)                                             gCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCa                                                                                                <  8:110454/149‑1 (MQ=255)                                                            aGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCaa                                                                                   <  5:162422/147‑1 (MQ=255)                                                            aGCATCGACAGCCCGCGCCCGCTGGGCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCaa                                                                                   <  6:299757/147‑1 (MQ=255)                                                            aGCATCGACAGCCCGCGCCAGCTGGTCATTGGAAGGCTTCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCTCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGc                                                                                 <  7:134854/149‑1 (MQ=255)                                                              cATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCt                                                                               >  3:193542/1‑149 (MQ=255)                                                               aTCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTc                                                                              >  8:102636/1‑149 (MQ=255)                                                                   aCAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTa                                                                          >  4:89220/1‑149 (MQ=255)                                                                        ccgcgcccgcT‑GTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGc                                                                    >  6:179217/1‑149 (MQ=255)                                                                              ccgcTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCt                                                               >  5:291556/1‑149 (MQ=255)                                                                                    gTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGc                                                         >  4:126881/1‑149 (MQ=255)                                                                                          ggCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGc                                                    <  6:158956/148‑1 (MQ=255)                                                                                           gCAGCCTGCTCAGCGCCCTGGTGCAGGGTGAGTCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCGACGCCCTGc                                                    >  1:219365/1‑148 (MQ=255)                                                                                           gCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCGACGCCCTGCaa                                                  >  7:218316/1‑149 (MQ=255)                                                                                                   cTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCAc                                          <  4:193542/149‑1 (MQ=255)                                                                                                       gcgcCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCGACGCCCTGCAACAGCTCACCcgc                                      <  2:219365/149‑1 (MQ=255)                                                                                                       gcgcCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCGACGCCCTGCAACAGCTCACCcgc                                      <  8:218316/149‑1 (MQ=255)                                                                                                            cTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGAccc                                 <  7:102636/149‑1 (MQ=255)                                                                                                                     ggTGA‑TCAGTCATTGGGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGa                        <  6:16572/149‑1 (MQ=255)                                                                                                                       tGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATg                      <  3:126881/149‑1 (MQ=255)                                                                                                                         a‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATgg                     <  5:42318/148‑1 (MQ=255)                                                                                                                           tCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGc                    <  3:89220/148‑1 (MQ=255)                                                                                                                             aGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGAc                 <  6:291556/149‑1 (MQ=255)                                                                                                                                 aTTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTCCAACAGCTCACCCGCCACCCAGGCGTCGATGCCGACCTgc             >  6:12942/1‑149 (MQ=255)                                                                                                                                    gcgcTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTgcgc           <  2:270888/148‑1 (MQ=255)                                                                                                                                    gcgcTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTgcgc           <  5:12942/148‑1 (MQ=255)                                                                                                                                         gcctGCCGCTGAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGa                        >  8:274115/4‑130 (MQ=255)                                                                                                                                         gcctGCCGCTGAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGa                        <  7:274115/127‑1 (MQ=255)                                                                                                                                         gcctGCCGCTGAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAgc     >  1:163317/4‑149 (MQ=255)                                                                                                                                          cctGCCGCTGAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAgcg    >  6:17814/3‑149 (MQ=255)                                                                                                                                           gTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAgcg    <  3:346796/148‑1 (MQ=255)                                                                                                                                            tGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGcc  >  7:338914/1‑149 (MQ=255)                                                                                                                                             gCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGAc                 <  2:207736/133‑1 (MQ=255)                                                                                                                                             gCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGAc                 >  1:207736/1‑133 (MQ=255)                                                                                                                                                   |                                                                                                                                              CGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCCCGCTGGTCATTGGCAGCCTGCTCAGCGCCCTGGTGCAGGGTGA‑TCAGTCATTGCGCTTCGTGCCGCTCAGCGGCAGGCAGGACCAACTCAGCATTGCCCAACAGCTGCTGACCCAGCAAAACCGCCAAGCCTCGCTACCCGGCCTGCTCAACGCCCTGCAACAGCTCACCCGCGACCCAGGCGTCGATGGCGACCTGCGCGCCAGCGCC  >  NC_002947/4918029‑4918315

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