Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 222,817 C→T 8.0% Y224Y (TAC→TAT)  PP_0170 → ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947222,8170CT 8.0% 143.5 / 2.4 50Y224Y (TAC→TAT) PP_0170ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (25/21);  new base T (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.00e-01

CGACAAATACACCTTGGCCGTGCCGCACCCGGACGCCACCGCCGCACTGTTGGCGGGCGGCACCGAGCTGAACGGGCACTTCTCCAACCCGCCGTTCCAGGACCAGGTGCTGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTACGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTACTGAGAAGTTCCGCAAGGACAACCCCAAGACCTACAAGGCCTTCGTCGATGCCCTGGCCGAAGCCGCAGACTTCGCCCAGAAAGACAAGGCAGCTGCCGCCGACACC  >  NC_002947/222670‑222964
                                                                                                                                                   |                                                                                                                                                   
cGACAAATACACCTTGGCCGTGCCGCACCCGGACGCCACCGCCGCACTGTTGGCGGGCGGCACCGAGCTGAACGGGCACTTCTCCAACCCGCCGTTCCAGGACCAGGTGCTGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTAc                                                                                                                                                     >  1:72934/1‑148 (MQ=255)
 gACAAATACACCTTGGCCGTGCCGCACCCGGACGCCACCGCCGCACTGTTGGCGGGCGGCACCGAGCTGAACGGGCACTTCTCCAACCCGCCGTTCCAGGACCAGGTGCTGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTACGa                                                                                                                                                   >  3:64490/1‑149 (MQ=255)
   cAAATACACCTTGGCCGTGCCGCACCCGGACGCCACCGCCGCACTGTTGGCGGGCGGCACCGAGCTGAACGGGCACTTCTCCAACCCGCCGTTCCAGGACCAGGTGCTGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTACGAcc                                                                                                                                                 <  2:72934/149‑1 (MQ=255)
      aTACACCTTGGCCGTGCCGCCCCCGGACGCCACCGCCGCACTGTTGGCGGGCGGCACCGAGCTGAACGGGCACTTCTCCAACCCGCCGTTCCAGGACCAGGTGCTGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTACGACCTa                                                                                                                                               >  8:219209/1‑148 (MQ=255)
      aTACACCTTGGCCGTGCCGCACCCGGACGCCACCGCCGCACTGTTGGCGGGCGGCACCGAGCTGAACGGGCACTTCTCCAACCCGCCGTTCCAGGACCAGGTGCTGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTACGACCTAc                                                                                                                                              >  6:223979/1‑149 (MQ=255)
      aTACACCTTGGCCGTGCCGCACCCGGACGCCACCGCCGCACTGTTGGCGGGCGGCACCGAGCTGAACGGGCACTTCTCCAACCCGCCGTTCCAGGACCAGGTGCTGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTACGACCTAc                                                                                                                                              <  7:79892/149‑1 (MQ=255)
            cTTGGCCGTGCCGCACCCGGACGCCACCGCCGCACTGTTGGCGGGCGGCACCGAGCTGAACGGGCACTTCTCCAACCCGCCGTTCCAGGACCAGGTGCTGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTAc                                                                                                                                                     <  2:50433/136‑1 (MQ=255)
            cTTGGCCGTGCCGCACCCGGACGCCACCGCCGCACTGTTGGCGGGCGGCACCGAGCTGAACGGGCACTTCTCCAACCCGCCGTTCCAGGACCAGGTGCTGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTAc                                                                                                                                                     >  1:50433/1‑136 (MQ=255)
                gCCGTGCCGCACCCGGACGCCACCGCCGCACTGTTGGCGGGCGGCACCGAGCTGAACGGGCACTTCTCCAACCCGCCGTTCCAGGACCAGGTGCTGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTACGACCTACTCGGCCCCaa                                                                                                                                    >  1:293738/1‑149 (MQ=255)
                          aCCCGGACGCCACCGCCGCACTGTTGGCGGGCGGCACCGAGCTGAACGGGCACTTCTCCAACCCGCCGTTCCAGGACCAGGTGCTGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTACGACCTACTCGGCCCCAACTCGCCGAcc                                                                                                                          <  4:64490/149‑1 (MQ=255)
                            ccGGACGCCACCGCCGCACTGTTGGCGGGCGGCACCGAGCTGAACGGGCACTTCTCCAACCCGCCGTTCCAGGACCAGGTGCTGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTACGACCTACTCGGCCCCAACTCGCCGACCtt                                                                                                                        >  1:149066/1‑149 (MQ=255)
                                    cACCGCCGCACTGTTGGCGGGCGGCACCGAGCTGAACGGGCACTTCTCCAACCCGCCGTTCCAGGACCAGGTGCTGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTACGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTg                                                                                                                >  8:26497/1‑149 (MQ=255)
                                     aCCGCCGCACTGTTGGCGGGCGGCACCGAGCTGAACGGGCACTTCTCCAACCCGCCGTTCCAGGACCAGGTGCTGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTATGACCTa                                                                                                                                               <  2:66072/117‑1 (MQ=255)
                                     aCCGCCGCACTGTTGGCGGGCGGCACCGAGCTGAACGGGCACTTCTCCAACCCGCCGTTCCAGGACCAGGTGCTGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTATGACCTa                                                                                                                                               >  1:66072/1‑117 (MQ=255)
                                       cgccgcACTGTTGGCGGGCGGCACCGAGCTGAACGGGCACTTCTCCAACCCGCCGTTCCAGGACCAGGTGCTGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTACGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTa                                                                                                             >  4:188547/1‑149 (MQ=255)
                                         ccgcACTGTTGGCGGGCGGCACCGAGCTGAACGGGCACTTCTCCAACCCGCCGTTCCAGGACCAGGTGCTGGCGAACAATGATGTGCACGTGGTGCTCAACAGCTACGACCTACTCGGCCACAACTCGCCGACCTTGCTGTTTGCTACt                                                                                                           >  6:28105/1‑149 (MQ=255)
                                         ccgcACTGTTGGCGGGCGGCACCGAGCTGAACGGGCACTTCTCCAACCCGCCGTTCCAGGACCAGGTGCTGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTACGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTACt                                                                                                           <  3:188547/149‑1 (MQ=255)
                                           gcACTGTTGGCGGGCGGCACCGAGCTGAACGGGCACTTCTCCAACCCGCCGTTCCAGGACCAGGTGCTGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTACGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTACTga                                                                                                         <  4:44546/149‑1 (MQ=255)
                                             aCTGTTGGCGGGCGGCACCGAGCTGAACGGGCACTTCTCCAACCCGCCGTTCCAGGACCAGGTGCTGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTACGACCTACTCGGCCCCAACTCGCCGACCTTg                                                                                                                       <  8:148888/133‑1 (MQ=255)
                                             aCTGTTGGCGGGCGGCACCGAGCTGAACGGGCACTTCTCCAACCCGCCGTTCCAGGACCAGGTGCTGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTACGACCTACTCGGCCCCAACTCGCCGACCTTg                                                                                                                       >  7:148888/1‑133 (MQ=255)
                                              cTGTTGGCGGGCGGCACCGAGCTGAACGGGCACTTCTCCAACCCGCCGTTCCAGGACCAGGTGCTGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTACGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTACTGAGaa                                                                                                      <  5:223979/149‑1 (MQ=255)
                                                 ttGGCGGGCGGCACCGAGCTGAACGGGCACTTCTCCAACCCGCCGTTCCAGGACCAGGTGCTGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTACGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTACTGAGAAGtt                                                                                                   <  2:293738/149‑1 (MQ=255)
                                                            cACCGAGCTGAACGGGCACTTCTCCAACCCGCCGTTCCAGGACCAGGTGCTGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTACGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTACTGAGAAGTTCCGCAAGGAc                                                                                         >  6:186725/1‑148 (MQ=255)
                                                             aCCGAGCTGAACGGGCACTTCTCCAACCCGCCGTTCCAGGACCAGGTGCTGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTACGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTACTGAGAAGTTCCGCAAGGACaa                                                                                       >  3:71763/1‑149 (MQ=255)
                                                                           gCACTTCTCCAACCCGCCGTTCCAGGACCAGGTGCTGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTACGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTACTGAGAAGTTCCGCAAGGACAACCCCAAGACCTAc                                                                          >  2:324797/1‑148 (MQ=255)
                                                                            cACTTCTCCAACCCGCCGTTCCAGGACCAGGTGCTGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTACGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTACTGAGAAGTTCCGCAAGGACAACCCCAAGACCTACaa                                                                        <  4:71763/149‑1 (MQ=255)
                                                                             aCTTCTCCAACCCGCCGTTCCAGGACCAGGTGCTGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTACGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTACTGAGAAGTTCCGCAAGGACAACCCCAAGACCTACAAg                                                                       >  2:182573/1‑149 (MQ=255)
                                                                               ttCTCCAACCCGCCGTTCCAGGACCAGGTGCTGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTACGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTACTGAGAAGTTCCGCAAGGACAACCCCAAGACCTACAAGGc                                                                     >  1:67018/1‑149 (MQ=255)
                                                                                 ctcCAACCCGCCGTTCCAGGACCAGGTGCTGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTACGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTACTGAGAAGTTCCGCAAGGACAACCCCAAGACCTACAAg                                                                       >  7:267015/1‑145 (MQ=255)
                                                                                 ctcCAACCCGCCGTTCCAGGACCAGGTGCTGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTACGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTACTGAGAAGTTCCGCAAGGACAACCCCAAGACCTACAAg                                                                       <  8:267015/145‑1 (MQ=255)
                                                                                          gccgTTCCAGGACCAGGTGCTGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTACGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTACTGAGAAGTTCCGCAAGGACAACCCCAAGACATAGAAGGCATtcggagctg                                                          >  6:251864/1‑143 (MQ=255)
                                                                                              ttCCAGGACCAGGTGCTGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTACGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTACTGAGAAGTTCCGCAAGGACAACCCCAAGACCTACAAGGCCTtcgtcgct                                                           <  6:57994/144‑3 (MQ=255)
                                                                                              ttCCAGGACCAGGTGCTGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTACGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTACTGAGAAGTTCCGCAAGGACAACCCCAAGACCTACAAGGCCTtcgtcgct                                                           >  5:57994/1‑142 (MQ=255)
                                                                                                 cAGGACCAGGTGCTGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTACGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTACTGAGAAGTTCCGCAAGGACAACCCCAAGACCTACAAGGCCTTCGTCGATGCCCTgg                                                    <  5:251864/148‑1 (MQ=255)
                                                                                                      ccAGGTGCTGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTACGCCCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTACTGAGAAGTTCCGCAAGGACAACCCCAAGACCTACAAGGCCTTCGTCGATGCCCTGGCCGAAg                                              >  5:63328/1‑149 (MQ=255)
                                                                                                           tGCTGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTACGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTACTGAGAAGTTCCGCAAGGACAACCCCAAGACCTACAAGGCCTTCGTCGATGCCCTGGCCGAAGCCGCa                                         <  7:219209/149‑1 (MQ=255)
                                                                                                              tGGCCAAGAAAGATGTGCACGTGGTGCTCAACAGCTACGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTACTGAGAAGTTCCGCAAGGACAACCCCAAGACCTACAAGGCCTTCGTCGAt                                                           <  7:244959/128‑1 (MQ=255)
                                                                                                              tGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTACGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTACTGAGAAGTTCCGCAAGGACAACCCCAAGACCTACAAGGCCTTCGTCGATGCCCTGGCCGAAGCCGCAGAc                                      >  7:123127/1‑149 (MQ=255)
                                                                                                                  caacaaAGATGTGCACGTGGTGCTCAACAGCTACGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTACTGAGAAGTTCCGCAAGGACAACCCCAAGACCTACAAGGCCTTCGTCGATGCCCTGGCCGAAGCCGCAGACTTCg                                  >  6:174449/1‑149 (MQ=255)
                                                                                                                   aacaaAGATGTGCACGTGGTGCTCAACAGCTACGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTACTGAGAAGTTCCGCAAGGACAACCCCAAGACCTACAAGGCCTTCGTCGATGCCCTGGCCGAAGCCGCAGACTTCGc                                 <  1:182573/149‑1 (MQ=255)
                                                                                                                   aacaaAGATGTGCACGTGGTGCTCAACAGCTACGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTACTGAGAAGTTCCGCAAGGACAACCCCAAGACCTACAAGGCCTTCGTCGATGCCCTGGCCGAAGCCGCAGACTTCGc                                 >  3:119583/1‑149 (MQ=255)
                                                                                                                           tgtgCACGTGGTGCTCAACAGCTACGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTACTGAGAAGTTCCGCAAGGACAACCCCAAGACCTACAAGGCCTTCGTCGATGCCCTGGCCGAAGCCGCAGACTTCGCCCAGaaa                          <  8:123127/148‑1 (MQ=255)
                                                                                                                            gtgCACGTGGTGCTCAACAGCTACGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTACTGAGAAGTTCCGCAAGGACAACCCCAAGACCTACAAGGCCTTCGTCGATGCCCTGGCCGAAGc                                             >  6:317979/1‑128 (MQ=255)
                                                                                                                               cACGTGGTGCTCAACAGCTACGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTACTGAGAAGTTCCGCAAGGACAACCCCAAGACCTACAAGGCCTTCGTCGATGCCCTGGCCGAAGCCGCAGACTTCGCCCAGAAAGACaa                     <  5:186725/149‑1 (MQ=255)
                                                                                                                               cACGTGGTGCTCAACAGCTACGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTACTGAGAAGTTCCGCAAGGACAACCCCAAGACCTACAAGGCCTTCGTCGATGCCCTGGCCGAAGCCGCAGACTTCGCCCAGAAAGACaa                     <  5:28105/149‑1 (MQ=255)
                                                                                                                               cACGTGGTGCTCAACAGCTACGACATACTAGGCCCCAAATCGCCGACCTTGCTGTTTGCTACTGAGAAGTTCCGCAAGGACAACCCCAAGACCTACAAGTCCTTCGTCGATGCCCTGGCCGAAGCCGCAGACTTCGCCCAGAAAGACaa                     >  1:187894/1‑149 (MQ=255)
                                                                                                                                 cGTGGTGCTCAACAGCTACGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTACTGAGAAGTTCCGCAAGGACAACCCCAAGACCTACAAGGCCTTCGTCGATGCCCTGGCCGAAGCCGCAGACTTCGCCCAGAAAGACAAgg                   >  7:197267/1‑149 (MQ=255)
                                                                                                                                  gtggtgCTCAACAGCTACGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTACTGAGAAGTTCCGCAAGGACAACCCCAAGACCTACAAGGCCTTCGTCGATGCCCTGGCCGAAGCCGCAGACTTCGCCCAGAAAGACAAGGc                  <  4:119583/149‑1 (MQ=255)
                                                                                                                                   tggtgCTCAACAGCTACGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTACTGAGAAGTTCCGCAAGGACAACCCCAAGACCTACAAGGCCTTCGTCGATGCCCTGGCCGAAGCCGCAGACTTCGCCCAGAAAGACAAGGCa                 <  8:197267/149‑1 (MQ=255)
                                                                                                                                           aaCAGCTATGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTACTGAGAAGTTCCGCAAGGACAACCCCAAGACCTACAAGGCCTTCGTCGATGCCCTGGCCGAAGCCGCAGACTTCGCCCAGAAAg                         >  4:319415/1‑133 (MQ=255)
                                                                                                                                           aaCAGCTATGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTACTGAGAAGTTCCGCAAGGACAACCCCAAGACCTACAAGGCCTTCGTCGATGCCCTGGCCGAAGCCGCAGACTTCGCCCAGAAAg                         <  3:319415/133‑1 (MQ=255)
                                                                                                                                           aaCAGCTACGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTACTGAGAAGTTCCGc                                                                                               >  5:263100/1‑63 (MQ=255)
                                                                                                                                           aaCAGCTACGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTACTGAGAAGTTCCGc                                                                                               <  6:263100/63‑1 (MQ=255)
                                                                                                                                                   cGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTACTGAGAAGTTCCGCAAGGACAACCCCAAGACCTACAAGGCCTTCGTCGATGCCCTGGCCGAAGCCGCAGACTTCGCCCAGAAAGACAAGGCAACTGCCGCCGACAcc  <  6:63328/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CGACAAATACACCTTGGCCGTGCCGCACCCGGACGCCACCGCCGCACTGTTGGCGGGCGGCACCGAGCTGAACGGGCACTTCTCCAACCCGCCGTTCCAGGACCAGGTGCTGGCCAACAAAGATGTGCACGTGGTGCTCAACAGCTACGACCTACTCGGCCCCAACTCGCCGACCTTGCTGTTTGCTACTGAGAAGTTCCGCAAGGACAACCCCAAGACCTACAAGGCCTTCGTCGATGCCCTGGCCGAAGCCGCAGACTTCGCCCAGAAAGACAAGGCAGCTGCCGCCGACACC  >  NC_002947/222670‑222964

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