Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,210,934 G→C 8.0% G109R (GGC→CGC)  PP_1954 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,210,9340GC 8.0% 151.8 / 4.7 50G109R (GGC→CGC) PP_1954hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (25/21);  new base C (2/2);  total (27/23)
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.17e-01

CGCTGTTTTCAAAGTTCTCTCGCTTACCAAGGCTTTCACCACGACTCTGGTGCTGGAAGCAATCGACCTCGGCTTATTGTCCCTGACTACCCGAGTTGTTGATGTAATCCCGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAGGCCTTTCGCCTACTCCACAGCCTCTTCCGTATGACCAACTGCACAATTTTAGCGAGACTATCCAGGCAATCTGCGCCCTGGATGTAGTTGGCGAGC  >  NC_002947/2210792‑2211078
                                                                                                                                              |                                                                                                                                                
cgaggatttCAAAGTGCTCTCGCCTACCAAGGCTTTCACCACGACTCTGGTGCTGGAAGCAATCGACCTCGGCTTATTGTCCCAGACTACCCGAGTTGTTGATGTAATCCCGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGGa                                                                                                                                            <  5:279822/143‑1 (MQ=255)
       ttCAAAGTTCTCTCGCTTACCAAGGCTTTCACCACGACTCTGGTGCTGGAAGCAATCGACCTCGGCTTATTGTCCCTGACTACCCGAGTTGTTGGTGTAATCCCGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGGAAAGAGc                                                                                                                                      <  4:47329/148‑1 (MQ=255)
                         aCCAAGGCTTTCACCACGACTCTGGTGCTGGAAGCAATCGACCTCGGCTTATTGTCCCTGACTACCCGAGTTGTTGATGTAATCCCGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATc                                                                                                                    <  7:240913/148‑1 (MQ=255)
                           cAAGGCTTTCACCACGACTCTGGTGCTGGAAGCAATCGACCTCGGCTTATTGTCCCTGACTACCCGAGTTGTTGATGTAATCCCGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGt                                                                                                                 >  4:161225/1‑149 (MQ=255)
                           cAAGGCTTTCACCACGACTCTGGTGCTGGAAGCAATCGACCTCGGCTTATAATCCCTGACTACCCGAGTTGTTGATGTAATCCCGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGt                                                                                                                 >  1:255530/1‑149 (MQ=255)
                            aaGGCTTTCACCACGACTCTTGTGCTGGAAGCAATCGACCTCGGCTTATTGTCCCTGACTACCCGAGTTGTTGATGTAATCCCGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGtt                                                                                                                <  2:93457/149‑1 (MQ=255)
                                       cacGACTCTGGTGCTGGAAGCAATCGACCTCGGCTTATTGTCCCTGACTACCCGAGTTGTTGATGTAATCCCGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACcgcg                                                                                                     >  3:304373/1‑149 (MQ=255)
                                          gACTCTGGTGCTGGAAGCAATCGACCTCGGCTTATTGTCCCTGACTACCCGAGTTGTTGATGTAATCCCGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAg                                                                                                  <  4:304373/149‑1 (MQ=255)
                                               tGGTGCTGGAAGCAATCGACCTCGGCTTATTGTCCCTGACTACCCGAGTTGTTGATGTAATCCCGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTCGTCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAgg                                                                                                 <  8:310868/145‑1 (MQ=255)
                                               tGGTGCTGGAAGCAATCGACCTCGGCTTATTGTCCCTGACTACCCGAGTTGTTGATGTAATCCCGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTCGTCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAgg                                                                                                 <  8:306594/145‑1 (MQ=255)
                                               tGGTGCTGGAAGCAATCGACCTCGGCTTATTGTCCCTGACTACCCGAGTTGTTGATGTAATCCCGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTCGTCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAgg                                                                                                 >  7:310868/1‑145 (MQ=255)
                                               tGGTGCTGGAAGCAATCGACCTCGGCTTATTGTCCCTGACTACCCGAGTTGTTGATGTAATCCCGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTCGTCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAgg                                                                                                 >  7:306594/1‑145 (MQ=255)
                                                     tGGAAGCAATCGACCTCGGCTTATTGTCCCTGACTACCCGAGTTGTTGATGTAATACCGGAGATTTGGGGGCGCGACCTGTTCAGCAGTGGCCGGAAAGCGCGTGTCAATGTCGGTAATCTGTTGACCCACCGCGCAGGCCTTTCGCCt                                                                                       >  5:83999/1‑149 (MQ=255)
                                                       gAAGCAATCGACCTCGGCTTATTGTCCCTGACTACCCGAGTTGTTGATGTAATCCCGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAGGCCTTTCGCCTAc                                                                                     >  1:220884/1‑149 (MQ=255)
                                                       gAAGCAATCGACCTCGGCTTATTGTCCCTGACTACCCGAGTTGTTGATGTAATCCCGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAGGCCTTTCGCCTAc                                                                                     >  6:75421/1‑149 (MQ=255)
                                                        aaGCAATCGACCTCGGCTTATTGTCCCTGACTACCCGAGTTGTTGATGTAATCCCGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAGGCCTTTCGCCTACt                                                                                    <  4:292083/149‑1 (MQ=255)
                                                                 aCCTCGGCTTATTGTCCCTGACTACCCGAGTTGTTGATGTAATCCCGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAGGCCTTTc                                                                                           >  6:318027/1‑133 (MQ=255)
                                                                 aCCTCGGCTTATTGTCCCTGACTACCCGAGTTGTTGATGTAATCCCGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAGGCCTTTc                                                                                           <  5:318027/133‑1 (MQ=255)
                                                                   cTCGGCTTATTGTCCCTGACTACCCGAGTTGTTGATGTAATCCCGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAGGCCTTTCGCCTACTCCACAGCctct                                                                         >  7:228763/1‑149 (MQ=255)
                                                                   cTCGGCTTATTGTCCCTGACTAACCGAGTTGTTGATGTAATCCCGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAGGCCTTTCGCCTACTCCACAGCctct                                                                         >  7:211166/1‑149 (MQ=255)
                                                                      ggCTTATTGTCCCTGACTACCCGAGTTGTTGATGTAATCCCGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAGGCCTTTCGCCTACTCCACAGCCTCtt                                                                        <  3:161225/147‑1 (MQ=255)
                                                                             tGTCCCTGACTACCCGAGTTGTTGATGTAATCCCGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAGGCCTTTCGCCTACTCCACAGCCTCTTCCGTATGAc                                                               >  7:134048/1‑149 (MQ=255)
                                                                                  cTGACTACCCGAGTTGTTGATGTAATCCCGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAGGCCTTTCGCCTACTCCACAGCCTCTTCCGTATGACCAAc                                                           >  8:206617/1‑148 (MQ=255)
                                                                                  cTGACTACCCGAGTTGTTGATGTAATCCCGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAGGCCTTTCGCCTACTCCACAGCCTCTTCCGTATGACCAAc                                                           <  7:206617/148‑1 (MQ=255)
                                                                                   tGACTACCCGAGTTGTTGATGTAATCCCGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAGGCCTTTCGCCTACTCCACAGCCTCTTCCGTATGACCAACTg                                                         >  5:254017/1‑149 (MQ=255)
                                                                                     aCTACCCGAGTTGTTGATGTAATCCCGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAGGCCTTTCGCCTACTCCACAGCCTCTTCCGTATGACCAACTGca                                                       <  5:75421/149‑1 (MQ=255)
                                                                                         cccGAGTTGTTGATGTAATCCCGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAGGCCTTTCGCCTACTCCACAGCCTCTTCCGTATGACCAACTGCACAAt                                                   >  3:179104/1‑149 (MQ=255)
                                                                                           cGAGTTGTTGATGTAATCCCGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAGGCCTTTCGCCTACTCCACAGCCTCTTCCGTATGACCAACTGCACAAttt                                                 >  8:102398/1‑149 (MQ=255)
                                                                                           cGAGTTGTTGATGTAATCCCGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAGGCCTTTCGCCTACTCCACAGCCTCTTCCGTATGACCAACTGAACAAttt                                                 >  5:21031/1‑149 (MQ=255)
                                                                                                 gttgATGTAATCCCGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAGGCCTTTCGCCTACTCCACAGCCTCTTCCGTATGACCAACTGCACAATTTTAGc                                             >  5:141484/1‑147 (MQ=255)
                                                                                                    gATGTAATCCCGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAGGCCTTTCGCCTACTCCACAGCCTCTTCCGTATGACCAACTGCACAATTTTa                                               <  5:90192/142‑1 (MQ=255)
                                                                                                    gATGTAATCCCGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAGGCCTTTCGCCTACTCCACAGCCTCTTCCGTATGACCAACTGCACAATTTTa                                               >  6:90192/1‑142 (MQ=255)
                                                                                                     aTGTAATGCCGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAGGCCTTTCGCCTACTCCACAGCCTCTTCCGTATGa                                                                <  8:247667/124‑1 (MQ=255)
                                                                                                     aTGTAATCCCGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAGGCCTTTCGCCTACTCCACAGCCTCTTCCGTATGACCAACTGCACAATTTTAGCGAGAc                                        >  4:281521/1‑148 (MQ=255)
                                                                                                          aTCCCGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAGGCCTTTCGCCTACTCCACAGCCTCTTCCGTATGACCAACTGCACAATTTTAGCGAGACTATCCa                                  <  6:141484/149‑1 (MQ=255)
                                                                                                             ccGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGTAAAGAGCGTGTCAATGTTGGTCATCTGTTGACCCACCGCGCAGGCCTTTCGCCTACTCCACAGCCTCTTCCGTATGACCAACTGCACAAttt                                                 <  6:210881/131‑1 (MQ=255)
                                                                                                             ccGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTTGGTCATCTGTTGACCCACCGCGCAGGCCTTTCGCCTACTCCACAGCCTCTTCCGTATGACCAACTGCACAAttt                                                 >  5:210881/1‑131 (MQ=255)
                                                                                                              cGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAGGCCTTTCGCCTACTCCACAGCCTCTTCCGTATGACCAACTGCACAATTTTAGCGAGACTATCCAGGc                               >  3:249332/1‑148 (MQ=255)
                                                                                                              cGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAGGCCTTTCGCCTACTCCACAGCCTCTTCCGTATGACCAACTGCACAATTTTAGCGAGACTATCCAGGc                               >  3:220880/1‑148 (MQ=255)
                                                                                                                gAGTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAGGCCTTTCGCCTACTCCACAGCCTCTTCCGTATGACCAACTGCACAATTTTAGCGAGACTATCCAGGCAAt                            <  8:134048/149‑1 (MQ=255)
                                                                                                                  gTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAGGCCTTTCGCCTACTCCACAGCCTCTTCCGTATGACCAACTGCACAATTTTAGCGAGACTATCCAGGCAATCt                          <  4:220880/149‑1 (MQ=255)
                                                                                                                   ttttGGGGGCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAGGCCTTTCGCCTACTCCACAGCCTCTTCCGTATGACCAACTGCACAATTTTAGCGAGACTATCCAGGCAATCTg                         >  6:165858/1‑149 (MQ=255)
                                                                                                                       gggggCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAGGCCTTTCGCCTACTCCACAGCCTCTTCCGTATGACCAACTGCACAATTTTAGCGAGACTATCCAGGCAATCTgcgc                      <  2:255530/148‑1 (MQ=255)
                                                                                                                         gggCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAGGCCTTTCGCCTACTCCACAGCCTCTTCCGTATGACCAACTGCACAATTTTAGCGAGACTATCCAGGCAATCTGCGCCCt                   <  2:175420/149‑1 (MQ=255)
                                                                                                                           gcgcGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAGGCCTTTCGc                                                                                         >  3:282453/1‑77 (MQ=255)
                                                                                                                           gcgcGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAGGCCTTTCGc                                                                                         <  4:282453/77‑1 (MQ=255)
                                                                                                                                    gTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAGGCCTTTCGCCTACTCCACAGCCTCTTCCGTATGACCAACTGCACAATTTTAGCGAGACTATCCAGGCAATCTGCGCCCTg                  >  6:38999/1‑139 (MQ=255)
                                                                                                                                    gTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAGGCCTTTCGCCTACTCCACAGCCTCTTCCGTATGACCAACTGCACAATTTTAGCGAGACTATCCAGGCAATCTGCGCCCTg                  <  5:38999/139‑1 (MQ=255)
                                                                                                                                        cGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAGGCCTTTCGCCTACTCCACAGCCTCTTCCGTATGACCAACTGCACAATTTTAGCGAGACTATCCAGGCAATCTGCGCCCTGGATGTAGTTGGCGa    <  4:179104/149‑1 (MQ=255)
                                                                                                                                          cAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAGGCCTTTCGCCTACTCCACAGCCTCTTCCGTATGACCAACTGCACAATTTTAGCGAGACTATCCAGGCAATCTGCGCCCTGGATGTAGTTGGCGAGc  >  2:162327/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                
CGCTGTTTTCAAAGTTCTCTCGCTTACCAAGGCTTTCACCACGACTCTGGTGCTGGAAGCAATCGACCTCGGCTTATTGTCCCTGACTACCCGAGTTGTTGATGTAATCCCGGAGTTTTGGGGGCGCGACCTGTTCCGCAGTGGCCGGAAAGAGCGTGTCAATGTCGGTCATCTGTTGACCCACCGCGCAGGCCTTTCGCCTACTCCACAGCCTCTTCCGTATGACCAACTGCACAATTTTAGCGAGACTATCCAGGCAATCTGCGCCCTGGATGTAGTTGGCGAGC  >  NC_002947/2210792‑2211078

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