Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,877,824 C→T 8.0% V54V (GTG→GTA)  PP_5151 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,877,8240CT 8.0% 141.6 / 2.7 50V54V (GTG→GTA) PP_5151hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (23/22);  new base T (2/2);  total (25/25)
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.68e-01

GGCCATCCAGCGACAGCTCGATCAGCTCCGAACGGGCGTTGGTGACCGTGAAGAGGGTTTTGCGGTCGGGGTCGTAGGTCAGCGCCGAAATGTCATCATCCAGCCCTTCGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGATTGAACCAGCCACGCTCGAACAGGCGAAACTCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGC  >  NC_002947/5877694‑5877969
                                                                                                                                  |                                                                                                                                                 
ggCCATCCAGCGACAGCTCGATCAGCTCCGAACGGGCGTTGGTGACCGTGAAGAGGGTTTTGCGGTCGGGGTCGTAGGTCAGCGCCGAAATGTCATCATCCAGCCCTTCGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAgg                                                                                                                                 >  6:224533/1‑149 (MQ=255)
 gCCATCCAGCGACAGCTCGATCAGCTCCGAACGGGCGTTGGTGACCGTGAAGAGGGTTTTGCGGTCGGGGTCGTAGGTCAGCGCCGAAATGTCATCATCCAGCCCTTCGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGc                                                                                                                                <  5:344347/149‑1 (MQ=255)
       ctgCGACAGCTCGATCAGCTCCGAACTGGCGTTTGTGACCGTTAAGAGGGTTTTGCGGTCGGGGTCGTAGGTCAGCGCCGAAATGTGATCATCCAGCCCTTCGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAAtg                                                                                                                           <  7:171829/146‑1 (MQ=255)
              aGCTCGATCAGCTCCGAACGGGCGTTGGTGACCGTGAAGAGGGTTTTGCGGTCGGGGTCGTAGGTCAGCGCCGAAATGTCATCATCCAGCCCTTCGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTc                                                                                                                   <  7:186046/149‑1 (MQ=255)
              aGCTCGATCAGCTCCGAACGGGCGTTGGTGACCGTGAAGAGGGTTTTGCGGTCGGGGTCGTAGGTCAGCGCCGAAATGTCATCATCCAGCCCTTCGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTc                                                                                                                   <  4:101486/149‑1 (MQ=255)
                 tCGATCAGCTCCGAACGGGCGTTGGTGACCGTGAAGAGGGTTTTGCGGTCGGGGTCGTAGGTCAGCGCCGAAATGTCATCATCCAGCCCTTCGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCca                                                                                                                 >  4:305794/1‑148 (MQ=255)
                   gATCAGCTCCGAACGGGCGTTGGTGACCGTGAAGAGGGTTTTGCGGTCGGGGTCGTAGGTCAGCGCCGAAATGTCATCATCCAGCCCTTCGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCcaca                                                                                                               <  4:310195/148‑1 (MQ=255)
                     tCAGCTCCGAACGGGCGTTGGTGACCGTGAAGAGGGTTTTGCGGTCGGGGTCGTAGGTCAGCGCCGAAATGTCATCATCCAGCCCTTCGATCGGTTGCGCCTCCAGCACTACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGc                                                                                                            >  1:335409/1‑149 (MQ=255)
                     tCAGCTCCGAACGGGCGTTGGTGACCGTGAAGAGGGTTTTGCGGTCGGGGTCGTAGGTCAGCGCCGAAATGTCATCATCCAGCCCTTCGATCGGTTGCGCCTCCAGCACTACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGc                                                                                                            >  1:340815/1‑149 (MQ=255)
                       aGCTCCGAACGGGCGTTCGTGACCGTGAAGAGGGTTTTGCGGTCCGGGTCGTAGGTCAGCGCCGAAATGTCAGCATCCAGCCCTTCGCTCGGGTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCt                                                                                                           <  2:253671/148‑1 (MQ=255)
                             gAACGGGCGTTGGTGACCGTGAAGAGGGTTTTGCGGTCGGGGTCGTAGGTCAGCGCCGAAATGTCATCATCCAGCCCTTCGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGccacca                                                                                                    <  3:305794/149‑1 (MQ=255)
                                   gCGTTGGTGACCGTGAAGAGGGTTTTGCGGTCGGGGTCGTAGGTCAGCGCCGAAATGTCATCATCCAGCCCTTCGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCtt                                                                                              <  2:165558/149‑1 (MQ=255)
                                    cGTTGGTGACCGTGAAGAGGGTTTTGCGGTCGGGGTCGTAGGTCAGCGCCGAAATGTCATCATCCAGCCCTTCGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTg                                                                                             <  7:244200/149‑1 (MQ=255)
                                            aCCGTGAAGAGGGTTTTGCGGTCGGGGTCGTAGGTCAGCGCCGAAATGTCATCATCCAGCCCTTCGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGATTGaa                                                                                     <  3:261329/149‑1 (MQ=255)
                                               gTGAAGAGGGTTTTGCGGTCGGGGTCGTAGGTCAGCGCCGAAATGTCATCATCCAGCCCTTCGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGATTGAACCa                                                                                  <  3:360277/149‑1 (MQ=255)
                                                  aaGAGGGTTTTGCGGTCGGGGTCGTAGGTCAGCGCCGAAATGTCATCATCCAGCCCTTCGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGATTGAACCAGc                                                                                <  5:224533/148‑1 (MQ=255)
                                                    gagGGTTTTGCGGTCGGGGTCGTAGGTCAGCGCCGAAATGTCATCATCCAGCCCTTCGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGATTGAACCAGCAAc                                                                             >  3:30163/1‑149 (MQ=255)
                                                       ggTTTTGCGGTCGGGGTCGTAGGTCAGCGCCGAAATGTCATCATCCAGCCCTTCGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGATTGAACCAGCCACGCt                                                                          >  5:286261/1‑149 (MQ=255)
                                                            tGCGGTCGGGGTCGTAGGTCAGCGCCGAAATGTCATCATCCAGCCCTTCGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGATTGAACCAGCCACGCTCGaa                                                                      <  5:288139/148‑1 (MQ=255)
                                                               ggTCGGGGTCGTAGGTCAGCGCCGAAATGTCATCATCCAGCCCTTCGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGATTGAACCAGCCACGCTCGAACAgg                                                                  <  3:328006/149‑1 (MQ=255)
                                                                 tCGGGGTCGTAGGTCAGCGCCGAAATGTCATCATCCAGCCCTTCGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGATTGAACCAGCCACGCTCGAACAGGc                                                                 >  2:14348/1‑148 (MQ=255)
                                                                     ggTCGTAGGTCAGCGCCGAAATGTCATCATCCAGCCCTTCGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGATTGAACCAGCCACGCTCGAACAGGCGAAAc                                                            >  8:90442/1‑149 (MQ=255)
                                                                            ggTCAGCGCCGAAATGTCATCATCCAGCCCTTCGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGATTGAACCAGCCACGCTCGAACAGGCGAAACTCCTGc                                                      <  2:254657/148‑1 (MQ=255)
                                                                             gTCAGCGCCGAAATGTCATCATCCAGCCCTTCGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGAtt                                                                                        >  3:227813/1‑113 (MQ=255)
                                                                             gTCAGCGCCGAAATGTCATCATCCAGCCCTTCGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGAtt                                                                                        <  4:227813/113‑1 (MQ=255)
                                                                                              atcatcCAGCCCTTCGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGATTGAACCAGCCACGCTCGAACAGGCGAAACTCCTGCCCGGCCCAGcc                                           >  6:2226/1‑141 (MQ=255)
                                                                                              atcatcCAGCCCTTCGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGATTGAACCAGCCACGCTCGAACAGGCGAAACTCCTGCCCGGCCCAGcc                                           <  5:2226/141‑1 (MQ=255)
                                                                                              atcatcCAGCCCTTCGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGATTGAACCAGCCACGCTCGAACAGGCGAAACTCCTGCCCGGCCCAGCCAagcag                                     >  1:184335/1‑147 (MQ=255)
                                                                                               tcatcCAGCCCTTCGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGATTGAACCAGCCACGCTCGAACAGGCGAAACTCCTGCCCGGCCCAGCCAagcagcag                                  >  4:152386/1‑149 (MQ=255)
                                                                                               tcatcCAGCCCTTCGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGATTGAACCAGCCACGCTCGAACAGGCGAAACTCCTGCCCGGCCCAGCCAagcagcag                                  <  6:286261/149‑1 (MQ=255)
                                                                                                catcCAGCCCTTCGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGATTGAACCAGCCACGCTCGAACAGGCGAAACTCCTGCCCg                                                   >  4:241900/1‑131 (MQ=255)
                                                                                                catcCAGCCCTTCGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGATTGAACCAGCCACGCTCGAACAGGCGAAACTCCTGCCCg                                                   <  3:241900/131‑1 (MQ=255)
                                                                                                  tcCAGCCCTTCGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGATTGAACCAGCCACGCTCGAACAGGCGAAACTCCTGCCCGGCCCAGCCAagcagcag                                  >  8:303426/1‑146 (MQ=255)
                                                                                                  tcCAGCCCTTCGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGATTGAACCAGCCACGCTCGAACAGGCGAAACTCCTGCCCGGCCCAGCCAagcagcag                                  >  3:225500/1‑146 (MQ=255)
                                                                                                   ccAGCCCTTCGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGATTGAACCAGCCACGCTCGAACAGGCGAAACTCCTGCCCGGCCCAGCCAAgcagcagca                                >  4:102008/1‑147 (MQ=255)
                                                                                                       cccTTCGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGATTGAACCAGCCACGCTCGAACAGGCGAAACTCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGc                            >  1:95286/1‑147 (MQ=255)
                                                                                                        ccTTCGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGATTGAACCAGCCACGCTCGAACAGGCGAAACTCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAg                         <  3:102008/149‑1 (MQ=255)
                                                                                                          ttCGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGATTGAACCAGCCACGCTCGAACAGGCGAAACTCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAggg                       >  2:207982/1‑149 (MQ=255)
                                                                                                            cGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGATTGAACCAGCCACGCTCGAACAGGCGAAACTCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAggg                       >  4:326180/1‑147 (MQ=255)
                                                                                                            cGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGATTGAACCAGCCACGCTCGAACAGGCGAAACTCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGc                      <  7:90442/148‑1 (MQ=255)
                                                                                                                  gTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGATTGAACCAGCCACGCTCGAACAGGCGAAACTCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAg                >  1:253074/1‑148 (MQ=255)
                                                                                                                     gcgcCTCCAGCACTACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGATTGAACCAGCCACGCTCGAACAGGCGAAACTCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACa            <  2:340815/149‑1 (MQ=255)
                                                                                                                      cgcCTCCAGCACTACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGATTGAACCAGCCACGCTCGAACAGGCGAAACTCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACa            <  2:335409/148‑1 (MQ=255)
                                                                                                                       gcCTCCAGCACCACCCGGTAGCCATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGATTGAACCAGCCACGCTCGAACAGGCGAAACTCCTGCCCGGCCCAGCCAAGCAGCCGCAGAGCAAGGGCAGCCAGAACAAg          >  1:310706/1‑149 (MQ=255)
                                                                                                                        ccTCCAGCACACCCCGGTATCGATCCAGGCCAATGCTGCGCTACACAGGCTGCCACCAGGTCTTGAGATTGAACCAGCCACGCTCGAACAGGCGAAACTCCTGCCCTGCCCAGCCAAGCAGCAGCAGAGCAATGGCAGCCAGAACAAgg         <  7:92623/149‑1 (MQ=255)
                                                                                                                         cTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGATTGAACCAGCCACGCTCGAACAGGCGAAACTCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAgg         >  6:202517/1‑148 (MQ=255)
                                                                                                                          tCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGATTGAACCAGCCACGCTCGAACAGGCGAAACTCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGaa       >  1:9587/1‑149 (MQ=255)
                                                                                                                            cAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGATTGAACCAGCCACGCTCGAACAGGCGAAACTCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGaa       >  4:168255/1‑147 (MQ=255)
                                                                                                                            cAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGATTGAACCAGCCACGCTCGAACAGGCGAAACTCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGaa       <  1:14348/147‑1 (MQ=255)
                                                                                                                              gCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGATTGAACCAGCCACGCTCGAACAGGCGAAACTCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGAAGCCAGAACAAGGAACAgg   >  3:222859/1‑149 (MQ=255)
                                                                                                                               caccacCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGATTGAACCAGCCACGCTCGAACAGGCGAAACTCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGc  >  8:145650/1‑149 (MQ=255)
                                                                                                                                  |                                                                                                                                                 
GGCCATCCAGCGACAGCTCGATCAGCTCCGAACGGGCGTTGGTGACCGTGAAGAGGGTTTTGCGGTCGGGGTCGTAGGTCAGCGCCGAAATGTCATCATCCAGCCCTTCGATCGGTTGCGCCTCCAGCACCACCCGGTAGCGATCCAGGCCAATGCTGCGCTCCACAGGCTGCCACCAGGTCTTGAGATTGAACCAGCCACGCTCGAACAGGCGAAACTCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGC  >  NC_002947/5877694‑5877969

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