Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,542,684 C→T 73.5% S48S (AGC→AGT)  PP_2235 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,542,6840CT73.5% 74.3 / 28.8 49S48S (AGC→AGT) PP_2235hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (7/6);  new base T (22/14);  total (29/20)
Fisher's exact test for biased strand distribution p-value = 7.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

CACCTCATGCCTTGTGGTTCACCCAGATGATCGAATACGAAGTACCCAGTATGTACCAGCAAGCATTGGCTCAGGCCCTGGTGGAGCGTAGTGAAGCATTGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGCATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTACCTGCAATGGCAGTCACGCCAGGCCTGGGCTGCGGCGGCCGTGCACTTCGTCGAGGAATCGTTTCTTGACCTGCTGGGCCAGCATCAGGCACGCGGTGTCAACT  >  NC_002947/2542551‑2542826
                                                                                                                                     |                                                                                                                                              
cacCTCATGCCTTGTGGTTCACCCAGATGATCGAATACGAAGTACCCAGTATGTACCAGCAAGCATTGGCTCAGGCCCTGGTGGAGCATAGTGAAGCATTGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGCATCCAGGCCAGcgac                                                                                                                                 >  7:163839/1‑149 (MQ=255)
  ccTCATGCCTTGTGGTTCACCCAGATGATCGAATACGAAGTACCCAGTATGTACCAGCAAGCATTGGCTCAGGCCCTGGTGGAGCGTAGTGAAGCATTGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGTATCCAGGCCAGcgacga                                                                                                                               <  8:117526/149‑1 (MQ=255)
     cATGCCTTGTGGTTCACCCAGATGATCGAATACGAAGTACCCAGTATGTACCAGCAAGCATTGGCTCAGGCCCTGGTGGAGCGTAGTGAAGCATTGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGTATCCAGGCCAGCGACGATgg                                                                                                                            >  6:158518/1‑149 (MQ=255)
        gCCTTGTGGTTCACCCAGATGATCGAATACGAAGTACCCAGTATGTACCAGCAAGCATTGGCTCAGGCCCTGGTGGAGCGTAGTGAAGCATTGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGTATCCAGGCCAGCGACGATGGCTg                                                                                                                         >  1:186294/1‑149 (MQ=255)
          cTTGTGGTTCACCCAGATGATCGAATACGAAGTACCCAGTATGTACCAGCAAGCATTGGCTCAGGCCCTGGTGGAGCGTAGTGAAGCATTGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGTATCCAGGCCAGCGACGATGGCTGcc                                                                                                                       >  5:39971/1‑149 (MQ=255)
            tgtgGTTCACCCAGATGATCGAATACGAAGTACCCAGTATGTACCAGCAAGCATTGGCTCAGGCCCTGGTGGAGCGTAGTGAAGCATTGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGTATCCAGGCCAGCGACGATGGCTGCCgg                                                                                                                     >  4:357329/1‑149 (MQ=255)
                    aCCCAGATGATCGAATACGAAGTACCCAGTATGTACCAGCAAGCATTGGCTCAGGCCCTGGTGGAGCATAGTGAAGCATTGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGCATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCa                                                                                                             <  8:163839/149‑1 (MQ=255)
                     cccAGATGATCGAATACGAAGTACCCAGTATGTACCAGCAAGCATTGGCTCAGGCCCTGGTGGAGCGTAGTGAAGCATTGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGTATCCAGGCCAGCGACGATgg                                                                                                                            <  4:237111/133‑1 (MQ=255)
                     cccAGATGATCGAATACGAAGTACCCAGTATGTACCAGCAAGCATTGGCTCAGGCCCTGGTGGAGCGTAGTGAAGCATTGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGTATCCAGGCCAGCGACGATgg                                                                                                                            >  3:237111/1‑133 (MQ=255)
                      ccAGATGATCGAATACGAAGTACCCAGTATGTACCAGCAAGCATTGGCTCAGGCCCTGGTGGAGCGTAGTGAAGCATTGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGTATCCAGGCCAGCGACGATGGCTGCCGGgtgt                                                                                                                 <  4:32720/143‑1 (MQ=255)
                      ccAGATGATCGAATACGAAGTACCCAGTATGTACCAGCAAGCATTGGCTCAGGCCCTGGTGGAGCGTAGTGAAGCATTGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGTATCCAGGCCAGCGACGATGGCTGCCGGgtgt                                                                                                                 >  3:32720/1‑143 (MQ=255)
                        aGATGATCGAATACGAAGTACCCAGTATGTACCATCAAGCATTGGCTCAGGCCCTGGTGGAGCGTAGTGAAGCATTGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGTATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTAc                                                                                                         <  6:39971/149‑1 (MQ=255)
                           tgatCGAATACGAAGTACCCAGTATGTACCAGCAAGCATTGGCTCAGGCCCTGGTGGAGCGTAGTGAAGCATTGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGTATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTACCTg                                                                                                      <  5:158518/149‑1 (MQ=255)
                                gAATACGAAGTACCCAGTATGTACCAGCAAGCATTGGCTCAGGCCCTGGTGGAGCGTAGTGAAGCATTGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGCATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTACCTGCAATg                                                                                                 >  5:144543/1‑149 (MQ=255)
                                   tACGAAGTACCCAGTATGTACCAGCAAGCATTGGCTCAGGCCCTGGTGGAGCGTAGTGAAGCATTGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGTATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTACCTGCAATGGCa                                                                                              <  2:25889/149‑1 (MQ=255)
                                     cGAAGTACCCAGTATGTACCAGCAAGCATTGGCTCAGGCCCTGGTGGAGCGTAGTGAAGCATTGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGCATCCAGGCCAGCGACGATGGCTGCCGGGTGtt                                                                                                                <  7:73533/129‑1 (MQ=255)
                                     cGAAGTACCCAGTATGTACCAGCAAGCATTGGCTCAGGCCCTGGTGGAGCGTAGTGAAGCATTGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGCATCCAGGCCAGCGACGATGGCTGCCGGGTGtt                                                                                                                >  8:73533/1‑129 (MQ=255)
                                             ccAGTATGTACCAGCAAGCATTGGCTCAGGCCCTGGTGGAGCGTAGTGAAGCATTGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGTATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTACCTGCAATGGCAGTCACGCCAg                                                                                    >  6:348297/1‑149 (MQ=255)
                                             ccAGTATGTACCAGCAAGCATTGGCTCAGGCCCTGGTGGAGCGTAGTGAAGCATTGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGTATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTACCTGCAATGGCAGTCACGCCAg                                                                                    >  2:215936/1‑149 (MQ=255)
                                             ccAGTATGTACCAGCAAGCATTGGCTCAGGCCCTGGTGGAGCGTAGTGAAGCATTGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGCATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTAc                                                                                                         <  3:119034/128‑1 (MQ=255)
                                             ccAGTATGTACCAGCAAGCATTGGCTCAGGCCCTGGTGGAGCGTAGTGAAGCATTGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGCATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTAc                                                                                                         >  4:119034/1‑128 (MQ=255)
                                                    gTACCAGCAAGCATTGGCTCAGGCCCTGGTGGAGCGTAGTGAAGCATTGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGTATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTACCTGCAATGGCAGTCACGCCAGGCCTggg                                                                             >  2:314080/1‑149 (MQ=255)
                                                        cAGCAAGCATTGGCTCAGGCCCTGGTGGAGCGTAGTGAAGCATTGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGCATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTACCTGCAATGGCAGTCACGCCAGGCCTGGGCTgc                                                                         >  1:277534/1‑149 (MQ=255)
                                                                aTTGGCTCAGGCCCTGGTGGAGCGTAGTGAAGCATTGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGCATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTACCTGCAATGGCAGTCACGCCAGGCCTGGGCTGCGGCGGCCg                                                                 >  7:96529/1‑149 (MQ=255)
                                                                    gCTCAGGCCCTGGTGGAGCGTAGTGAAGCATTGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGTATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTACCTGCAATGGCAGTCACGCCAGGCCTGGGCTGCGGCGGCCGTGCa                                                             <  5:348297/149‑1 (MQ=255)
                                                                      tCAGGCCCTGGTGGAGCGTAGTGAAGCATTGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGCATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTACCTGCAATGGCAGTCACGCCAGGCCTGGGCTGCGGCGGCCGTGCACt                                                           >  1:336050/1‑149 (MQ=255)
                                                                          gCCCTGGTGGAGCGTAGTGAAGCATTGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGTATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTACCTGCAATGGCAGTCACGCCAGGCCTGGGCTGCGGCGGCCGTGCACTtcg                                                        >  2:139600/1‑148 (MQ=255)
                                                                            ccTGGTGGAGCGTAGTGAAGCATTGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGTATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTACCTGCAATGGCAGTCACGCCAGGCCTGGGCTGCGGCGGCCGTGCACTtcgtc                                                      >  5:116759/1‑148 (MQ=255)
                                                                             cTGGTGGAGCGTAGTGAAGCATTGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGTATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTACCTGCAATGGCAGTCACGCCAGGCCTGGGCTGCGGCGGCCGTGCACTTCGTCGa                                                    <  4:156371/149‑1 (MQ=255)
                                                                             cTGGTGGAGCGTAGTGAAGCATTGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGCATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTACCTGCAATGGCAGTCACGCCAGGCCTGGGCTGCGGCGGCCGTGCACTTCGTCGa                                                    <  2:277534/149‑1 (MQ=255)
                                                                                gtggAGCGTAGTGAAGCATTGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGTATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTACCTGCAATGGCAGTCACGCCAGGCCTGGGCTGCGGCGGCCGTGCACTTCGTCGAg                                                   <  3:357329/147‑1 (MQ=255)
                                                                                     gCGTAGTGAAGCATTGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGTATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTACCTGCAATGGCAGTCACGCCAGGCCTGGGCTGCGGCGGCCGTGCACTTCGTCGAGGAATCGt                                            >  3:102977/1‑149 (MQ=255)
                                                                                       gTAGTGAAGCATTGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGTATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTACCTGCAATGGCAGTCACGCCAGGCCTGGGCTGCGGCGGCCGTGCACTTCGTCGAGGAATCGttt                                          <  1:215936/149‑1 (MQ=255)
                                                                                          gTGAAGCATTGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGTATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTACCTGCAATGGCAGTCACGCCAGGCCTGGGCTGCGGCGGCCGTGCACTTCGTCGAGGAATCGTTTCtt                                       <  1:314080/149‑1 (MQ=255)
                                                                                             aaGCATTGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGTATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTACCTGCAATGGCAGTCACGCCAGGCCTGGGCTGCGGCGGCCGTGCACTTCGTCGAGGAATCGTTTCTTGAc                                    >  4:542/1‑149 (MQ=255)
                                                                                                  ttGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGTATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTACCTGCAATGGCAGTCACGCCAGGCCTGGGCTGCGGCGGCCGTGCACTTCGTCGAGGAATCGTTTCTTGACctgct                               >  3:108183/1‑149 (MQ=255)
                                                                                                       tgctCGTTGCGAAGGGCTGCAGGGTGTGAGTATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTACCTGCAATGGCAGTCACGCCAGGCCTGGGCTGCGGCGGCCGTGCACTTCGTCGAGGAATCGTTTCTTGACCTGCTGGGc                           >  4:36802/1‑148 (MQ=255)
                                                                                                            gTTGCGAAGGGCTGCAGGGTGTGAGTATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTACCTGCAATGGCAGTCACGCCAGGCCTGGGCTGCGGCGGCCGTGCACTTCGTCGAGGAATCGTTTCTTGACCTGCTGGGCCAGCAt                     <  2:186294/149‑1 (MQ=255)
                                                                                                                 gAAGGGCTGCAGGGTGTGAGTATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTACCTGCAATGGCAGTCACGCCAGGCCTGGGCTGCGGCGGCCGTGCACTTCGTCGAGGAATCGTTTCTTGACCTGCTGGGCCAGCATCAgg                 >  7:71371/1‑148 (MQ=255)
                                                                                                                   aGGGCTGCAGGGTGTGAGTATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTACCTGCAATGGCAGTCACGCCAGGCCTGGGCTGCGGCGGCCGTGCACTTCGTCGAGGAATCGTTTCTTGACCTGCTGGGCCAGCATCAGGCAc              >  2:207375/1‑149 (MQ=255)
                                                                                                                   aGGGCTGCAGGGTGTGAGTATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTACCTGCAATGGCAGTCACGCCAGGCCTGGGCTGCGGCGGCCGTGCACTTCGTCGAGGAATCGTTTCTTGACCTGCTGGGCCAGCATCAGGCAc              >  8:297411/1‑149 (MQ=255)
                                                                                                                      gCTGCAGGGTGTGAGTATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTACCTGCAATGGCAGTCACGCCAGGCCTGGGCTGCGGCGGCCGTGCACTTCGTCGAGGAATCGTTTCTTGACCTGCTGGGCCAGCATCAGGCAcgcg           >  8:100230/1‑149 (MQ=255)
                                                                                                                        tGCAGGGTGTGAGCATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTACCTGCAATGGCAGTCACGCCAGGCCTGGGCTGCGGCGGCCGTGCACTTCGTCGAGGAATCGTTTCTTGACCTGCTGGGCCAGCATCAGGCACGCGgt         <  8:96529/149‑1 (MQ=255)
                                                                                                                          cAGGGTGTGAGTATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTACCTGCAATGGCAGTCACGCCAGGCCTGGGCTGCGGCGGCCGTGCACTTCGTCGa                                                    <  7:105730/104‑1 (MQ=255)
                                                                                                                          cAGGGTGTGAGTATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTACCTGCAATGGCAGTCACGCCAGGCCTGGGCTGCGGCGGCCGTGCACTTCGTCGa                                                    >  8:105730/1‑104 (MQ=255)
                                                                                                                              gtgtGAGCATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTACCTGCAATGGCAGTCACGCCAGGCCTGGGCTGCGGCGGCCGTGCACTTCGTCGAGGAATCGTTTCTTGACCTGCTGGGCCAGCATCAGGCACGCAGTGTCAAc   <  2:336050/149‑1 (MQ=255)
                                                                                                                               tgtgAGTATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTACCTGCAATGGCAGTCACGCCAGGCCTGGGCTGCGGCGGCCGTGCACTTCGTCGAGGAATCGTTTCTTGACCTGCTGGGCCAGCATCAGGCACGCGGTGTCAACt  >  7:37984/1‑149 (MQ=255)
                                                                                                                                  gAGTATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTACCTGCAATGGCAGTCACGCCAGGCCTGGGCTGCGGCGGCCGTGCACTTCGTCGAGGAATCGTTTCTTGACCTGCTGGGCCAGCATCAGGCACGCGgtgt       >  7:227516/1‑141 (MQ=255)
                                                                                                                                  gAGTATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTACCTGCAATGGCAGTCACGCCAGGCCTGGGCTGCGGCGGCCGTGCACTTCGTCGAGGAATCGTTTCTTGACCTGCTGGGCCAGCATCAGGCACGCGgtgt       <  8:227516/141‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                              
CACCTCATGCCTTGTGGTTCACCCAGATGATCGAATACGAAGTACCCAGTATGTACCAGCAAGCATTGGCTCAGGCCCTGGTGGAGCGTAGTGAAGCATTGGCTGCTCGTTGCGAAGGGCTGCAGGGTGTGAGCATCCAGGCCAGCGACGATGGCTGCCGGGTGTTGCAGTACCTGCAATGGCAGTCACGCCAGGCCTGGGCTGCGGCGGCCGTGCACTTCGTCGAGGAATCGTTTCTTGACCTGCTGGGCCAGCATCAGGCACGCGGTGTCAACT  >  NC_002947/2542551‑2542826

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