Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,050,935	T→G	6.7%	F256C (TTC→TGC)	PP_2663 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,050,935	0	T	G	6.7%	173.4 / 2.4	56	F256C (TTC→TGC)	PP_2663	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (24/28);  new base G (2/2);  total (26/30)
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

ATCACATCCTGCCGCGCCAGGAAAAGCTGGTGGTAACCGGTGCCGACAGCGCCCAGCGCCGGGTTTTCGAACTGAATGCCGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCATGATCAG‑TACTACGTGCGGGCGATCCAGCAAATGCATCCTGACCTCAGCTTGAGCTTCTACTGCGCGGTGGAGAACGGCATCGTGCTCACGGCCATGACCCCGGGCCCGCTG  >  NC_002947/3050788‑3051077                                                                                                                                                    |                                                                                                                                                aTCACATCCTGCCGCGCCAGGAAAAGCTGGTGGTAACCGGTGCCGACAGCGCCCAGCGCCGGGTTTTCGAACTGAATGCCGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTc                                                                                                                                                <  4:167953/149‑1 (MQ=255)    acaTCCTGCCGCGCCAGGAAAAGCTGGTGGTAACCGGTGCCGACAGCGCCCAGCGCCGGGTTTTCGAACTGAATGCCGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTcgcc                                                                                                                                             >  4:266425/1‑149 (MQ=255)      aTCCTGCCGCGCCAGGAAAAGCTGGTGGTAACCGGTGCCGACAGCGCCCAGCGCCGGGTTTTCGAACTGAATGCCGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTcgccgc                                                                                                                                           >  2:308378/1‑149 (MQ=255)      aTCCTGCCGCGCCAGGAAAAGCTGGTGGTAACCGGTGCCGACAGCGCCCAGCGCCGGGTTTTCGAACTGAATGCCGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTcgccgc                                                                                                                                           >  6:360137/1‑149 (MQ=255)                  cAGGAAAAGCTGGTGGTAACCGGTGCCGACAGCGCCCAGCGCCGGGTTTTCGAACTGAATGCCGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGc                                                                                                                               >  7:186624/1‑149 (MQ=255)                      aaaaGCTGGTGGTAACCGGTGCCGACAGCGCCCAGCGCCGGGTTTTCGAACTGAATGCCGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTg                                                                                                                           <  4:110066/149‑1 (MQ=255)                      aaaaGCTGGTGGTAACCGGTGCCGACAGCGCCCAGCGCCGGGTTTTCGAACTGAATGCCGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTg                                                                                                                           <  1:77013/149‑1 (MQ=255)                          gCTGGTGGTAACCGGTGCCGACAGCGCCCAGCGCCGGGTTTTCGAACTGAATGCCGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCa                                                                                                                       >  2:253580/1‑149 (MQ=255)                              gtggtAACCGGTGCCGACAGCGCCCAGCGCCGGGTTTTCGATCTGAATGCCGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCAtt                                                                                                                     <  1:116208/147‑1 (MQ=255)                              gtggtAACCGGTGCCGACAGCGCCCAGCGCCGGGTTTTCGAACTGAATGCCGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCa                                                                                                                   >  8:61805/1‑149 (MQ=255)                              gtggtAACCGGTGCCGACAGCGCCCAGCGCCGGGTTTTCGAACTGAATGCCGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCACGGCCATCGGGTGGGCGGCGAACACCCGGTGATCGAGCGCagc                                                                                                                     <  7:61870/147‑3 (MQ=25)                                 gtAACCGGTGCCGACAGCGCCCAGCGCCGGGTTTTCGAACTGAATGCCGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCATg                                                                                                                 <  3:127816/148‑1 (MQ=255)                                    aCCGGTGCCGACAGCGCCCAGCGCCGGGTTTTCGAACTGAATGCCGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCATGATCa                                                                                                             <  3:200126/149‑1 (MQ=255)                                      cGGTGCCGACAGCGCCCAGCGCCGGGTTTTCGAACTGAATGCCGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCATGATCAg                                                                                                            >  6:176915/1‑148 (MQ=255)                                      cGGTGCCGACAGCGCCCAGCGCCGGGTTTTCGAACTGAATGCCGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCATGATCAg                                                                                                            >  5:243102/1‑148 (MQ=255)                                        gTGCCGACAGCGCCCAGCGCCGGGTTTTCGAACTGAATGCCGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTc                                                                                                                                                <  7:134810/110‑1 (MQ=255)                                        gTGCCGACAGCGCCCAGCGCCGGGTTTTCGAACTGAATGCCGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTc                                                                                                                                                >  8:134810/1‑110 (MQ=255)                                        gTGCCGACAGCGCCCAGCGCCGGGTTTTCGAACTGAATGCCGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCATGATCAG‑ta                                                                                                         >  6:141536/1‑148 (MQ=255)                                         tGCCGACAGCGCCCAGCGCCGGGTTTTCGAACTGAATGCCGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCATGATCAG‑tac                                                                                                        >  5:257439/1‑148 (MQ=255)                                            cGACAGCGCCCAGCGCCGGGTTTTCGAACTGAATGCCGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCATGATCAG‑TACTACg                                                                                                    >  3:69004/1‑149 (MQ=255)                                            cGACAGCGCCCAGCGCCGGGTTTTCGAACTGAATGCCGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCATGATCAG‑TACTACg                                                                                                    >  6:209914/1‑149 (MQ=255)                                              aCAGCGCCCAGCGCCGGGTTTTCGAACTGAATGCCGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCATGATCAG‑TACTACGTg                                                                                                  >  5:247166/1‑149 (MQ=255)                                               cagcgcccagcgccGGGTTTTCGAACTGAATGCCGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCATGATCAG‑TACTACGTGc                                                                                                 <  6:257439/149‑1 (MQ=255)                                                 gcgcccagcgccGGGTTTTCGAACTGTATGCCGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCGCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGGGCGCATTCATGATCAG‑TACTACGTGctg                                                                                               <  1:243936/149‑3 (MQ=255)                                                  cgcccagcgccGGGTTTTCGAACTGAATGCCGAGCCGGCGGCCGAGGAATATGCCCGCCACATCTGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCATGATCAG‑TACTACGTGCggg                                                                                              <  6:243102/149‑1 (MQ=255)                                                     ccagcgccGGGTTTTCGAACTGAATGCCGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCATGATCAG‑TACTACGTGCGGGCGa                                                                                           >  5:342153/1‑149 (MQ=255)                                                        gcgccGGGTTTTCGAACTGAATGCCGAGCCGGCGGGCGAGGAATATGCCCGCCACATCGGTGAGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCATGATCAG‑TACTACGTGCGGGCGATcc                                                                                        <  1:327003/149‑1 (MQ=255)                                                           ccGGGTTTTCGAACTGAATGCCGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCATCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCATGATCAG‑TACTACGTGCGGGCGATCCAGc                                                                                     <  1:308378/149‑1 (MQ=255)                                                           ccGGGTTTTCGAACTGAATGCCGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCATGATCAG‑TACTACGTGCGGGCGATCCAGc                                                                                     <  5:360137/149‑1 (MQ=255)                                                              ggTTTTCGAACTGAATGCCGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCATGATCAG‑TACTACGTGCGGGCGATCCAGCaaa                                                                                  <  5:209914/149‑1 (MQ=255)                                                                    cGAACTGAATGCCGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCATGATCAG‑TACTACGTGCGGGCGATCCAGCAAATGCAt                                                                             <  1:253580/148‑1 (MQ=255)                                                                         tGAATGCCGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCATGATCAG‑TACTACGTGCGGGCGATCCAGCAAATGCATCCTGAc                                                                       <  5:141536/149‑1 (MQ=255)                                                                          gAATGCCGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGt                                                                                                                            <  4:28232/96‑1 (MQ=255)                                                                          gAATGCCGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGt                                                                                                                            >  3:28232/1‑96 (MQ=255)                                                                           aaTGCCGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCATGATCAG‑TACTACGTGCGGGCGATCCAGCAAATGCATCCTGAc                                                                       <  8:186624/147‑1 (MQ=255)                                                                           aaTGCCGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCATGATCAG‑TACTACGTGCGGGCGATCCAGCAAATGCATCCTGACCt                                                                     >  8:260594/1‑149 (MQ=255)                                                                               ccGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCATGATCAG‑TACTACGTGCGGGCGATCCAGCAAATGCATCCTGACCTCAg                                                                  <  6:342153/148‑1 (MQ=255)                                                                               ccGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCATGATCAG‑TACTACGTGCGGGCGATCCAGCAAATGCATCCTGACCTCAg                                                                  <  7:61805/148‑1 (MQ=255)                                                                                 gAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCATGATCAG‑TACTACGTGCGGGCGATCCAGCAAATGCATCCTGACCTCAGCtt                                                               <  6:247166/149‑1 (MQ=255)                                                                                 gAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCATGATCAG‑TACTACGTGCGGGCGATCCAGCAAATGCATCCTGACCTCAGCtt                                                               >  7:262574/1‑149 (MQ=255)                                                                                  aGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCATGATCAG‑TACTACGTGCGGGCGATCCAGCAAATGCATCCTGACCTCAGCTTg                                                              <  5:258980/149‑1 (MQ=255)                                                                                   gCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCATGATCAG‑TACTACGTGCGGGCGATCCAGCAAATGCATCCTGACCTCAGCTTGa                                                             <  1:40101/149‑1 (MQ=255)                                                                                     cggcggCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCATGATCAG‑TACTACGTGCGGGCGATCCAGCAAATGCATCCTGACCTCAGCTTGAGc                                                           <  7:260594/149‑1 (MQ=255)                                                                                            cGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCATGATCAG‑TACTACGTGCGGGCGATCCAGCAAATGCATCCTGACCTCAGCTTGAGCTTCTACt                                                    >  5:102804/1‑149 (MQ=255)                                                                                             gAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTAATGATCAGAGATAACGTGCGGGCGATCCAGCAAATGCATCCTGACCTCAGCTTGAGCTTGTACt                                                    >  1:311728/1‑149 (MQ=255)                                                                                                gAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCATGATCAG‑TACTACGTGCGGGCGATCCAGCAAATGCATCCTGACCTCAGCTTGAGCTTCTACTgcgc                                                >  1:295102/1‑149 (MQ=255)                                                                                                  atatGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCATGATCAG‑TACTACGTGCGGGCGATCCAGCAAATGCATCCTGACCTCAGCTTGAGCTTCTACTGCGCgg                                              <  4:69004/149‑1 (MQ=255)                                                                                                   tatGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGAGCGCATTCATGATCAG‑TACTACGTGCGGGCGATCCAGCAAATGCATCCTGACCTCAGCTTGAGCTACTACTGCGCGGt                                             >  5:3851/1‑149 (MQ=255)                                                                                                     tGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCATGATCAG‑TACTACGTGCGGGCGATCCAGCAAATGCATCCTGACCTCAGCTTGAGCTTCTACTGCGCGGTg                                            <  6:169753/148‑1 (MQ=255)                                                                                                     tGCCCGCCACAGCGGTGTGCCCGTGGCTGCGCTCGAGCACCGGGTGTGCGCCGCCCACCCGCTGGCCGTGCGCATTCATGATCAG‑TACTACGTGCGGGCGATCCAGCAAATGCATCCTGACCTCAGCTTGAGCTTCTACTGCGCGGTgg                                           <  6:188690/149‑1 (MQ=255)                                                                                                        ccGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCATGATCAG‑TACTACGTGCGGGCGATCCAGCAAATGCATCCTGACCTCAGCTTGAGCTTCTACTGCGCGGTGgag                                         >  7:133766/1‑148 (MQ=255)                                                                                                                      tgCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCATGATCAG‑TACTACGTGCGGGCGATCCAGCAAATGCATCCTGACCTCAGCTTGAGCTTCTACTGCGCGGTGGAGAACGGCATCGTGCt                           <  5:176915/148‑1 (MQ=255)                                                                                                                             gaatGCGCACGGCCAGCGGGTGGGCGGCGAACACCCGGTGATCGAGCGCATTCATGATCAG‑TACTACGTGCGGGCGATCCAGCAAATGCATCCTGACCTCAGCTTGAGCTTCTACTGCGCGGTGGAGAACGGCATCGTGCTCACGGCCa                   >  3:71599/4‑149 (MQ=25)                                                                                                                               atGCGCACGGCCAGCGGGTGGGCGGCGAACACCCGGTGATCGAGCGCATTCATGATCAG‑TACTACGTGCGGGCGATCCAGCAAATGCATCCTGACCTCAGCTTGAGCTTCTACTGCGCGGTGGAGAACGGCATCGTGCTCACGGCCAt                  >  6:145232/2‑148 (MQ=25)                                                                                                                                   gcTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCATGATCAG‑TACTACGTGCGGGCGATCCAGCAAATGCATCCTGACCTCAGCTTGAGCTTCTACTGCGCGGTGGAGAACGGCATCGTACTCACGGCCATGAccc             >  3:99287/1‑149 (MQ=255)                                                                                                                                              gggTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCATGATCAG‑TACTACGTGCGGGCGATCCAGCAAATGCATCCTGACCTCAGCTTGAGCTTCTACTGCGCGGTGGAGAACGGCATCGTGCTCACGGCCATGACCCCGGGCCCgctg  <  8:133766/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                ATCACATCCTGCCGCGCCAGGAAAAGCTGGTGGTAACCGGTGCCGACAGCGCCCAGCGCCGGGTTTTCGAACTGAATGCCGAGCCGGCGGCCGAGGAATATGCCCGCCACATCGGTGTGCCCGTGGCTGCGCTCGATCACCGGGTGTTCGCCGCCCACCCGCTGGCCGTGCGCATTCATGATCAG‑TACTACGTGCGGGCGATCCAGCAAATGCATCCTGACCTCAGCTTGAGCTTCTACTGCGCGGTGGAGAACGGCATCGTGCTCACGGCCATGACCCCGGGCCCGCTG  >  NC_002947/3050788‑3051077

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