Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,443,932 T→C 8.9% R897R (CGT→CGC)  PP_3062 → tail length determination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,443,9320TC 8.9% 120.5 / 3.3 46R897R (CGT→CGC) PP_3062tail length determination protein
Reads supporting (aligned to +/- strand):  ref base T (20/21);  new base C (2/2);  total (23/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 = 6.62e-01

TGGTGCGGCCACCGGCCAGCCCTGAGTCGCCGGCAGATCCTGCAGCACCGGAGGCACAACCCGAGCCGTCGGTGCCCAAGCTGGGCGATACGGTGCGGGAGGTGGCCAGCCCGGCGCCGGTGGAGCCAGCGGTGTCGTATGACCCGCGTGACCCAGCGTCGAAAGATCCGTTCCTGCTGCCGGCGCTGACCGCCAACAAGGTGCGGTTCCCTGGTGCTGGACTGGTACGGCCGCAGGCCGAGCCCGAGCCGCCGACACCCAAACTGGGGGCTGCGGTGCGGGAAGCAGCACCGGC  >  NC_002947/3443784‑3444078
                                                                                                                                                    |                                                                                                                                                  
tGGTGCGGCCACCGGCCAGCCCTGAGTCGCCGGCAGATCCTGCAGCACCGGAGGCACAACCCGAGCCGTCGGTGCCCAAGCTGGGCGATACGGTGCGGGAGGTGGCCAGCCCGGCGCCGGTGGAGCCAGCGGTGTCGTATGACCCGCGt                                                                                                                                                    >  5:202010/1‑149 (MQ=255)
        ccACCGGCCAGCCCTGAGTCGCCGGCAGATCCTGCAGCACCGGAGGCACAACCCGAGCCGTCGGTGCCCAAGCTGGGCGATACGGTGCGGGAGGTGGCCAGCCCGGCGCCGGTGGAGCCAGCGGTGTCGTATGACCCGCGTGACCCAGc                                                                                                                                            >  6:79802/1‑149 (MQ=255)
         cACCGGCCAGCCCTGAGTCGCCGGCAGATCCTGCAGCACCGGAGGCACAACCCGAGCCGTCGGTGCCCAAGCTGGGCGATACGGTGCGGGAGGTGGCCAGCCCGGCGCCGGTGGAGCCAGCGGTGTCGTATGACCCGCGTGACCCAGCg                                                                                                                                           >  7:255845/1‑149 (MQ=255)
               ccAGCCCTGAGTCGCCGGCAGATCCTGCAGCACCGGAGGCACAACCCGAGCCGTCGGTGCCCAAGCTGGGCGATACGGTGCGGGAGGTGGCCAGCCCGGCGCCGGTGGAGCCAGCGGTGTCGTATGACCCGCGTGACCCAGCGTCGaaa                                                                                                                                     <  2:106283/149‑1 (MQ=255)
                 aGCCCTGAGTCGCCGGCAGATCCTGCAGCACCGGAGGCACAACCCGAGCCGTCGGTGCCCAAGCTGGGCGATACGGTGCGGGAGGTGGCCAGCCCGGCGCCGGTGGAGCCAGCGGTGTCGTATGACCCGCGTGACCCAGCGTCGa                                                                                                                                       >  8:253259/1‑145 (MQ=255)
                  gCCCTGAGTCGCCGGCAGATCCTGCAGCACCGGAGGCACAACCCGAGCCGTCGGTGCCCAAGCTGGGCGATACGGTGCGGGAGGTGGCCAGCCCGGCGCCGGTGGAGCCAGCGGTGTCGTATGACCCGCGTGACCCAGCGTCGAAAGAt                                                                                                                                  <  1:241770/149‑1 (MQ=255)
                    ccTGAGTCGCCGGCAGATCCTGCAGCACCGGAGGCACAACCCGAGCCGTCGGTGCCCAAGCTGGGCGATACGGTGCGGGAGGTGGCCAGCCCGGCGCCGGTGGAGCCAGCGGTGTCGTATGACCCGCGTGACCCAGCGTCGAAAGATcc                                                                                                                                <  4:114409/149‑1 (MQ=255)
                             ccGGCAGATCCTGCAGCACCGGAGGCACAACCCGAGCCGTCGGTGCCCAAGCTGGGCGATACGGTGCGGGAGGTGGCCAGCCCGGCGCCGGTGGAGCCAGCGGTGTCGTATGACCCGCGTGACCCAGCGTCGAAAGATCCGTTCctgct                                                                                                                       <  1:74455/149‑1 (MQ=255)
                              cGGCAGATCCTGCAGCACCGGAGGCACAACCCGAGCCGTCGGTGCCCAAGCTGGGCGATACGGTGCGGGAGGTGGCCAGCCCGGCGCCGGTGGAGCCAGCGGTGTCGTATGACCCGCGTGACCCAGCGTCGAAAGATCCGTTCctgctg                                                                                                                      >  4:296395/1‑149 (MQ=255)
                                  aGATCCTGCAGCACCGGAGGCACAACCCGAGCCGTCGGTGCCCAAGCTGGGCGATACGGTGCGGGAGGTGGCCAGCCCGGCGCCGGTGGAGCCAGCGGTGTCGTATGACCCGCGTGACCCAGCGTCGAAAGATCCGTTCCTGATGCCgg                                                                                                                  >  7:273360/1‑149 (MQ=255)
                                    aTCCTGCAGCACCGGAGGCACAACCCGAGCCGTCGGTGCCCAAGCTGGGCGATACGGTGCGGGAGGTGGCCAGCCCGGCGCCGGTGGAGCCAGCGGTGTCGTATGACCCGCGTGACCCAGCGTCGAAAGATCCGTTCCTGCTGCCGgcg                                                                                                                >  5:198973/1‑149 (MQ=255)
                                          cagcaCCGGAGGCACAACCCGAGCCGTCGGTGCCCAAGCTGGGCGATACGGTGCGGGAGGTGGCCAGCCCGGCGCCGGTGGAGCCAGCGGTGTCGTATGACCCGCGTGACCCAGCGTCGAAAGATCCGTTCCTGCTGCCGGCGCTGAcc                                                                                                          <  7:253259/149‑1 (MQ=255)
                                                       acaACCCGAGCCGTCGGTGCCCAAGCTGGGCGATACGGTGCGGGAGGTGGCCAGCCCGGCGCCGGTGGAGCCAGCGGTGTCGTATGACCCGCGTGACCCAGCGTCGAAAGATCCGTTCCTGCTGCCGGCGCTGACCGCCAACAAGGTGc                                                                                             >  4:113810/1‑149 (MQ=255)
                                                       acaACCCGAGCCGTCGGTGCCCAAGCTGGGCGATACGGTGCGGGAGGTGGCCAGCCCGGCGCCGGTGGAGCCAGCGGTGTCGTATGACCCGCGTGACCCAGCGTCGAAAGATCCGTTCCTGCTGCCGGCGCTGACCGCCAACAAGGTGc                                                                                             >  8:178722/1‑149 (MQ=255)
                                                       acaACCCGAGCCGTCGGTGCCCAAGCTGGGCGATACGGTGCGGGAGGTGGCCAGCCCGGCGCCGGTGGAGCCAGCGGTGTCGTATGACCCGCGTGACCCAGCGTCGAAAGATCCGTTCCTGCTGCCGGCGCTGACCGCCAACAAGGTGc                                                                                             >  4:70753/1‑149 (MQ=255)
                                                       acaACCCGAGCCGTCGGTGCCCAAGCTGGGCGATACGGTGCGGGAGGTGGCCAGCCCGGCGCCGGTGGAGCCAGCGGTGTCGTATGACCCGCGTGACCCAGCGTCGAAAGATCCGTTCCTGCTGCCGGCGCTGACCGCCAACAAGGTGc                                                                                             >  3:31590/1‑149 (MQ=255)
                                                         aaCCCGAGCCGTCGGTGCCCAAGCTGGGCGATACGGTGCTGGAGGTGGCCAGCCCGGCGCCGGTGGAGCCAGCGGTGTCGTATGACCCGCGTGACCCAGCGTCGAAAGATCCGTTCCTGCTGCCGGCGCTGACCGCCAACAAGGTGCgg                                                                                           <  5:181708/149‑1 (MQ=255)
                                                         aaCCCGAGCCGTCGGTGCCCAAGCTGGGCGATACGGTGCGGGAGGTGGCCAGCCCGGCGCCGGTGGAGCCAGCGGTGTCGTATGACCCGCGTGACCCAGCGTCGAAAGATCCGTTCCTGATGCCCGCCATGACAGCGCACAAgagcctc                                                                                           >  8:77450/1‑143 (MQ=37)
                                                         aaCCCGAGCCGTCGGTGCCCAAGCTGGGCGATACGGTGCGGGAGGTGGCCAGCCCGGCGCCGGTGGAGCCAGCGGTGTCGTATGAACCGCGTGACCCAGCGTCGAAAGATCCGTTCCTGATGCCCGCNATGACAGCAAACAAGNTGCgg                                                                                           >  8:82021/1‑149 (MQ=37)
                                                           cccGAGCCGTCGGTGCCCAAGCTGGGCGATACGGTGCGGGAGGTGGCCAGCCCGGCGCCGGTGGAGCCAGCGGTGTCGTATGACCCGCGTGACCCAGCGTCGAAAGATCCGTTCCTGCTGCCGGCGCTGACCGCCAACAAGGTGCGGtt                                                                                         <  5:79802/149‑1 (MQ=255)
                                                               aGCCGTCGGTGCCCAAGCTGGGCGATACGGTGCGGGAGGTGGCCAGCCCGGCGCCGGTGGAGCCAGCGGTGTCGTATGACCCGCGTGACCCAGCGTCGAAAGATCCGTTCCTGCTGCCGGCGCTGACCGCCAACAAGGTGCGGTTCCCt                                                                                     <  6:202010/149‑1 (MQ=255)
                                                                      ggTGCCCAAGCTGGGCGATACGGTGCGGGAGGTGGCCAGCCCGGCGCCGGTGGAGCCAGCGGTGTCGTATGACCCGCGTGACCCAGCGTCGAAAGATCCGTTCCTGCTGCCGGCGCTGACCGCCAACAAGGTGCGGTTCCCTGGTGCTg                                                                              <  3:70753/149‑1 (MQ=255)
                                                                           ccAAGCTGGGCGATACGGTGCGGGAGGTGGCCAGCCCGGCGCCGGTGGAGCCAGCGGTGTCGTATGACCCGCGTGACCCAGCGTCGAAAGATCCGTTCCTGCTGCCGGCGCTGACCGCCAACAAGGTGCGGTTCCCTGGTGCTGGACTg                                                                         >  8:74374/1‑149 (MQ=255)
                                                                            cAAGCTGGGCGATACGGTGCGGGAGGTGGCCAGCCCGGCGCCGGTGGAGCCAGCGGTGTCGTATGACCCGCGTGACCCAGCGTCGAAAGATCCGTTCCTGCTGCCGGCGCTGACCGCCAACAAGGTGCGGTTCCCTGGTGCTGGACTg                                                                         <  7:156539/148‑1 (MQ=255)
                                                                            cAAGCTGGGCGATACGGTGCGGGAGGTGGCCAGCCCGGCGCCGGTGGAGCCAGCGCTGTCGTATGACCCGCGTGACCCAGCGTCGAAAGATCCGTTCCTGCTGCCGGCGCTGACCGCCAACAAGGTGCGGTTCCCTGCTGCTGGACTg                                                                         >  8:156539/1‑148 (MQ=255)
                                                                               gCTGGGCGATACGGTGCGGGAGGTGGCCAGCCCGGCGCCGGTGGAGCCAGCGGTGTCGTATGACCCGCGTGACCCAGCGTCGAAAGATCCGTTCCTGCTGCCGGCGCTGACCGCCAACAAGGTGCGGTTCCCTGGTGCTGGACTGGTAc                                                                     >  7:89137/1‑149 (MQ=255)
                                                                                   ggCGATACGGTGCGGGAGGTGGCCAGCCCGGCGCCGGTGGAGCCAGCGGTGTCGTATGACCCGCGTGACCCAGCGTCGAAAGATCCGTTCCTGCTGCCGGCGCTGACCGCCAACAAGGTGCGGTTCCCTGGTGCTGGACTGGTACgg                                                                   <  7:74374/147‑1 (MQ=255)
                                                                                     cGATACGGTGCGGGAGGTGGCCAGCCCGGCGCGGGGGGAGCCAGCGGTGTCGTATGACCCGCGTGACCCAGCGTCGAAAGATCCGTTCCTGCTGCCGGCGCTGACCGCCAACAAGGGGCGGTTCCCTGGTGCTGTACTGGTACGGCCGc                                                               <  6:266253/149‑1 (MQ=255)
                                                                                       aTACGGTGCGGGAGGTGGCCAGCCCGGCGCCGGTGGAGCCAGCGGTGTCGTATGACCCGCGTGACCCAGCGTCGAAAGATCCGTTCCTGCTGCCGGCGCTGACCGCCAACAAGGTGCGGTTCCCTGGTGCTGGACTGGTACGGCCGCAg                                                             >  6:48478/1‑149 (MQ=255)
                                                                                          cGGTGCGGGAGGTGGCCAGCCCGGCGCCGGTGGAGCCAGCGGTGTCGTATGACCCGCGTGACCCAGCGTCGAAAGATCCGTTCCTGCTGCCGGCGCTGACCGCCAACAAGGTGCGGTTCCCTGGTGCTGGACTGGTACGGCCGCAGGc                                                           <  6:198973/148‑1 (MQ=255)
                                                                                               cGGGAGGTGGCCAGCCCGGCGCCGGTGGAGCCAGCGGTGTCGTATGACCCGCGTGACCCAGCGTCGAAAGATCCGTTCCTGCTGCCGGCGCTGACCGCCAACAAGGTGCGGTTCCCTGGTGCTGGACTGGTACGGCCGCAGGccgagcc                                                     <  7:178722/149‑1 (MQ=255)
                                                                                                   aGGTGGCCAGCCCGGCGCCGGTGGAGCCAGCGGTGTCGTATGACCCGCGTGACCCAGCGTCGAAAGATCCGTTCCTGCTGCCGGCGCTGACCGCCAACAAGGTGCGGTTCCCTGGTGCTGGACTGGTACGGCCGCAGGccgagcccgag                                                 <  8:89137/149‑1 (MQ=255)
                                                                                                        gCCAGCCCGGCGCCGGTGGAGCCAGCGGTGTCGTATGACCCGCGCGACCCAGCGTCGAAAGATCCGTTCCTGCTGCCGGCGCTGACCGCCAACAAGGTGCGGTTCCCTGGTGCTGGACTGGTACGGCCGCAGGCCGAGCCCGAgccgcc                                            >  3:242496/1‑149 (MQ=255)
                                                                                                        gCCAGCCCGGCGCCGGTGGAGCCAGCGGTGTCGTATGACCCGCGCGACCCAGCGTCGAAAGATCCGTTCCTGCTGCCGGCGCTGACCGCCAACAAGGTGCGGTTCCCTGGTGCTGGACTGGTACGGCCGCAGGCCGAGCCCGAgccgcc                                            >  3:238224/1‑149 (MQ=255)
                                                                                                               cGGCGCCGGTGGAGCCAGCGGTGTCGTATGACCCGCGTGACCCAGCGTCGAAAGATCCGTTCCTGCTGCCGGCGCTGACCGCCAACAAGGTGCGGTTCCCTGGTGCTGGACTGGTACGGCCGCAGGCCGAGCCCGAGCCGCCGACAcc                                      >  8:182431/1‑148 (MQ=255)
                                                                                                                  cgcCGGTGGAGCCAGCGGTGTCGTATGACCCGCGTGACCCAGCGTCGAAAGATCCGTTCCTGCTGCCGGCGCTGACCGCCAACAAGGTGCGGTTCCCTGGTGCTGGACTGGTACGGCCGCAGGCCGAGCCCGAGCCGCCGACACCCaaa                                  <  5:48478/149‑1 (MQ=255)
                                                                                                                       gTGGAGCCAGCGGTGTCGTATGACCCGCGTGACCCAGCGTCTAAAGATCCGTTCCGGCTGCCGGCGCTGACCGCCAACAAGGTGCGGTTCCCTGGTGCTGTACTGGTACGGCCGCAGGCCGAGCCCGAGCCGCCGACACCCAAACTggg                             <  8:255845/149‑1 (MQ=255)
                                                                                                                             ccAGCGGTGTCGTATGACCCGCGTGACCCAGCGTCGAAAGATCCGTTCCTGCTGCCGGCGCTGACCGCCAACAAGGTGCGGTTCCCTGGTGCTGGACTGGTACGGCCGCAGGCCGAGCCCGAGCCGCCGACACCCAAACTGGGGGCtgc                       >  4:283426/1‑149 (MQ=255)
                                                                                                                                         tATGACCCGCGCGACCCAGCGTCGAAAGATCCGTTCCTGCTGCCGGCGCTTACCGCCAACAAGGTGCGGTTCCCTGGTGCTGGACTGGTACGGCCGCAGGCCGAGCCCGAGCCGCCGACACCCAAACTGGGGGCTGCGGTGCGGGAagc           <  4:242496/149‑1 (MQ=255)
                                                                                                                                         tATGACCCGCGCGACCCAGCGTCGAAAGATCCGTTCCTGCTGCCGGCGCTGACCGCCAACAAGGTGCGGTTCCCTGGTGCTGGACTGGTACGGCCGCAGGCCGAGCCCGAGCCGCCGACACCCAAACTGGGGGCTGCGGTGCGGGAagc           <  4:238224/149‑1 (MQ=255)
                                                                                                                                             aCCCGCGTGACCCAGCGTCGAAAGATCCGTTCCTGCTGCCGGCGCTGACCGCCAACAAGGTGCGGTTCCCTGGTGCTGGACTGGTACGGCCGCAGGCCGAGCCCGAGCCGCCGACACCCAAACTGGGGGCTGCGGTGCGGGAAgcagca       <  3:113810/149‑1 (MQ=255)
                                                                                                                                               ccGCGTGACCCAGCGTCGAAAGATCCGTTCCTGCTGCCGGCGCTGACCGCCAACAAGGTGCGGTTCCCTGGTGCTGGACTGGTACGGCCGCAGGCCGAGCCCGAGCCGCCGACACCCAAACTGGGGGCTGCGGTGCGGGAAGCAGCAcc     >  6:153577/1‑149 (MQ=255)
                                                                                                                                               ccGCGTGACCCAGCGTCGAAAGATCCGTTCCTGCTGCCGGCGCTGACCGCCAACAAGGTGCGGTTCCCTGGTGCTGGACTGGTACGGCCGCAGGCCGAGCCCGAGCCGCCGACACCCAAACTGGGGGCTGCGGTGCGGGAAGCAGCAcc     >  6:148948/1‑149 (MQ=255)
                                                                                                                                               ccGCGAGACCCAGCGTCGAAAGATCCGTTCCTGCTGCCGGCGCTGACCGCCAACAAGGTGCGGTTCCCTGGAGCTGGACTGGTACGGCCGCAGGCAGAGCCCGAGCCGCCGTCACACCACCTGGAGGAGGAGCTCCGGCAAGCAGCAcc     >  5:43451/1‑149 (MQ=21)
                                                                                                                                                cgcgTGACCCAGCGTCGAAAGATCCGTTCCTGCTGCCGGCGCTGACCGCCAACAAGGTGCGGTTCCCTGGTGCTGGACTGGTACGGCCGCAGGCCGAGCCCGAGCCGCCGACACCCAAACTGGGGGCTGCGGTGCGGGAAGCAGCACCg    <  7:182431/149‑1 (MQ=255)
                                                                                                                                                  cgTGACCCAGCGTCGAAAGATCCGTTCCTGCTGCCGGCGCTGACCGCCAACAAGGTGCGGTTCCCTGGTGCTGGACTGGTACGGCCGCAGGCCGAGCCCGAGCCGCCGACACCCAAACTGGGGGCTGCGGTGCGGGAAGCAGCACCGGc  <  5:153577/149‑1 (MQ=255)
                                                                                                                                                  cgTGACCCAGCGTCGAAAGATCCGTTCCTGCTGCCGGCGCTGACCGCCAACAAGGTGCGGTTCCCTGGTGCTGGACTGGTACGGCCGCAGGCCGAGCCCGAGCCGCCGACACCCAAACTGGGGGCTGCGGTGCGGGAAGCAGCACCGGc  <  5:148948/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TGGTGCGGCCACCGGCCAGCCCTGAGTCGCCGGCAGATCCTGCAGCACCGGAGGCACAACCCGAGCCGTCGGTGCCCAAGCTGGGCGATACGGTGCGGGAGGTGGCCAGCCCGGCGCCGGTGGAGCCAGCGGTGTCGTATGACCCGCGTGACCCAGCGTCGAAAGATCCGTTCCTGCTGCCGGCGCTGACCGCCAACAAGGTGCGGTTCCCTGGTGCTGGACTGGTACGGCCGCAGGCCGAGCCCGAGCCGCCGACACCCAAACTGGGGGCTGCGGTGCGGGAAGCAGCACCGGC  >  NC_002947/3443784‑3444078

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