Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,119,836 A→T 8.7% D140E (GAT→GAA)  PP_4506 ← nitroreductase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,119,8360AT 8.7% 152.7 / 4.4 57D140E (GAT→GAA) PP_4506nitroreductase family protein
Reads supporting (aligned to +/- strand):  ref base A (29/23);  new base T (2/3);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

GGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGCGGCGTACCGACATACAGGTAGCCGATCAGCTCCTCGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGGCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAGAATACCGTGCGCCGCACAGCCCGCCGCCAGGCGCTGCTCGGACTTGGGCACCTTGAAGTGGTCCTGCAACCTGGCCACCACCAC  >  NC_002947/5119696‑5119982
                                                                                                                                            |                                                                                                                                                  
ggTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGCGGCGTACCGACATACAGGTAGCCGATCAGCTCCTCGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGGCATCGAAGGc                                                                                                                                            <  1:366538/149‑1 (MQ=255)
 gTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGCGGCGTACCGACATACAGGTAGCCGATCAGCTCCTCGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGGCATCGAAGGcc                                                                                                                                           <  7:118278/149‑1 (MQ=255)
   tACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGCGGCGTACCGACATACAGGTAGCCGATCAGCTCCTCGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGGCATCGAAGGCCAt                                                                                                                                         >  2:109452/1‑149 (MQ=255)
    aCTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGCGGCGTACCGACATACAGGTAGCCGATCAGCTCCTCGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGGCATCGAAGGCCAt                                                                                                                                         <  1:88255/148‑1 (MQ=255)
    aCTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGCGGCGTACCGACATACAGGTAGCCGATCAGCTCCTCGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGGCATCGAAGGCCATg                                                                                                                                        >  7:111661/1‑149 (MQ=255)
     cTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGCGGCGTACCGACATACAGGTAGCCGATCAGCTCCTCGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGGCATCGAAGGCCATGt                                                                                                                                       <  1:244846/149‑1 (MQ=255)
        aCCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGCGGCGTACCGACATACAGGTAGCCGATCAGCTCCTCGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGGCATCGAAGGCCATGTCAc                                                                                                                                    <  3:204641/149‑1 (MQ=255)
        aCCACAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGCGGCGTACCGACATACAGGTAGCCGATCAGCTCCTCGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGCCTTCGATGCCCATGtgc                                                                                                                                     <  3:239269/148‑3 (MQ=255)
                    aaCGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGCGGCGTACCGACATACAGGTAGCCGATCAGCTCCTCGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGGCATCGAAGGCCATGTCACCGGTACGCcaca                                                                                                                        >  3:112778/1‑149 (MQ=255)
                     aCGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGCGGCGTACCGACATACAGGTAGCCGATCAGCTCCTCGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGGCATCGAAGGCCATGTCACCGGTACGCCacac                                                                                                                       >  1:154607/1‑149 (MQ=255)
                     aCGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGCGGCGTACCGACATACAGGTAGCCGATCAGCTCCTCGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGGCATCGAAGGCCATGTCACCGGTACGCCacac                                                                                                                       >  1:368614/1‑149 (MQ=255)
                           tCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGCGGCGTACCGACATACAGGTAGCCGATCAGCTCCTCGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGGCATCGAAGGCCATGTCACCGGTACGCCACACCGCGcc                                                                                                                 >  2:43290/1‑149 (MQ=255)
                                 gTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGCGGCGTACCGACATACAGGTAGCCGATCAGCTCCTCGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGGCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATccc                                                                                                           >  3:361579/1‑149 (MQ=255)
                                 gTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGCGGCGTACCGACATACAGGTAGCCGATCAGCTCCTCGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGCCTTCGATGCCCATGtg                                                                                                                                      >  1:359707/1‑121 (MQ=255)
                                 gTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGCGGCGTACCGACATACAGGTAGCCGATCAGCTCCTCGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGCCTTCGATGCCCATGtg                                                                                                                                      <  2:359707/122‑2 (MQ=255)
                                            cGGCGCCGTCCGCGGTTCGGACAGCGGCGTACCGACATACAGGTAGCCGCACTGCTCCACGTACTCGGCCAACCCCAACCCCTTGTGCACATGGGCATCGAAGGCCAAGTCACCGGTACGCCACACCGCGCCGCTCCCCTGCGCGTgcg                                                                                                >  1:97181/1‑149 (MQ=255)
                                              gcgcCGTCCGCGGTTCGGTCAGCGGCGTACCGACATACAGGTAGCCGATCAGCTCCTCGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGGCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGGGCgca                                                                                              <  5:115507/149‑1 (MQ=255)
                                                    tCCGCGGTTCGGTCAGCGGCGTACCGACATACAGGTAGCCGATCAGCTCCTCGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGGCATCGAAGGCCATGTCACCGGTACGACACACCGCGCCGATCCCCTGCGCGTGCGCAGCAAt                                                                                         >  7:335142/1‑148 (MQ=255)
                                                     ccGCGGTTCGGTCAGCCGCGTACCGTCCTACTGGTAGCCGATCAGCTCCGCGTTGTCGGCCATCCCCAACCCCTTGTGCAGATGGGAATCGATGGCCATGTCAGCGGTACGCCACACCGCGCCGCTCCCCTGCGCGTGCGCAGCAATct                                                                                       <  2:97181/149‑2 (MQ=255)
                                                      cgcgGTTCGGTCAGCGGCGTACCGACATACAGGTAGCCGATCAGCTCCTCGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGGCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAg                                                                                      <  4:133606/149‑1 (MQ=255)
                                                       gcgGTTCGGTCAGCGGCGTACCGACATACAGGTAGCCGATCAGCTCCTCGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGGCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAGa                                                                                     <  1:43290/149‑1 (MQ=255)
                                                         ggTTCGGTCAGCGGCGTACCGACATACAGGTAGCCGATCAGCTCCTCGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGGCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAGAAt                                                                                   <  3:189945/149‑1 (MQ=255)
                                                               gTCAGCGGCGTACCGACATACAGGTAGCCGATCAGCTCCTCGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGGCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAGAATACCGTg                                                                             >  6:42427/1‑149 (MQ=255)
                                                                  aGCGGCGTACCGACATACAGGTAGCCGATCAGCTCCTCGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGGCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAGAATACCGTgcgc                                                                          >  7:72658/1‑149 (MQ=255)
                                                                    cggcgTACCGACATACAGGTAGCCGATCAGCTCCTCGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGGCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAGAATACCGTGcgccg                                                                        <  8:111661/149‑1 (MQ=255)
                                                                            cGACATACAGGTAGCCGATCAGCTCCTCGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGGCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAGAATACCGTGCGCCGCACAGccc                                                                >  6:345888/1‑149 (MQ=255)
                                                                            cGACATACAGGTAGCCGATCAGCTCCTCGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGGCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAGAATACCGTGCGCCGCACAGcc                                                                 <  6:378639/148‑1 (MQ=255)
                                                                            cGACATACAGGTAGCCGATCAGCTCCTCGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGGCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAGAATACCGTGCGCCGCACAGcc                                                                 >  5:378639/1‑148 (MQ=255)
                                                                             gACATACAGGTAGCCGATCAGCTCCTCGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGGCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAGAATACCGTGCGCCGCACAGCccg                                                               >  6:340075/1‑149 (MQ=255)
                                                                                    aGGTAGCCGATCAGCTCCTCGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGCCTTCGATGCCCATGt                                                                                                                                       >  3:327704/1‑70 (MQ=38)
                                                                                    aGGTAGCCGATCAGCTCCTCGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGCCTTCGATGCCCATGt                                                                                                                                       <  4:327704/70‑1 (MQ=38)
                                                                                      gTAGCCGATCAGCTCCTCGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGGCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAGAATACCGTGCGCCGCACAGCCCGCCGCCAGgc                                                      >  5:364476/1‑149 (MQ=255)
                                                                                      gTAGCCGATCAGCTCCTCGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGGCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAGAATACCGTGCGCCGCACAGCCCGCCGCCAGgc                                                      <  8:52191/149‑1 (MQ=255)
                                                                                         gCCGATCAGCTCCTCGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGGCAACGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAGAAAACCGAGCGCCGCACAGCCCCCCGCCCGGCgct                                                   >  8:361164/1‑149 (MQ=255)
                                                                                             aTCAGCTCCTCGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGGCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAGAATACCGTGCGCCGCACAGCCCGCCGCCAGGCGCTGCTc                                               <  2:368614/149‑1 (MQ=255)
                                                                                                aGCTCCTCGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGGCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAGAATACCGTGCGCCGCACAGCCCGCCGCCAGGCGCTGCTCGGa                                            <  7:361164/149‑1 (MQ=255)
                                                                                                   tcctcGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGGCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAGAATACCGTGCGCCGCACAGCCCGCCGCCAGGCGCTGCTCGGACtt                                         >  1:59053/1‑149 (MQ=255)
                                                                                                    cctcGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGGCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAGAATACCGTGCGCCGCACAGCCCGCCGCCAGGCGCTGCTCGGACTTg                                        <  5:42427/149‑1 (MQ=255)
                                                                                                              ggCCAACCCCAAGCCCTTGTGCACATGGGCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAGAATACCGTGCGCCGCACAGCCCGCCGCCAGGCGCTGCTCGGACTTGGGCACCTTg                               <  2:82189/148‑1 (MQ=255)
                                                                                                              ggCCAACCCCAACCCCTTGTGCACATGGGCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAGAATACCGTGCGCCGCACAGCCCGCCGCCAGGCGCTGCTCGGACTTGGGCACCtt                                <  4:112778/147‑1 (MQ=255)
                                                                                                                ccaaccccaaccCCTTGTGCACATGGGCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAGAATACCGTGCGCCGCACAGCCCGCCGCCAGGCGCTGCTCGGACTTGGGCACCTTGAAg                            >  1:261773/1‑149 (MQ=255)
                                                                                                                     cccaaccccTTGTGCACATGGGCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAGAATACCGTGCGCCGCACAGCCCGCCGCCAGGCGCTGCTCGGACTTGGGCACCTTGAAg                            <  5:233358/144‑1 (MQ=255)
                                                                                                                     cccaaccccTTGTGCACATGGGCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAGAATACCGTGCGCCGCACAGCCCGCCGCCAGGCGCTGCTCGGACTTGGGCACATTGAAg                            >  6:233358/1‑144 (MQ=255)
                                                                                                                       caaccccTTGTGCACATGGGCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAGAATACCGTGCGCCGCACAGCCCGCCGCCAGGCGCTGCTCGGACTTGGGCACCTTGAAGTGGTCCt                     <  6:364476/149‑1 (MQ=255)
                                                                                                                           cccTTGTGCACATGGGCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAGAATACCGTGCGCCGCACAGCCCGCCGCCAGGCGCTGCTCGGACTTGGGCACCTTGAAGTGGTCCTGCaa                 >  3:109258/1‑149 (MQ=255)
                                                                                                                                   cacaTGGGCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAGAATACCGTGCGCCGCACAGCCCGCCGCCAGgc                                                      <  6:180425/104‑1 (MQ=255)
                                                                                                                                   cacaTGGGCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAGAATACCGTGCGCCGCACAGCCCGCCGCCAGgc                                                      >  5:180425/1‑104 (MQ=255)
                                                                                                                                    acaTGGGCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAGAATACCGTGCGCCGCACAGCCCGCCGCCAGGCGCTGCTCGGACTTGGGCACCTTGAAGTGGTCCTGCAACCTGGccac        >  5:198604/1‑149 (MQ=255)
                                                                                                                                       tGGGCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAGAATACCGTGCGCCGCACAGCCCGCCGCCAGGCGCTGCTCGGACTTGGGCACCTTGAAGTGGTCCTGCAACCTGGccaccac     >  8:19472/1‑149 (MQ=255)
                                                                                                                                        gggCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAGAATACCGTGCGCCGCACAGCCCGCCGCCAGGCGCTGCTCGGACTTGGGCACCTTGaa                             <  1:236896/124‑1 (MQ=255)
                                                                                                                                        gggCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAGAATACCGTGCGCCGCACAGCCCGCCGCCAGGCGCTGCTCGGACTTGGGCACCTTGaa                             >  2:236896/1‑124 (MQ=255)
                                                                                                                                          gCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAGAATACCGTGCGCCGCACAGCCCGCCGCCAGGCGCTGCTCGGACTTGGGCACCTTGAAGTGGTCCTGCAACCTGGccac        <  8:302533/143‑1 (MQ=255)
                                                                                                                                          gCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAGAATACCGTGCGCCGCACAGCCCGCCGCCAGGCGCTGCTCGGACTTGGGCACCTTGAAGTGGTCCTGCAACCTGGCcaccaccac  >  1:184484/1‑149 (MQ=255)
                                                                                                                                          gCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAGAATACCGTGCGCCGCACAGCCCGCCGCCAGGCGCTGCTCGGACTTGGGCACCTTGAAGTGGTCCTGCAACCTGGCcaccaccac  >  6:303956/1‑149 (MQ=255)
                                                                                                                                          gCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAGAATACCGTGCGCCGCACAGCCCGCCGCCAGGCGCTGCTCGGACTTGGGCACCTTGAAGTGGTCCTGCAACCTGGCcaccaccac  >  2:170959/1‑149 (MQ=255)
                                                                                                                                          gCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAGAATACCGTGCGCCGCACAGCCCGCCGCCAGGCGCTGCTCGGACTTGGGCACCTTGAAGTGGTCCTGCAACCTGGCcaccaccac  >  2:307981/1‑149 (MQ=255)
                                                                                                                                          gCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAGAATACCGTGCGCCGCACAGCCCGCCGCAAGGCGCTGCTCGGACTTGGGCACCTTGAAGTGGTCCTGCAACCTGGccac        >  7:302533/1‑143 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
GGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGCGGCGTACCGACATACAGGTAGCCGATCAGCTCCTCGTTCTCGGCCAACCCCAACCCCTTGTGCACATGGGCATCGAAGGCCATGTCACCGGTACGCCACACCGCGCCGATCCCCTGCGCGTGCGCAGCAATCAGAATACCGTGCGCCGCACAGCCCGCCGCCAGGCGCTGCTCGGACTTGGGCACCTTGAAGTGGTCCTGCAACCTGGCCACCACCAC  >  NC_002947/5119696‑5119982

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