Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,247,865 G→C 11.8% R282G (CGG→GGG)  PP_4623 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,247,8650GC11.8% 98.5 / 5.9 34R282G (CGG→GGG) PP_4623hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (15/15);  new base C (2/2);  total (17/17)
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.32e-01

ATAGGCTTGCAGTTCCTCGCCCTGCATGAACCGCCGCCATTCGGGCAGCAGGCTGGTGCGGTACGCTCGCACATGCCGCGTCGCACGTTTCTGCAGCGTAGCCACCAGAGGCTCAAGCCGAGCCGCACTGTCGAGCTGGGCCAGACGCCGGGCGTCATCGACCTTGACCCGCAGGCTGCGCCGGGTGGCTTCAAGCTGGATAGCCTGCAGATCACGGATGAAATACCTGTGAAGTGGCACCTGCAGCGGCACGAAATCGAGTGTCAAACTCTCCGTAGACAACGCCGCCCGCT  >  NC_002947/5247716‑5248008
                                                                                                                                                     |                                                                                                                                               
aTAGGCTTGCAGTTCCTCGCCCTGCATGAACCGCCGCCATTCGGGCAGCAGGCTGGTGCGGTACGCTCGCACATGCCGCGTCGCACGTTTCTGCAGCGTAGCCACCAGAGGCTCAAGCCGAGCCGCACTGTCGAGCTGGGCCAGACGCCg                                                                                                                                                 <  1:177175/150‑1 (MQ=255)
 tAGGCTTGCAGTTCCTCGCCCTGCATGAACCGCCGCCATTCGGGCAGCAGGCTGGTGCGGTACGCTCGCACATGCCGCGTCGCACGTTTCTGCAGCGTAGCCACCAGAGGCTCAAGCCGAGCCGCACTGTCGAGCTGGGCCAGACGCCgg                                                                                                                                                <  2:508334/150‑1 (MQ=255)
                  gCCCTGCATGAACCGCCGCCATTAACGCGTTAGGCTGGTGCGGTACGCTCGCACATGCCGCGTCGCACGTTTCTGCAGCGTAGCCACCAGAGGCTCAAGCCGAGCCGCACTGTCGAGCTGGGCCAGACGCCGGGCGTCATCGACCTTGAc                                                                                                                               >  1:420270/1‑150 (MQ=255)
                        cATGAACCGCCGCCATTCGGGCAGCAGGCTGGTGCGGTACGCTCGCACATGCCGCGTCGCACGTTTCTGCAGCGTAGCCACCAGAGGCTCAAGCCGAGCCGCACTGTCGAGCTGGGCCAGACGCCGGGCGTCATCGACCTTGACCCGCAg                                                                                                                         <  2:592008/150‑1 (MQ=255)
                              ccgccgCCATTCGGGCAGCAGGCTGGTGCGGTACGCTCGCACATGCCGCGTCGCACGTTTCTGCAGCGTAGCCACCAGAGGCTCAAGCCGAGCCGCACTGTCGAGCTGGGCCAGACGCCCGGCGTctggc                                                                                                                                       <  2:229546/130‑5 (MQ=255)
                              ccgccgCCATTCGGGCAGCAGGCTGGTGCGGTACGCTCGCACATGCCGCGTCGCACGTTTCTGCAGCGTAGCCACCAGAGGCTCAAGCCGAGCCGCACTGTCGAGCTGGGCCAGACGCCCGGCGTctggc                                                                                                                                       >  1:229546/1‑126 (MQ=255)
                                       ttCGGGCAGCAGGCTGGTGCGGTACGCTCGCACATGCCGCGTCGCACGTTTCTGCAGCGTAGCCACCAGAGGCTCAAGCCGAGCCGCACTGTCGAGCTGGGCCAGACGCCGGGCGTCATCGACCTTGACCCGCAGGCTGCGCCGGGTGGc                                                                                                          <  1:534966/150‑1 (MQ=255)
                                                            gTACGCTCGCACATGCCGCGTCGCACGTTTCTGCAGCGTAGCCACCAGAGGCTCAAGCCGAGCCGCACTGTCGAGCTGGGCCAGACGCCGGGCGTCATCGACCTTGACCCGCAGGCTGCGCCGGGTGGCTTCAAGCTGGATAGCCTGCAg                                                                                     >  2:509824/1‑150 (MQ=255)
                                                                  tCGCACATGCCGCGTCGCACGTTTCTGCAGCGTAGCCACCAGAGGCTCAAGCCGAGCCGCACTGTCGAGCTGGGCCAGACGCCGGGCGTCATCGACCTTGACCCGCAGGCTGCGCCGGGTGGCTTCAAGCTGGATAGCCTGCAGATCAct                                                                               <  2:423669/150‑2 (MQ=255)
                                                                  tCGCACATGCCGCGTCGCACGTTTCTGCAGCGTAGCCACCAGAGGCTCAAGCCGAGCCGCACTGTCGAGCTGGGCCAGACGCCGGGCGTCATCGACCTTGACCCGCAGGCTGCGCCGGGTGGCTTCAAGCTGGATAGCCTGCAGATCACg                                                                               >  2:20439/1‑150 (MQ=255)
                                                                    gCACATGCCGCGTCGCACGTTTCTGCAGCGTAGCCACCAGAGGCTCAAGCCGAGCCGCACTGTCGAGCTGGGCCAGACGCCGGGCGTCATCGACCTTGACCCGCAGGCTGCGCCGGGTGGCTTCAAGCTGGATAGCCTGCAGATCACGGa                                                                             <  2:216202/150‑1 (MQ=255)
                                                                    gCACATGCCGCGTCGCACGTTTCTGCAGCGTAGCCACCAGAGGCTCAAGCCGAGCCGCACTGTCGAGCTGGGCCAGACGCCGGGCGTCATCGACCTTGACCCGCAGGCTGCGCCGGGTGGCTTCAAGCTGGATAGCCTGCAGATCACGGa                                                                             <  2:420270/150‑1 (MQ=255)
                                                                     cacaTGCCGCGTCGCACGTTTCTGCAGCGTAGCCACCAGAGGCTCAAGCCGAGCCGCACTGTCGAGCTGGGCCAGACGCCGGGCGTCATCGACCTTGACCCGCAGGCTGCGCCGGGTGGCTTCAAGCTGGATAGCCTGCAGATCACGGAt                                                                            <  1:361054/150‑1 (MQ=255)
                                                                      acaTGCCGCGTCGCACGTTTCTGCAGCGTAGCCACCAGAGGCTCAAGCCGAGCCGCACTGTCGAGCTGGGCCAGACGCCGGGCGTCATCGAc                                                                                                                                     <  2:320059/92‑1 (MQ=255)
                                                                      acaTGCCGCGTCGCACGTTTCTGCAGCGTAGCCACCAGAGGCTCAAGCCGAGCCGCACTGTCGAGCTGGGCCAGACGCCGGGCGTCATCGACCTTGAc                                                                                                                               >  1:5857/1‑98 (MQ=255)
                                                                      acaTGCCGCGTCGCACGTTTCTGCAGCGTAGCCACCAGAGGCTCAAGCCGAGCCGCACTGTCGAGCTGGGCCAGACGCCGGGCGTCATCGACCTTGAc                                                                                                                               <  2:5857/98‑1 (MQ=255)
                                                                      acaTGCCGCGTCGCACGTTTCTGCAGCGTAGCCACCAGAGGCTCAAGCCGAGCCGCACTGTCGAGCTGGGCCAGACGCCGGCCGGCATCGAc                                                                                                                                     >  1:320059/1‑92 (MQ=255)
                                                                        aTGCCGCGTCGCACGTTTCTGCAGCGTAGCCACCAGAGGCTCAAGCCGAGCCGCACTGTCGAGCTGGGCCAGACGCCGGGCGTCATCGACCTTGACCCGCAGGCTGCGCCGGGTGGCTTCAAGCTGGATAGCCTGCAGATCACGGATGaa                                                                         <  2:158035/150‑1 (MQ=255)
                                                                        aTGCCGCGTCGCACGTTTCTGCAGCGTAGCCACCAGAGGCTCAAGCCGAGCCGCACTGTCGAGCTGGGCCAGACGCCGGGCGTCATCGACCTTGACCCGCAGGCTGCGCCGGGTGGCTTCAAGCTGGATAGCCTGCAGATCACGGATGaa                                                                         >  1:366388/1‑150 (MQ=255)
                                                                                 cGCACGTTTCTGCAGCGTAGCCACCAGAGGCTCAAGCCGAGCCGCACTGTCGAGCTGGGCCAGACGCCGGGCGTCATCGACCTTGACCCGCAGGCTGCGCCGGGTGGCTTCAAGCTGGATAGCCTGCAGATCACGGATGAAATACCtgtg                                                                <  1:20439/150‑1 (MQ=255)
                                                                                                        ccaGAGGCTCAAGCCGAGCCGCACTGTCGAGCTGGGCCAGACGCCGGGCGTCATCGACCTTGACCCGCAGGCTGCGCCGGGTGGCTTCAAGCTGGATAGCCTGCAGATCACGGATGAAATACCTGTGAAGTGGCACCTGCAGCGGCACGa                                         >  1:213347/1‑150 (MQ=255)
                                                                                                           gagGCTCAAGCCGAGCCGCACTGTCGAGCTGGGCCAGACGCCGGGCGTCATCGACCTTGACCCGCAGGCTGCGCCGGGTGGCTTCAAGCTGGATAGCCTGCAGATCACGGATGAAATACCTGTGAAGTg                                                           >  1:260866/1‑129 (MQ=255)
                                                                                                           gagGCTCAAGCCGAGCCGCACTGTCGAGCTGGGCCAGACGCCGGGCGTCATCGACCTTGACCCGCAGGCTGCGCCGGGTGGCTTCAAGCTGGATAGCCTGCAGATCACGGATGAAATACCTGTGAAGTg                                                           <  2:260866/129‑1 (MQ=255)
                                                                                                              gCTCAAGCCGAGCCGCACTGTCGAGCTGGGCCAGACGCCGGGCGTCATCGACCTTGACCCGCAGGCTGCGCCGGGTGGCTTCAAGCTGGATAGCCTGCAGATCACGGATGAAATACCTGTGAAGTGGCACCTGCAGCGGCACGAAATCGa                                   <  2:269440/150‑1 (MQ=255)
                                                                                                                tCAAGCCGAGCCGCACTGTCGAGCTGGGCCAGACGCCGGGCGTCATCGACCTTGACCCGCAGGCTGCGCCGGGTGGCTTCAAGGTGGATAGCCTGCAGATCACGGATGAAATACCTGTGAAGTg                                                           <  1:588639/124‑1 (MQ=255)
                                                                                                                tCAAGCAGAGCCGCACTGTCGAGCTGGGCCAGACGCCGGGCGTCATCGACCTTGACCCGCAGGCTGCGCCGGGTGGCTTCAAGCTGGATAGCCTGCAGATCACGGATGAAATACCTGTGAAGTg                                                           >  2:588639/1‑124 (MQ=255)
                                                                                                                                tGTCGAGCTGGGCCAGACGCCGGGCGTCATCGACCt                                                                                                                                   <  2:110737/36‑1 (MQ=255)
                                                                                                                                tGTCGAGCTGGGCCAGACGCCGGGCGTCATCGACCt                                                                                                                                   >  1:110737/1‑36 (MQ=255)
                                                                                                                                tGTCGAGCTGGGCCAGACGCCGGGCGTCATCGACCTTGACCCGCAGGCTGCGCCGGGTGGCTTCAAGCTGGATAGCCTGCAGATCACGGATGAAATACCTGTGAAGTGGCACCTGCAGCGGCACGAAATAGAGTGTCAAACTCTCCGTAg                 >  2:108948/1‑150 (MQ=255)
                                                                                                                                       cTGGGCCAGACGCCGGGAGTCATCGACCTTGACCCGCAGGCTGCGCCGGGTGGCTTCAAGCTGGATAGCCTGCAGATCACGGATGAAATACCTGTGAAGTGGCACATGCAGGGGCACGAAAATGAGCGAGAAACTCTCCGTAGACAAcgc          >  2:314888/1‑150 (MQ=255)
                                                                                                                                             cAGACGCCGGGCGTCATCGACCTTGACCCGCAGGCTGCGCCGGGTGGCTTCAAGCTGGATAGCCTGCAGATCACGGATGAAATACCTGTGAAGTGGCACCTGCAGCGGCACGAAATCGAGTGTCAAACTCTCCGTAGACAACGccgcccg    >  1:55038/1‑150 (MQ=255)
                                                                                                                                             atgACGCCCGGCGTCATCGACCTTGACCCGCAGGCTGCGCCGGGTGGCTTCAAGCTGGATAGCCTGCAGa                                                                                    >  1:158517/3‑70 (MQ=255)
                                                                                                                                             atgACGCCCGGCGTCATCGACCTTGACCCGCAGGCTGCGCCGGGTGGCTTCAAGCTGGATAGCCTGCAGa                                                                                    <  2:158517/68‑1 (MQ=255)
                                                                                                                                               gACGCCGGGCGTCATCGACCTTGACCCGCAGGCTGCGCCGGGTGGCTTCAAGCTGGATAGCCTGCAGATCACGGATGAAATACCTGTGAAGTGGCGCCTGCAGCGGCACGAAATCGAGTGTCAAACTCTCCGTAGACAACGccgccc     >  1:66249/1‑147 (MQ=255)
                                                                                                                                               gACGCCGGGCGTCATCGACCTTGACCCGCAGGCTGCGCCGGGTGGCTTCAAGCTGGATAGCCTGCAGATCACGGATGAAATACCTGTGAAGTGGCGCCTGCAGCGGCACGAAATCGAGTGTCAAACTCTCCGTAGACAACGccgccc     <  2:66249/147‑1 (MQ=255)
                                                                                                                                               gACGCCGGGCGTCATCGACCTTGACCCGCAGGCTGCGCCGGGTGGCTTCAAGCTGGATAGCCTGCAGATCACGGATGAAATACCTGTGAAGTGGCACCTGCAGCGGCACGAAATCGAGTGTCAAACTCTCCGTAGACAACGCCGCCCGct  >  1:18203/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
ATAGGCTTGCAGTTCCTCGCCCTGCATGAACCGCCGCCATTCGGGCAGCAGGCTGGTGCGGTACGCTCGCACATGCCGCGTCGCACGTTTCTGCAGCGTAGCCACCAGAGGCTCAAGCCGAGCCGCACTGTCGAGCTGGGCCAGACGCCGGGCGTCATCGACCTTGACCCGCAGGCTGCGCCGGGTGGCTTCAAGCTGGATAGCCTGCAGATCACGGATGAAATACCTGTGAAGTGGCACCTGCAGCGGCACGAAATCGAGTGTCAAACTCTCCGTAGACAACGCCGCCCGCT  >  NC_002947/5247716‑5248008

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