Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,260,386 G→A 8.1% S792S (AGC→AGT)  PP_3734 ← ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,260,3860GA 8.1% 136.7 / 2.1 49S792S (AGC→AGT) PP_3734ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (19/26);  new base A (2/2);  total (21/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.54e-01

TTTGCGGTGCAGGCTGATCACAGCCGCACAACAACAACCCGCTCAGCAACAGCCACGCACTAAGCTTCATCGGACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGTGAGCATGCCCAGCACGGTCAACTGCAGCAGTTGCACGGGGAACACATACAGCGGTAAACGCCAGCCGAAGGCCTGCACGTTCACCACTGCGACCAGGCACCAGGCCAGCAACAGCCCAAGCGGAATCGCCAGCAGTACGGT  >  NC_002947/4260243‑4260532
                                                                                                                                               |                                                                                                                                                  
tttGCGGTGCAGGCTGATCACAGCCGCACAACAACAACCCGCTCAGCAACAGCCACGCACTAAGCTTCATCGGACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGt                                                                                                                                               >  6:325032/1‑149 (MQ=255)
   gCGGTGCAGGCTGATCACAGCCGCACAACAACAACCCGCTCAGCAACAGCCACGCACTAAGCTTCATCGGACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGTGa                                                                                                                                             <  8:239848/148‑1 (MQ=255)
    cGGTGCAGGCTGATCACAGCCGCACAACAACAACCCGCTCAGCACCAGCCACGCACTAAGCTTCATCGGACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGTGAGc                                                                                                                                           <  7:302864/149‑1 (MQ=255)
    cGGTGCAGGCTGATCACAGCCGCACAACAACAACCCGCTCAGCAACAGCCACGCACTAAGCTTCATCGGACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCCTGCTggcctgc                                                                                                                                           <  8:284816/149‑6 (MQ=255)
       tGCAGGCTGATCACAGCCGCACAACAACAACCCGCTCAGCAACAGCCACGCACTAAGCTTCATCGGACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGTGAGCATg                                                                                                                                        >  2:78188/1‑149 (MQ=255)
             cTGATCACAGCCGCACAACAACAACCCGCTCAGCAACAGCCACGCACTAAGCTTCATCGGACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGTGAGCATGCCCAGc                                                                                                                                  <  1:208104/149‑1 (MQ=255)
                aTCACAGCCGCACAACAACAACCCGCTCAGCAACAGCCACGCACTAAGCTTCATCGGACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGACCAGCAGGCTGGTGAGCATGCCCAGCACg                                                                                                                               <  6:325109/149‑1 (MQ=255)
                 tCACAGCCGCACAACAACAACCCGCTCAGCAACAGCCACGCACTAAGCTTCATCGGACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGTGAGCATGCCCAGCACgg                                                                                                                              >  2:125780/1‑149 (MQ=255)
                 tCACAGCCGCACAACAACAACCCGCTCAGCAACAGCCACGCACTAAGCTTCATCGGACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGACTGGTGAGCATGCCCAGCACgg                                                                                                                              >  1:83660/1‑149 (MQ=255)
                 tCACAGCCGCACAACAACAACCCGCTCAGCAACAGCCACGCACTAAGCTTCATCGGACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGACTGGTGAGCATGCCCAGCACg                                                                                                                               >  3:65748/1‑148 (MQ=255)
                    caGCCGCACAACAACAACCCGCTCAGCAACAGCCACGCACTAAGCTTCATCGGACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGTGAGCATGCCCAGCACg                                                                                                                               <  7:87947/145‑1 (MQ=255)
                    caGCCGCACAACAACAACCCGCTCAGCAACAGCCACGCACTAAGCTTCATCGGACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGTGAGCATGCCCAGCACg                                                                                                                               >  8:87947/1‑145 (MQ=255)
                     aGCCGCACAACAACAACCCGCTCAGCAACAGCCACGCACTAAGCTTCATCGGACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGTGAGCATGCCCAGCACGGTCaa                                                                                                                          >  5:239851/1‑149 (MQ=255)
                        cGCACAACAACAACCCGCTCAGCAACAGCCACGCACTAAGCTTCATCGGACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGTGAGCATGCCCAGCACGGTCAACTg                                                                                                                       <  1:78188/149‑1 (MQ=255)
                             aacaacaacCCGCTCAGCAACAGCCACGCACTAAGCTTCATCGGACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGt                                                                                                                                               <  8:205464/120‑1 (MQ=255)
                             aacaacaacCCGCTCAGCAACAGCCACGCACTAAGCTTCATCGGACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGt                                                                                                                                               >  7:205464/1‑120 (MQ=255)
                                aacaacCCGCTCAGCAACAGCCACGCACTAAGCTTCATCGGACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGTGAGCATGCCCAGCACGGTCAACTGCAGCAGtt                                                                                                               <  7:110655/149‑1 (MQ=255)
                                   aacCCGCTCAGCAACAGCCACGCTCTAAGCTTCATCGGACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGTGAGCATGCCCAGCACGGTCAACTGCAGCAGTTGCa                                                                                                            <  7:63742/149‑1 (MQ=255)
                                   aacCCGCTCAGCAACAGCCACGCACTAAGCTTCATCGGACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGTGAGCATGCCCAGCACGGTCAACTGCAGCAGTTGCa                                                                                                            >  3:308278/1‑149 (MQ=255)
                                    acCCGCTCAGCAACAGCCACGCACTAAGCTTCATCGGACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGTGAGCATGCCCAGCACGGTCAACTGCAGCAGTTGCAc                                                                                                           <  4:308278/149‑1 (MQ=255)
                                     cccGCTCAGCAACAGCCACGCACTAAGCTTCATCGGACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGTGAGCATGCCCAGCACGGTCAACTGCAGCAGTTGCACg                                                                                                          >  6:3099/1‑149 (MQ=255)
                                      ccGCTCAGCAACAGCCACGCACTAAGCTTCATCGGACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAgg                                                                                                                                                    <  7:279935/106‑1 (MQ=255)
                                      ccGCTCAGCAACAGCCACGCACTAAGCTTCATCGGACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAgg                                                                                                                                                    >  8:279935/1‑106 (MQ=255)
                                         cTCAGCAACAGCCACGCACTAAGCTTCATCGGACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGTGAGCATGCCCAGCACGGTCAACTGCAGCAGTTGCACGGGGa                                                                                                      >  7:91120/1‑149 (MQ=255)
                                          tCAGCAACAGCCACGCACTAAGCTTCATCGGACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGTGAGCATGCCCAGCACGGTCAACTACAGCAGTTGCACGGGGaa                                                                                                     >  3:236770/1‑149 (MQ=255)
                                             gCAACAGCCACGCACTAAGCTTCATCGGACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGTGAGCATGCCCAGCACGGTCAACTGCAGCAGTTGCACGGGGAacac                                                                                                  <  5:102132/149‑1 (MQ=255)
                                              cAACAGCCACGCACTAAGCTTCATCGGACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGTGAGCATGCCCAGCACGGTCAACTGCAGCAGTTGCACGGGGAAcaca                                                                                                 <  2:204740/149‑1 (MQ=255)
                                                        gCACTAAGCTTCATCGGACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGTGAGCATGCCCAGCACGGTCAACTGCAGCAGTTGCACGGGGAACACATACAg                                                                                            >  7:18377/1‑144 (MQ=255)
                                                        gCACTAAGCTTCATCGGACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGTGAGCATGCCCAGCACGGTCAACTGCAGCAGTTGCACGGGGAACACATACAg                                                                                            <  8:18377/144‑1 (MQ=255)
                                                          aCTAAGCTTCATCGGACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGTGAGCATGCCCAGCACGGTCAACTGCAGCAGTTGCACGGGGAACACATACAGCGGTaaa                                                                                     <  1:125780/149‑1 (MQ=255)
                                                            tAAGCTTCATCGGACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGTGAGCATGCCCAGCACGGTCAACTGCAGCAGTTGCACGGGGAACACATACAGCGGTAAACg                                                                                   >  7:271920/1‑149 (MQ=255)
                                                                   cATCGGACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGTGAGCATGCCCAGCACGGTCAACTGCAGCAGTTGCACGGGGAACACATACAGCGGTAAACGCCAGc                                                                              <  5:325032/147‑1 (MQ=255)
                                                                    aTCGGACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGTGAGCATGCCCAGCACGGTCAACTGCAGCAGTTGCACGGGGAACACATACAGCGGTAAACGCCAGCCg                                                                            >  4:154142/1‑148 (MQ=255)
                                                                    aTCGGACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGACTGGTGAGCATGCCCAGCACGGTCAACTGCAGCAGTTGCACGGGGAACACATACAGCGGTAAACGCCAGCCg                                                                            <  4:65748/148‑1 (MQ=255)
                                                                    aTCGGACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGACTGGTGAGCATGCCCAGCACGGTCAACTGCAGCAGTTGCACGGGGAACACATACAGCGGTAAACGCCAGCCg                                                                            <  2:83660/148‑1 (MQ=255)
                                                                       ggACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGTGAGCATGCCCAGCACGGTCAACTGCAGCAGTTGCACGGGGAACACATACAGCGGTAAACGCCAGCCGAAgg                                                                        >  8:274261/1‑149 (MQ=255)
                                                                         aCAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGTGAGCATGCCCAGCACGGTCAACTGCAGCAGTTGCACGGGGAACACATACAGCGGTAAACGCCAGCCGAAGGcc                                                                      <  5:3099/149‑1 (MQ=255)
                                                                             cGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGTGAGCATGCCCAGCACGGTCAACTGCAGCAGTTGCACGGGGAACACATACAGCGGTAAACGCCAGCCGAAGGCCTGCa                                                                  >  8:85083/1‑149 (MQ=255)
                                                                                    aaCAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTg                                                                                                                                                 <  8:19786/63‑1 (MQ=255)
                                                                                                  cgggcggCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGTGAGCATGCCCAGCACGGTCAACTGCAGCAGTTGCACGGGGAACACATACAGCGGTAAACGCCAGCCGAAGGCCTGCACGTTCACCACTGCGaccaggc                                             <  3:154142/149‑1 (MQ=255)
                                                                                                      cggcgCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGTGAGCATGCCCAGCACGGTCAACTGCAGCAGTTGCACGGGGAACACATACAGCGGTAAACGCCAGCCGAAGGCCTGCACGTTCACCACTGCGaccaggcac                                           >  6:198589/1‑147 (MQ=255)
                                                                                                       ggcgCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGTGAGCATGCCCAGCACGGTCAACTGCAGCAGTTGCACGGGGAACACATACAGCGGTAAACGCCAGCCGAAGGCCTGCACGTTCACCACTGCGaccaggcaccag                                        <  8:271920/149‑1 (MQ=255)
                                                                                                            gcCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGTGAGCATGCCCAGCACGGTCAACTGCAGCAGTTGCACGGGGAACACATACAGCGGTAAACGCCAGCCGAAGGCCTGCACGTTCACCACTGCGACCAGGCACCAGGCCAg                                   <  6:291118/149‑1 (MQ=255)
                                                                                                             ccAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGTGAGCATGCCCAGCACGGTCAACTGCAGCAGTTGCACGGGGAACACATACAGCGGTAAACGCCAGCCGAAGGCCTGCACGTTCACCACTGCGACCAGGCACCAGGCCAGc                                  <  8:91120/149‑1 (MQ=255)
                                                                                                               aGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGTGAGCATGCCCAGCACGGTCAACTGCAGCAGTTGCACGGGGAACACATACAGCGGTAAACGCCAGCCGAAGGCCTGCACGTTCACCACTGCGACCAGGCACCAGGCCAGCaa                                <  1:83412/149‑1 (MQ=255)
                                                                                                                    ccACAGCGGCCAGGCGCAGGCCAGCAGGCTGGTGAGCATGCCCAGCACGGTCAACTGCAGCAGTTGCACGGGGAACACATACAGCGGTAAACGCCAGCCGAAGGCCTGCACGTTCACCACTGCGACCAGGCACCAGGCCAGCAACAGc                            <  7:274419/148‑1 (MQ=255)
                                                                                                                        aGCGTCCAGGCGCAGGCCAGCAGGCTGGTGAGCATGCCCAGCACGGTCAACTGCAGCAGTTGCACGGGGAACACATACAGCGGTAAACGCCAGCCGAAGGCCTGCACGTTCACCACTGCGACCAGGCACCAGGCCAGCAACAGCCCAAg                       <  6:239851/149‑1 (MQ=255)
                                                                                                                          cGGCCAGGCGCAGGCCAGCAGGCTGGTGAGCATGCCCAGCACGGTCAACTGCAGCAGTTGCACGGGGAACACATACAGCGGTAAACGCCAGCCGAAGGCCTGCACGTTCACCACTGCGACCAGGCACCAGGCCAGCAACAGCCCAAGCg                     >  7:57999/1‑149 (MQ=255)
                                                                                                                          cGGCCAGGCGCAGGCCAGCAGGCTGGTGAGCATGCCCAGCACGGTCAACTGCAGCAGTTGCACGGGGAACACATACAGCGGTAAACGCCAGCCGAAGGCCTGCACGTTCACCACTGCGACCAGGCACCAGGCCACCAACAGCCCAAGCg                     <  7:274261/149‑1 (MQ=255)
                                                                                                                                         cagcagGCTGGTGAGCATGCCCAGCACGGTCAACTGCAGCAGTTGCACGGGGAACACATACAGCGGTAAACGCCAGCCGAAGGCCTGCACGTTCACCACTGCGACCAGGCACCAGGCCAGCAACAGCCCAAGCGGAATCGCCAGCAGTa      >  4:185835/1‑149 (MQ=255)
                                                                                                                                             agGCTGGTGAGCATGCCCAGCACGGTCAACTGCAGCAGTTGCACGGGGAACACATACAGCGGTAAACGCCAGCCGAAGGCCTGCACGTTCACCACTGCGACCAGGCACCAGGCCAGCAACAGCCCAAGCGGAATCGCCAGCAGTACGGt  >  1:245297/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
TTTGCGGTGCAGGCTGATCACAGCCGCACAACAACAACCCGCTCAGCAACAGCCACGCACTAAGCTTCATCGGACAACGGTCTCAACAGTTCGCGCGGCGGGCGGCGCGCCAGCTGCCACAGCGGCCAGGCGCAGGCCAGCAGGCTGGTGAGCATGCCCAGCACGGTCAACTGCAGCAGTTGCACGGGGAACACATACAGCGGTAAACGCCAGCCGAAGGCCTGCACGTTCACCACTGCGACCAGGCACCAGGCCAGCAACAGCCCAAGCGGAATCGCCAGCAGTACGGT  >  NC_002947/4260243‑4260532

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