Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,335,984	G→T	10.2%	R122R (CGG→AGG)	PP_3810 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,335,984	0	G	T	10.2%	116.3 / 2.4	39	R122R (CGG→AGG)	PP_3810	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (18/17);  new base T (2/2);  total (20/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.77e-01

GCCTGGCGCAGGGTGGGTTTCAAGTAGACGTAGCAAGCGCTGTAAAAATTAACGGAATTGACGCGCCCCATGTGGCCACGGCCGATGCCGGTATTAACCCTCAGGTATCACGCAACAGGTCGTGCGTATCGAGCAACCGGAACGCAACCTGCGGATTACGCTCCAGCCCGCGCCGGATCGCCGCCGGGATAGACTGCCGGGTCTTGCGGCACAGCCCAGGCTGATCGTCCAGTTCGATGACGATGCCGCGCATGGTCCGTACCTCGTTGTAACCGGGGGCGATGTG  >  NC_002947/4335846‑4336131                                                                                                                                           |                                                                                                                                                    gcCTGGCGCAGGGTGGGTTTCAAGTAGACGTAGCAAGCGCTGTAAAAATTAACGGAATTGACGCGCCCCATGTGGCCACGGCCGATGCCGGTATTAACCCTCAGGTATCACGCAACAGGTCGTGCGTATCGAGCAACCGGAACGCAACCt                                                                                                                                          >  2:954632/1‑150 (MQ=255)   cTGGCGCAGGGTGGGTTTCAAGTAGACGTAGCAAGCGCTGTAAAAATTAACGGAATTGACGCGCCCCATGTGGCCACGGCCGATGCCGGTATTAACCCTCAGGTATCACGCAACAGGTCGTGCGTATCGAGCAACCGGAACGCAACCTGc                                                                                                                                        >  1:172687/1‑150 (MQ=255)   cTGGCGCAGGGTGGGTTTCAAGTAGACGTAGCAAGCGCTGTAAAAATTAACGGAATTGACGCGCCCCATGTGGCCACGGCCGATGCCGGTATTAACCCTCAGGTATCACGCAACAGGTCGTGCGTATCGAGCAACCGGAACGCAACCTGc                                                                                                                                        <  1:399916/150‑1 (MQ=255)          aGGGTGGGTTTCAAGTAGACGTAGCAAGCGCTGTAAAAATTAACGGAATTGACGCGCCCCATGTGGCCACGGCCGATGCCGGTATTAACCCTCAGGTATCACGCAACAGGTCGTGCGTATCGAGCAACCGGAACGCAACCTGCGGATTAc                                                                                                                                 >  2:281622/1‑150 (MQ=255)              tgggtTTCAAGTAGACGTAGCAAGCGCTGTAAAAATTAACGGAATTGACGCGCCCCATGTGGCCACGGCCGATGCCGGTATTAACCCTCAGGTATCACGCAACAGGTCGTGCGTATCGAGCAACCGGAACGCAACCTGCGGATTACGCTc                                                                                                                             <  1:215841/150‑1 (MQ=255)                      aaGTAGACGTAGCAAGCGCTGTAAAAATTAACGGAATTGACGCGCCCCATGTGGCCACGGCCGATGCCGGTATTAACCCTCAGGTATCACGCAACAGGTCGTGCGTATCGAGCAACCGGAACGc                                                                                                                                               >  1:913207/1‑124 (MQ=255)                      aaGTAGACGTAGCAAGCGCTGTAAAAATTAACGGAATTGACGCGCCCCATGTGGCCACGGCCGATGCCGGTATTAACCCTCAGGTATCACGCAACAGGTCGTGCGTATCGAGCAACCGGAACGc                                                                                                                                               <  2:913207/124‑1 (MQ=255)                          aGACGTAGCAAGCGCTGTAAAAATTAACGGAATTGACGCGCCCCATGTGGCCACGGCCGATGCCGGTATTAACCCTCAGGTATCACGCAACAGGTCGTGCGTATCGAGCAACCGGAACGCAACCTGCGGATTACGCTCCAGCCCGCGCCg                                                                                                                 <  2:451452/150‑1 (MQ=255)                            aCGTAGCAAGCGCTGTAAAAATTAACGGAATTGACGCGCCCCATGTGGCCGCGGCCGATGCCGGTATTAACCCTCAGGTATCACGCAACAGGTCGTGCGTATCGAGCAACCTGAACGCAACCGGCGGATTACGCTCCAGCCCGCGCCGGa                                                                                                               >  2:590964/1‑150 (MQ=255)                            aCGTAGCAAGCGCTGTAAAAATTAACGGAATTGACGCGCCCCATGTGGCCACGGCCGATGCCGGTATTAACCCTCAGGTATCACGCAACAGGTCGTGCGTATCGAGCAACCTGAACGCAACCGGCGGATTACGCTCCAGCCCGCGCCGGa                                                                                                               >  2:590778/1‑150 (MQ=255)                                                   aaCGGAATTGACGCGCCCCATGTGGCCACGGCCGATGCCGGTATTAACCCTCAGGTATCACGCAACAGGTCGTGCGTATCGAGCAACCGGAACGCAACCTGCGGATTACGCTCCAGCCCGCGCCGGATCGCCGCCGGGATAGACTGCCgg                                                                                        <  1:281622/150‑1 (MQ=255)                                                    aCGGAATTGACGCGCCCCATGTGGCCACGGCCGATGCCGGTATTAACCCTCAGGTATCACGCAACAGGTCGTGCGTATCGAGCAACCGGAACGCAACCTGc                                                                                                                                        <  2:74243/101‑1 (MQ=255)                                                    aCGGAATTGACGCGCCCCATGTGGCCACGGCCGATGCCGGTATTAACCCTCAGGTATCACGCAACAGGTCGTGCGTATCGAGCAACCGGAACGCAACCTGc                                                                                                                                        >  1:74243/1‑101 (MQ=255)                                                        aaTTGACGCGCCCCATGTGGCCACGGCCGATGCCGGTATTAACCCTCAGGTATCACGCAACAGGTCGTGCGTATCGAGCAACCGGAACGCAACCTGCGGATTACGCTCCAGCCCGCGCCGGATCGCCGCCGGGATAGACTGCCGGGTCtt                                                                                   <  1:954632/150‑1 (MQ=255)                                                         aTTGACGCGCCCCATGTGGCCACGGCCGATGCCGGTATTAACCCTCAGGTATCACGCAACAGGTCGTGCGTATCGAGCAACCGGAACGCAACCTGCGGATTACGCTCCAGCCCGCGCCGGATCGCCGCCGGGATAGACTGCCGGGTCTTg                                                                                  >  1:580580/1‑150 (MQ=255)                                                          ttGACGCGCCCCATGTGGCCACGGCCGATGCCGGTATTAACCCTCAGGTATCACGCAACAGGTCGTGCGTATCGAGCAACCGGAACGCAACCTGCGGATTACGCTCCAGCCCGCGCCGGATCGCCGCCGGGATAGACTGCCGGGTCTTGc                                                                                 <  1:525821/150‑1 (MQ=255)                                                                    ccATGTGGCCACGGCCGATGCCGGTATTAACCCTCAGGTATCACGCAACAGGTCGTGCGTATCGAGCAACCGGAACGCAACCTGCGGATTACGCTCCAGCCCGCGCCGGATCGCCGCCGGGATAGACTGCCGGGTCTTGCGGCACAGccc                                                                       >  2:327116/1‑150 (MQ=255)                                                                            ccACGGCCGATGCCGGTATTAACCCTCAGGTATCACGCAACAGGTCGTGCGTATCGAGCAACCGGAACGCAACCTGCGGATTACGCTCCAGCCCGCGCCGGATCGCCGCCGGGATAGACTGCCGGGTCTTGCGGCACAGCCCAGGCTgct                                                               >  1:244415/1‑148 (MQ=255)                                                                            ccACGGCCGATGCCGGTATTAACCCTCAGGTATCACGCAACAGGTCGTGCGTATCGAGCAACCGGAACGCAACCTGCGGATTACGCTCCAGCCCGCGCCGGATCGCCGCCGGGATAGACTGCCGGGTCTTGCGGCACAGCCCAGGCTGAt                                                               <  1:327116/150‑1 (MQ=255)                                                                             cACGGCCGATGCCGGTATTAACCCTCAGGTATCACGCAACAGGTCGTGCGTATCGAGCAACCGGAACGCAACCTGCGGATTACGCTCCAGCCCGCGCCGGATCGCCGCCGGGATAGACTGCCGGGTCTTGCGGCACAGCCCAGGCTGATc                                                              >  2:472049/1‑150 (MQ=255)                                                                                ggCCGATGCCGGTATTAACCCTCAGGTATCACGCAACAGGTCGTGCGTATCGAGCAACCTGAACGCAACCTGCGGATTACGCTCCAGCCCGCGCCGGATCGCCGCCGGGATAGACTGCCGGGTCTTGCGGCACAGCCCAGGCTGATCGTc                                                           <  1:590964/150‑1 (MQ=255)                                                                                ggCCGATGCCGGTATTAACCCTCAGGTATCACGCAACAGGTCGTGCGTATCGAGCAACCTGAACGCAACCTGCGGATTACGCTCCAGCCCGCGCCGGATCGCCGCCGGGATAGACTGCCGGGTCTTGCGGCACAGCCCAGGCTGATCGTc                                                           <  1:590778/150‑1 (MQ=255)                                                                                 gCCGATGCCGGTATTAACCCTCAGGTATCACGCAACGGTATGTGCGTATCGAGCAACCGGAACGCAACCTGCGGATTACGCTCCAGCCCGCGCCGGATCGCCGCCGGGATAGACTGCCGGGTCTTGCGGCACAGCCCAGGCTGATCGTcc                                                          <  1:57878/150‑1 (MQ=255)                                                                                       gCCGGTATTAACCCTCAGGTATCACGCAACAGGTCGTGCGTATCGAGCAACCGGAACGCAACCTGCGGATTACGCTCCAGCCCGCGCCGGATCGCCGCCGGGATAGACTGCCGGGTCTTGCGGCACAGCCCAGGCTGATCGTCCAGTTCg                                                    >  1:786618/1‑150 (MQ=255)                                                                                        ccGGTATTAACCCTCAGGTATCACGCAACAGGTCGTGCGTATCGAGCAACCGGAACGCAACCTGCGGATTACGCTCCAGCCCGCGCCGGATCGCCGCCGGGATAGACTGCCGGGTCTTGCGGCACAGCCCAGGCTGATCGTCCAGTTCGa                                                   <  1:973035/150‑1 (MQ=255)                                                                                                  cccTCAGGTATCACGCAACAGGTCGTGCGTATCGAGCAACCGGAACGCAACCTGCGGATTACGCTCCAGCCCGCGCCGGATCGCCGCCGGGATAGACTGCCGGGTCTTGCGGCACAGCCCAGGCTGATCGTCCAGTTCGATGACGATGcc                                         <  1:1010489/150‑1 (MQ=255)                                                                                                      cAGGTATCACGCAACAGGTCGACTCCAGNGAGCAACCGGAACGCAACCTGCGGATTACGCTCCAGCCCGCGCCGGATCGCCGCCGGGATAGACTGCCGGGTCTTGCGGCACAGCCCAGGCTGATCGTCCAGTTCGATGACGATGCCncgc                                     >  1:997035/1‑150 (MQ=255)                                                                                                       aGGTATCACGCAACAGGTCGTGCGTATCGAGCAACCGGAACGCAACCTGCGGATTACGCTCCAGCCCGCGCCGGATCGCCGCCGGGATAGACTGCCGGGTCTTGCGGCACAGCCCAGGCTGATCGTCCAGTTCGATGACGATGCCGCGCa                                    <  2:203437/150‑1 (MQ=255)                                                                                                          tATCACGCAACAGGTCGTGCGTATCGAGCAACCGGAACGCAACCTGCGGATTACGCTCCAGCCCGCGCCGGATCGCCGCCGGGATAGACTGCCGGGTCTTGCGGCACAGCCCAgg                                                                    >  1:883375/1‑115 (MQ=255)                                                                                                          tATCACGCAACAGGTCGTGCGTATCGAGCAACCGGAACGCAACCTGCGGATTACGCTCCAGCCCGCGCCGGATCGCCGCCGGGATAGACTGCCGGGTCTTGCGGCACAGCCCAgg                                                                    <  2:883375/115‑1 (MQ=255)                                                                                                          tATCACGCAACAGGTCGTGCGTATCGAGCAACCGGAACGCAACCTGCGGATTACGCTCCAGCCCGCGCCGGATCGCCGCCGGGATAGACTGCCGGGTCTTGCGGCACAGCCCAGGCTGATCGTCCAGTTCGATGACGATGCCGCGCATgg                                 >  1:767674/1‑150 (MQ=255)                                                                                                            tCACGCAACAGGTCGTGCGTATCGAGCAACCGGAACGCAACCTGCGGATTACGCTCCAGCCCGCGCCGGATCGCCGCCGGGATAGACTGCCGGGTCTTGCGGCACAGCCCAGGCTGATCGTCCAGTTCGATGACGATGCCGCGCATGGTc                               >  1:67306/1‑150 (MQ=255)                                                                                                              aCGCAACAGGTCGTGCGTATCGAGCAACCGGAACGCAACCTGCGGATTACGCTCCAGCCCGCGCCGGATCGCCGCCGGGATAGACTGCCGGGTCTTGCGGCACAGCCCAGGCTGATCGTCCAGTTCGATGACGATGCCGCGCATGGTCCg                             >  2:755712/1‑150 (MQ=255)                                                                                                                gCAACAGGTCGTGCGTATCGAGCAACCGGAACGCAACCTGCGGATTACGCTCCAGCCCGCGCCGGATCGCCGCCGGGATAGACTGCCGGGTCTTGCGGCACAGCCCAGGCTGATCGTCCAGTTCGATGACGATGCCGCGCATGGTCCGTa                           >  1:1085949/1‑150 (MQ=255)                                                                                                                               tATCGAGCAACCGGAACGCAACCTGCGGATTACGCTCCAGCCCGCGCCGGATCGCCGCCGGGATAGACTGCCGGGTCTTGCGGCACAGCCCAGGCTGATCGTCCAGTTCGATGACGATGCCGCGCATGGTCCGTACCTCGTTGTAACCgg            <  2:767674/150‑1 (MQ=255)                                                                                                                                  cGAGCAACCGGAACGCAACCTGCGGATTACGCTCCAGCCCGCGCCGGATCGCCGCCGGGATAGACTGCCGGGTCTTGCGGCACAGCCCAGGCTGATCGTCCAGTTCGATGACGATGCCGCGCATGGTCCGTACCTCGTTGTAACCggggg         <  2:1044434/150‑1 (MQ=255)                                                                                                                                  cGAGCAACCGGAACGCAACCTGCGGATTACGCTCCAGCCCGCGCCGGATCGCCGCCGGGATAGACTGCCGGGTCTTGCGGCACAGCCCAGGCTGATCGTCCAGTTCGATGACGATGCCGCGCATGGTCCGTACCTCGTTGTAACCggggg         >  1:707792/1‑150 (MQ=255)                                                                                                                                  cGAGCAACCGGAACGCAAAATACGGATTACGCTCCAGACCGCGCCAGATCGCCGCCGGTATAGACTGCCGGGTCTTGCGGCACAGCCCAGGCTGATCGTCCAGTTCGATGACGATGCCGCGCATGGTCCGTACCTCGTTGTAACCggggg         <  2:67306/150‑1 (MQ=255)                                                                                                                                         ccGGAACGCAACCTGCGGATTACGCTCCAGCCCGCGCCGGATCGCCGCCGGGATAGACTGCCGGGTCTTGCGGCACAGCCCAGGCTGATCGTCCAGTTCGATGACGATGCCGCGCATGGTCCGTACCTCGTTGTGACAGGGGGCGAtggg  >  2:905820/1‑148 (MQ=255)                                                                                                                                           |                                                                                                                                                    GCCTGGCGCAGGGTGGGTTTCAAGTAGACGTAGCAAGCGCTGTAAAAATTAACGGAATTGACGCGCCCCATGTGGCCACGGCCGATGCCGGTATTAACCCTCAGGTATCACGCAACAGGTCGTGCGTATCGAGCAACCGGAACGCAACCTGCGGATTACGCTCCAGCCCGCGCCGGATCGCCGCCGGGATAGACTGCCGGGTCTTGCGGCACAGCCCAGGCTGATCGTCCAGTTCGATGACGATGCCGCGCATGGTCCGTACCTCGTTGTAACCGGGGGCGATGTG  >  NC_002947/4335846‑4336131

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