Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,593,267 T→A 9.7% R101R (CGT→CGA)  PP_3171 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,593,2670TA 9.7% 128.6 / 5.7 52R101R (CGT→CGA) PP_3171hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (23/24);  new base A (3/2);  total (26/26)
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.10e-01

TGCCACGTCGCTGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCG  >  NC_002947/3593126‑3593415
                                                                                                                                             |                                                                                                                                                    
tGCCACGTCGCTGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGAcc                                                                                                                                               <  3:111708/149‑1 (MQ=255)
    aCGTCGCTGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGATACGAcctggg                                                                                                                                           <  1:158763/149‑1 (MQ=255)
         gCTGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGAcctgggcctgg                                                                                                                                      <  5:262331/149‑1 (MQ=255)
                aTCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTg                                                                                                                                <  8:218503/148‑1 (MQ=255)
                 tCAACGAATTGGCCAACGCTAGCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTgc                                                                                                                               <  6:121549/148‑1 (MQ=255)
                 tCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTgcg                                                                                                                              <  1:79408/149‑1 (MQ=255)
                      gAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAgg                                                                                                                         >  3:61453/1‑149 (MQ=255)
                      gAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGATACGAcctgggt                                                                                                                                          >  3:215071/1‑131 (MQ=255)
                      gAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGATACGAcctgggt                                                                                                                                          <  4:215071/132‑2 (MQ=255)
                           cGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGt                                                                                                                    >  4:178318/1‑149 (MQ=255)
                               aaCGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGGTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCgg                                                                                                                <  3:181313/149‑1 (MQ=255)
                                aCGCCACCCGACACATTCGTGTGTGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCggg                                                                                                               <  4:27717/149‑1 (MQ=255)
                                      cccGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGAcc                                                                                                         >  8:175551/1‑149 (MQ=255)
                                                  gtgtGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCAtt                                                                                              <  3:37925/148‑1 (MQ=255)
                                                     tgGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTgcgc                                                                                          >  1:86281/1‑149 (MQ=255)
                                                        gTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCgcggc                                                                                       <  3:178318/149‑1 (MQ=255)
                                                             ggCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTc                                                                                                     <  2:228920/130‑1 (MQ=255)
                                                             ggCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTc                                                                                                     >  1:228920/1‑130 (MQ=255)
                                                              gCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTgcgc                                                                                                                             <  2:132034/105‑1 (MQ=255)
                                                              gCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTgcgc                                                                                                                             >  1:132034/1‑105 (MQ=255)
                                                               cGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGgcgccgcgc                                                                                >  8:95195/1‑149 (MQ=255)
                                                                       tgctgcCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCCCGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCgg                                                                        >  2:157815/1‑149 (MQ=255)
                                                                              cacccATGCGCCGTTGGGGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTAGCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCg                                                                 <  7:43255/146‑1 (MQ=255)
                                                                              cAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCg                                                                 <  4:61453/149‑1 (MQ=255)
                                                                               aaCCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCgg                                                                >  2:160432/1‑149 (MQ=255)
                                                                                  cATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGc                                                               >  2:70874/1‑147 (MQ=255)
                                                                                        ccGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCCGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCgg                                                                        <  4:139760/132‑1 (MQ=255)
                                                                                        ccGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACAt                                                       >  1:261513/1‑149 (MQ=255)
                                                                                         cGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATc                                                      <  1:220553/149‑1 (MQ=255)
                                                                                              gtggtggCCTAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGTGATTTCGGCCAGGACATCGCCg                                                  <  7:39861/148‑1 (MQ=255)
                                                                                                     ccGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCTAGCTGCGt                                          <  7:175551/149‑1 (MQ=255)
                                                                                                      cGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGt                                          >  8:120553/1‑148 (MQ=255)
                                                                                                       gAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGa                                        >  3:11298/1‑149 (MQ=255)
                                                                                                          cAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTa                                     >  3:107479/1‑149 (MQ=255)
                                                                                                           aGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTAc                                    <  4:11298/149‑1 (MQ=255)
                                                                                                             tttGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGt                                          >  7:70066/1‑141 (MQ=255)
                                                                                                             tttGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGt                                          <  8:70066/141‑1 (MQ=255)
                                                                                                             ttgggCACCATAAATACATGCCACCCAGTCGGTACCGACCGGGGCCTGGGCCGTGCGCCAGGCAGGTCCGGGCCGACCGTCCGGGATTTGCGCGGCGCCGCGCCAGAGCGGGATTTCGTCCAGGACATCGCCGAGCTGCGTGATTAcc                                   <  1:157815/145‑1 (MQ=255)
                                                                                                                ggCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGc                               >  1:27063/1‑149 (MQ=255)
                                                                                                                  cACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGAtt                                      <  7:72341/140‑1 (MQ=255)
                                                                                                                  cACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGAtt                                      >  8:72341/1‑140 (MQ=255)
                                                                                                                   aCCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGTCATCGCCGAGCTGCGTGATTACCTGGCCGa                            <  1:160432/149‑1 (MQ=255)
                                                                                                                   aCCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGa                            >  6:191821/1‑149 (MQ=255)
                                                                                                                      cTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACaa                         >  3:59914/1‑149 (MQ=255)
                                                                                                                       tCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCGAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCg                                                                         <  1:202150/100‑1 (MQ=255)
                                                                                                                       tCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCg                                                                         >  2:202150/1‑100 (MQ=255)
                                                                                                                          aTACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGc                     <  8:118346/149‑1 (MQ=255)
                                                                                                                                     ccAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGt          >  3:31451/1‑149 (MQ=255)
                                                                                                                                     ccAGGTCGATACGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGt          >  5:200454/1‑149 (MQ=255)
                                                                                                                                     ccAGGTCGATACGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGt          >  8:120755/1‑149 (MQ=255)
                                                                                                                                      cAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTg         >  7:34858/1‑149 (MQ=255)
                                                                                                                                      cAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTg         >  8:241701/1‑149 (MQ=255)
                                                                                                                                             tATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCgcggcg  >  6:217632/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
TGCCACGTCGCTGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCG  >  NC_002947/3593126‑3593415

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