Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,422,811 T→A 8.7% A323A (GCT→GCA) ‡ PP_3038 → TraC domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,422,8110TA 8.7% 129.6 / 3.7 46A323A (GCT→GCA) ‡PP_3038TraC domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (23/19);  new base A (2/2);  total (25/21)
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.36e-01

CGGTGGTTGTGGGTGGCCAGGTAGTCGGCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGTTACAACGCCTTCAGCTCCAAGATCATGC  >  NC_002947/3422669‑3422953
                                                                                                                                              |                                                                                                                                                
cGGTGGTTGTGGGTGGCCAGGTAGTCGGCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAATGATCAAGCCACCGGCTACCGGt                                                                                                                                            >  4:185540/1‑149 (MQ=255)
  gtggtTGTGGGTGGCCAGGTAGTCGGCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCCAGCCACCGGTAGCCGGTgg                                                                                                                                          <  7:27682/149‑1 (MQ=255)
       tgtgGGTGGCCAGGTAGTCGGCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTgg                                                                                                                                      >  4:309022/1‑148 (MQ=255)
                       gTCGGCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAg                                                                                                                     <  3:309022/149‑1 (MQ=255)
                        tCGGCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAg                                                                                                                     >  3:54768/1‑148 (MQ=255)
                           gCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACt                                                                                                                 >  3:291270/1‑149 (MQ=255)
                              cGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACt                                                                                                                 >  5:276163/1‑146 (MQ=255)
                               gATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAgg                                                                                                                                  >  4:215997/1‑128 (MQ=255)
                               gATCTTCTCCGCTGAGCGTGAGGTTAAATGGACCGACTTGAACGTCCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGGCGGCATGTCACGGCAGTGATCAAGCCTCCGGCTACCGGTGGCTGGAAgg                                                                                                                                  <  3:215997/128‑1 (MQ=255)
                                     ctcGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAgtggtg                                                                                                       >  7:73989/1‑149 (MQ=255)
                                     ctcGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAgtggt                                                                                                        <  4:54768/148‑1 (MQ=255)
                                     ctcGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAATGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAgtg                                                                                                          <  3:185540/146‑1 (MQ=255)
                                        ggCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGccc                                                                                                    >  5:223543/1‑149 (MQ=255)
                                        ggCGGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCTCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGccc                                                                                                    <  5:186505/149‑1 (MQ=255)
                                           tGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGa                                                                                                 >  5:189326/1‑149 (MQ=255)
                                             aGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATc                                                                                               <  4:27985/149‑1 (MQ=255)
                                                     gTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGATGGGCTTGAGGGCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGGCAGGACTGAGAGTGGTGCCCTGATCGCCCatat                                                                                       <  2:292121/149‑1 (MQ=255)
                                                      ttAAATGGACCGACTTCAACGCCCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATg                                                                                      >  2:39697/1‑149 (MQ=255)
                                                        aaaTGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGc                                                                                     <  8:73989/148‑1 (MQ=255)
                                                         aaTGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGCCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATg                                                                                      >  3:106055/1‑146 (MQ=255)
                                                          aTGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAg                                                                                   >  2:84507/1‑148 (MQ=255)
                                                            ggACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGCCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGaa                                                                                 <  4:106055/148‑1 (MQ=255)
                                                              aCCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTACCCTGATCGCCCATATGCAGAACGt                                                                              >  5:191989/1‑149 (MQ=255)
                                                                  aCTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGaa                                                                          <  6:189326/149‑1 (MQ=255)
                                                                        aCGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATc                                                                    <  4:291270/149‑1 (MQ=255)
                                                                           aCCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTg                                                                 >  1:17648/1‑149 (MQ=255)
                                                                               gCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGcc                                                              >  1:285675/1‑148 (MQ=255)
                                                                               gCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGCGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGTCTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGcc                                                              >  2:91146/1‑148 (MQ=255)
                                                                                 aTTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGTAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAat                                                          <  2:17648/149‑1 (MQ=255)
                                                                                    gcgcTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAatgat                                                      >  6:62504/1‑149 (MQ=255)
                                                                                         gAAGGCGTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATCCAGAACGTCGAACTGATCCTGGCCTATGATGAgcg                                                 >  6:150351/1‑149 (MQ=255)
                                                                                                   aGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCgg                                        <  6:223543/148‑1 (MQ=255)
                                                                                                     gCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGt                                     <  1:84507/149‑1 (MQ=255)
                                                                                                         ttCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATc                                 >  5:61360/1‑149 (MQ=255)
                                                                                                            gCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGt                              <  5:62504/149‑1 (MQ=255)
                                                                                                                gCAACTCACGGGAGTGATGAAGCCACCGGCTACCGGTGGGTGGAATGACATGCTGGCCAGGCCTGCGAGTGGTGCCCTGCTCGCCCATATTCAGAACGTCGTACTGTTTCTGGCCAAATGATGAGCGCTGGAGCGGGGTGATCGGTTAc                           <  1:39697/149‑1 (MQ=255)
                                                                                                                          gCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGAAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGTTACAACGCCTTCAg                >  1:19146/1‑149 (MQ=255)
                                                                                                                               cttccaGCCACCGGTAGCCGGTGGCTTGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGTTACAACGCCTTCAGCTcc            >  3:114818/6‑148 (MQ=255)
                                                                                                                                aTCAAGCCACCGGCTACCGGTGGCTGGATGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAatg                                                        >  1:184172/1‑103 (MQ=255)
                                                                                                                                aTCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAatg                                                        <  2:184172/103‑1 (MQ=255)
                                                                                                                                  ccaGCCACCGGTAGCCGGTGGCTTGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGTTACAACGCCTTCAGCTCCaa          <  4:114818/145‑1 (MQ=255)
                                                                                                                                   aaGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATTCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGTTACAACGCCTTCAGCTCCAAGAt       <  6:276163/149‑1 (MQ=255)
                                                                                                                                    aGCCACCGGTAGCCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGTTACAACGCCTTCAGCTCCAAGAt       >  4:137005/1‑148 (MQ=255)
                                                                                                                                      ccACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGTTACAACGCCTTCAGCTCCAAGATCAt    >  1:200943/1‑149 (MQ=255)
                                                                                                                                       cACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCaa                                                           >  8:288672/1‑94 (MQ=255)
                                                                                                                                         ccGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGTCTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATCGATGAGCGCTGGAGCGGGGTGATCGGTTACAACGCCTTCAGCTCCAAGATCATgc  <  1:91146/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                
CGGTGGTTGTGGGTGGCCAGGTAGTCGGCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGTTACAACGCCTTCAGCTCCAAGATCATGC  >  NC_002947/3422669‑3422953

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: