Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,727,861 T→C 8.3% A84A (GCA→GCG)  PP_2389 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,727,8610TC 8.3% 131.7 / 2.0 48A84A (GCA→GCG) PP_2389hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (22/22);  new base C (2/2);  total (24/24)
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.66e-01

TGCCTTGAGGCCAAGGGTGCGGGTCGCGCGCATGCACGCCGCGCAGTACCTTCCAGGCCAAGGCCTGCTGCTCGGCCAGGCTGCTGCCGGCTGCGGTTTCGAGAAACACCACCAAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGTGCGACGGTGCGCACATAGTTGGCCAGGCCCTGCGCGATATCATCGGCCAGTTGCTCGGCGCGGGCGAAGAGAATGTGGCAGGTTTCGCCGTTTACCGCCTTGCGGCCGAACAGGCAGGGGAATTCGGGGTTGGCGA  >  NC_002947/2727713‑2727997
                                                                                                                                                    |                                                                                                                                        
tGCCTTGAGGCCAAGGGTGCGGGTCGCGCGCATGCACGCCGCGCAGTACCTTCCAGGCCAAGGCCTGCTGCTCGGCCAGGCTGCTGCCGGCTGCGGTTTCGAGAAACACCACCAAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGt                                                                                                                                          <  8:257314/149‑1 (MQ=255)
  ccTTGAGGCCAAGGGTGCGGGTCGCGCGCATGCACGCCGCGCAGTACCTTCCAGGCCAAGGCCTGCTGCTCGGCCAGGCTGCTGCCGGCTGCGGTTTCGAGAAACACCACCAAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGTGc                                                                                                                                        <  6:62100/149‑1 (MQ=255)
  ccTTGAGGCCAAGGGTGCGGGTCGCGCGCATGCACGCCGCGCAGTACCTTCCAGGCCAAGGCCTGCTGCTCGGCCAGGCGGCTGCCGGCGGCGGTTTCGAGAAACACCACCAACGGGCTCCCGATGCGCTGCTGGGGTGCGATCGGTGc                                                                                                                                        <  2:40397/149‑1 (MQ=255)
   cTTGAGGCCAAGGGTGCGGGTCGCGCGCATGCACGCCGCGCAGTACCTTCCAGGCCAAGGCCTGCTGCTCGGCCAGGCTGCTGCCGGCTGCGGTTTCGAGAAACACCACCAAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGTGc                                                                                                                                        >  2:193766/1‑148 (MQ=255)
       aGGCCAAGGGTGCGGGTCGCGCGCATGCACGCCGCGCAGTACCTTCCAGGCCAAGGCCTGCTGCTCGGCCAGGCTGCTGCCGGCTGCGGTTTCGAGAAACACCACCAAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGTGCGACg                                                                                                                                    <  3:99906/148‑1 (MQ=255)
        ggCCAAGGGTGCGGGTCGCGCGCATGCACGCCGCGCAGTACCTTCCAGGCCAAGGCCTGCTGCTCGGCCAGGCTGCTGCCGGCTGCGGTTTCGAGAAACACCACCAAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGTGCGACGGt                                                                                                                                  <  1:176678/149‑1 (MQ=255)
         gCCAAGGGTGCGGGTCGCGCGCATGCACGCCGCGCAGTACCTTCCAGGCCAAGGCCTGCTGCTCGGCCAGGCTGCTGCCGGCTGCGGTTTCGAGAAACACCACCAAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGTGCGACGGTg                                                                                                                                 >  1:235636/1‑149 (MQ=255)
            aaGGGTGCGGGTCGCGCGCATGCACGCCGCGCAGTACCTTCCAGGCCAAGGCCTGCTGCTCGGCCAGGCTGCTGCCGGCTGCGGTTTCGAGAAACACCACCAAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGTGCGACGGTgcgc                                                                                                                              <  4:7233/149‑1 (MQ=255)
            aaGGGTGCGGGTCGCGCGCATGCACGCCGCGCAGTACCTTCCAGGCCAAGGCCTGCTGCTCGGCCAGGCTGCTGCCGGCTGCGGTTTCGAGAAACACCACCAAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGTGCGACGGTgcg                                                                                                                               >  7:180358/1‑148 (MQ=255)
                    gggTCGCGCGCATGCACGCCGCGCAGTACCTTCCAGGCCAAGGCCTGCTGCTCGGCCAGGCTGCTGCCGGCTGCGGTTTCGAGAAACACCACCAAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGTGCGACGGTGCGCACATAGtt                                                                                                                      <  2:263855/149‑1 (MQ=255)
                           gcgcATGCACGCCGCGCAGTACCTTCCAGGCCAAGGCCTGCTGCTCGGCCAGGCTGCTGCCGGCTGCGGTTTCGAGAAACACCACCAAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGTGCGACGGTGCGCACATAGTTGGCCAgg                                                                                                               <  1:193766/149‑1 (MQ=255)
                                   aCGCCGCGCAGTACCTTCCAGGCCAAGGCCTGCTGCTCGGCCAGGCTGCTGCCGGCTGCGGTTTCGAGAAACACCACCAAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGTGCGACGGTGCGCACATAGTTGGCCAGGc                                                                                                              >  8:236563/1‑142 (MQ=255)
                                   aCGCCGCGCAGTACCTTCCAGGCCAAGGCCTGCTGCTCGGCCAGGCTGCTGCCGGCTGCGGTTTCGAGAAACACCACCAAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGTGCGACGGTGCGCACATAGTTGGCCAGGc                                                                                                              <  7:236563/142‑1 (MQ=255)
                                   aCGCCGCGCAGTACCTTCCAGGCCAAGGCCTGCTGCTCGGCCAGGCTGCTGCCGGCTGCGGTTTCGAGAAACACCACCAAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGTGCGACGGTGCGCACATAGTTGGCCAGGCCCTgcgc                                                                                                       >  7:88544/1‑149 (MQ=255)
                                   aCGCCGCGCAGTACCTTCCAGGCCAAGGCCTGCTGCTCGGCCAGGCTGCTGCCGGCTGCGGTTTCGAGAAACACCACCAAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGTGCGACGGTGCGCACATAGTTGGCCAGGCCCTgcgc                                                                                                       >  1:348497/1‑149 (MQ=255)
                                     gccgcGCAGTACCTTCCAGGCCAAGGCCTGCTGCTCGGCCAGGCTGCTGCCGGCTGCGGTTTCGAGAAACACCACCAAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGTGCGACGGTGCGCACATAGTTGGCCAGGCCCTGCGCGa                                                                                                     <  2:235636/149‑1 (MQ=255)
                                     gccgcGCAGTACCTTCCAGGCCAAGGCCTGCTGCTCGGCCAGGCTGCTGCCGGCTGCGGTTTCGAGAAACACCACCAAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGTGCGACGGTGCGCACATAGTTGGCCAGGCCCTGCGCGa                                                                                                     >  1:380597/1‑149 (MQ=255)
                                     gccgcGCAGTACCTTCCAGGCCAAGGCCTGCTGCTCGGCCAGGCTGCTGCCGGCTGCGGTTTCGAGAAACACCACCAAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGCGCGACGGTGCGCACATAGTTGGCCAGGCCCTGCGCGa                                                                                                     >  6:268298/1‑149 (MQ=255)
                                     gccgcGCAGTACCTTCCAGGCCAAGGCCTGCTGCTCGGCCAGGCTGCTGCCGGCTGCGGTTTCGAGAAACACCACCAAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGCGCGAAGGTGCGCACATAGTTGGCCAGGCCCTGCGCGa                                                                                                     >  6:273154/1‑149 (MQ=255)
                                       cgcgCAGTACCTTCCAGGCCAAGGCCTGCTGCTCGGCCAGGCTGCTGCCGGCTGCGGTTTCGAGAAACACCACCAAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGTGCGACGGTGCGCACATAGTTGGCCAGGCCCTGCGCGata                                                                                                   <  2:380597/149‑1 (MQ=255)
                                       cgatcAGTACCATCCATGCCAAGGCCGGCTGCGAGGCCAGGCGGCTGCCGTCGGCGGTATCGAGAAACACCACCACAGGGCTGCCGATGCTCTGCTTGGGTGCGAGCTGTGCGAGGTTGCGCACATAGTTGGCCAGGCCCTGCGCGata                                                                                                   <  4:379234/145‑1 (MQ=255)
                                            aGTACCTTCCAGGCCAAGGCCTGCTGCGCGGCCAGGCTGCTGCCGGCTGCGGTTTCGAGAAACACCACCAAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGTGCGACGGTGCGCACATAGTTGGCCAGGCCCTGCGCGATatcatc                                                                                              <  8:180358/149‑1 (MQ=255)
                                                    ccAGGCCAAGGCCTGCTGCTCGGCCAGGCTGCTGCCGGCTGCGGTTTCGAGAAACACCACCAAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGTGTGACGGTGCGCACATAGTTGGCCAGGCCCTGCGCGATATCATCGGCCAg                                                                                        <  7:212177/147‑1 (MQ=255)
                                                          cAAGGCCTGCTGCTCGGCCAGGCTGCTGCCGGCTGCGGTTTCGAGAAACACCACCAAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGTGCGACGGTGCGCACATAGTTGGCCAGGCCCTGCGCGATATCATCGGCCAGTTGCTCgg                                                                                >  7:307114/1‑149 (MQ=255)
                                                                ctgctgCTCGGCCAGGCTGCTGCCGGCTGCGGTTTCGAGAAACACCACCAAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGTGCGACGGTGCGCACATAGTTGGCCAGGCCCTGCGCGATATCATCGGCCAGTTGCTCGGCGCggg                                                                          >  5:71855/1‑149 (MQ=255)
                                                                     gctCGGCCAGGCTGCTGCCGGCTGCGGTTTCGAGAAACACCACCAAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGTGCGACGGTGCGCACATAGTTGGCCAGGCCCTGCGCGATATCATCGGCCAGTTGCTCGGCGCGGGCg                                                                        >  7:42161/1‑146 (MQ=255)
                                                                            cAGGCTGCTGCCGGCTGCGGTTTCGAGAAACACCACCAAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGTGCGACGGTGCGCACATAGTTGGCCAGGCCCTGCGCGATATCATCGGCCAGTTGCTCGGCGCGGGCGAAGAGAAtgt                                                              >  8:105868/1‑149 (MQ=255)
                                                                                 tgctgcCGGCTGCGGTTTCGAGAAACACCACCAAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGTGCGACGGTGCGCACATAGTTGGCCAGGCCCTGCGCGATATCATCGGCCAGTTCCTCGGAGCGGGCGAAGAGAATGTGGCAg                                                         >  5:311705/1‑149 (MQ=255)
                                                                                       cGGCTGCGGTTTCGAGAAACACCACCAAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGTGCGACGGTGCGCACATAGTTGGCCAGGCCCTGCGCGATATCATCGGCCAGTTGCTCGGCGCGGGCGAAGAGAATGTGGCAGGTTTCg                                                   >  6:17477/1‑149 (MQ=255)
                                                                                         gCTGCGGTTTCGAGAAACACCACCAAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGTGCGACGGTGCGCACATAGTTGGCCAGGCCCTGCGCGATATCATCGGCCAGTTGCTCTGCGCGGGCGAAGAGAATGTGGCAGGTTTCGcc                                                 <  5:63558/149‑1 (MQ=255)
                                                                                           tGCGGTTTCGAGAAACACCACCAAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGTGCGACGGTGCGCACATAGTTGGCCAGGCCCTGCGCGATATCATCGGCCAGTTGCTCGGCGCGGGCGAAGAGAATGTGGCAGGTTTCGCCGt                                               >  3:347283/1‑149 (MQ=255)
                                                                                              ggTTTCGAGAAACACCACCAAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGTGCGACGGTGCGCACATAGTTGGCCAGGCCCTGCGCGATATCATCGGCCAGTTGCTCGGCGCGGGCGAAGAGAATGTGGCAGGTTTCGCCGTTTa                                            <  8:88544/149‑1 (MQ=255)
                                                                                                tttCGAGAAACACCACCAAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGTGCGACGGTGCGCACATAGTTGGCCAGGCCCTGCGCGATATCATCGGCCAGTTGCTCGGCGCGGGCGAAGAGAATGTGGCAGGTTTCGCCGTTTAcc                                          <  7:105868/149‑1 (MQ=255)
                                                                                                tttCGAGAAACACCACCAAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGTGCGACGGTGCGCACATAGTTGGCCAGGCCCTGCGCGATATCATCGGCCAGTTGCTCGGCGCGGGCGAAGAGAATGTGGCAGGTTTCGCCGTTTAcc                                          >  5:249482/1‑149 (MQ=255)
                                                                                                tttCGAGAAACACCACCAAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGTGCGACGGTGCGCACATAGTTGGCCAGGCCCTGCGCGATATCATCGGCCAGTTGCTCGGCGCGGGCGAAGAGAATGTGGCAGGTTTCGCCGTTTAcc                                          >  2:174908/1‑149 (MQ=255)
                                                                                                   cGAGAAACACCACCAAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGTGCGACGGTGCGCACATAGTTGGCCAGGCCCTGCGCGATATCATCGGCCAGTTGCTCGGCGCGGGCGAAGAGAATGTGGCAGGTTTCGCCGTTTACCGcc                                       >  7:108098/1‑149 (MQ=255)
                                                                                                       aaaCACCACCAAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGTGCGACGCTGCGCACATAGTTGGCCAGGCCCTGCGCGATATCATCGGCCAGTTGCTCGGCGCGGGCGAAGAGAATGTGGCAGGTTTCGCCGTTTACCGCCTTGc                                   >  1:295249/1‑149 (MQ=255)
                                                                                                            ccaccaAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGTGCGACGGTGCGCACATAGTTGGCCAGGCCCTGCGCGATATCATCGGCCAGTTGCTCGGCGCGGGCGAAGAGAATGTGGCAGGTTTCGCCGTTTACCGCCTTGCGGcc                               >  2:244830/1‑148 (MQ=255)
                                                                                                             caccaAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGTGCGACGGTGCGCACATAGTTGGCCAGGCCCTGCGCGATATCATCGGCCAGTTGCTCGGCGCGGGCGAAGAGAATGTGGCAGGTTTCGCCGTTTACCGCCTTGCGGCCg                              >  6:308511/1‑148 (MQ=255)
                                                                                                             caccaAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGTGCGACGGTGCGCACATAGTTGGCCAGGCCCTGCGCGATATCATCGGCCAGTTGCTCGGCGCGGGCGAAGAGAATGTGGCAGGTTTCGCCGTTTACCGCCTTGCGGCCGa                             <  5:17477/149‑1 (MQ=255)
                                                                                                                       cTGCCGATGCGCTGCTTGGGTGCGATCGGCGCGACGGTGCGCACATAGTTGGCCAGGCCCTGCGCGATATCATCGTCCAGTTGCTCGGCGCGGGCGAAGAGAATGTGGCAGGTTTCGCCGTTTACCGCCTTGCGGCCGAACAGGCAggg                   <  5:273154/149‑1 (MQ=255)
                                                                                                                       cTGCCGATGCGCTGCTTGGGTGCGATCGGCGCGACGGTGCGCACATAGTTGGCCAGGCCCTGCGCGATATCATCGGCCAGTTGCTCGGCGCGGGCGAAGAGAATGTGGCAGGTTTCGCCGTTTACCGCCTTGCGGCCGAACAGGCAggg                   <  5:268298/149‑1 (MQ=255)
                                                                                                                                cgcTGCTTGGGTGCGATCGGTGCGACGGTGCGCACATAGTTGTCCAGGCCCTGCGCGATATCATCGGCCAGTTGCTCGGCGCGGGCGAAGAGAATGTGGCAGGTTTCGCCGTTTACCGCCTTGCGGCCGAACAGGCAGGGGAATTCgg           <  4:85389/148‑1 (MQ=255)
                                                                                                                                    gctTGGGTGCGATCGGTGCGACGGTGCGCACATAGTTGGCCAGGCCCTGCGCGATATCATCGGCCAGTTGCTCGGCGCGGGCGAAGAGAATGTGGCAGGTTTCGCCGTTTACCGCCTTGCGGc                                >  5:227568/1‑123 (MQ=255)
                                                                                                                                    gctTGGGTGCGATCGGTGCGACGGTGCGCACATAGTTGGCCAGGCCCTGCGCGATATCATCGGCCAGTTGCTCGGCGCGGGCGAAGAGAATGTGGCAGGTTTCGCCGTTTACCGCCTTGCGGc                                <  6:227568/123‑1 (MQ=255)
                                                                                                                                        gggTGCGATCGGTGCGACGGTGCGCACATAGTTGGCCAGGCCCTGCGCGATATCATCGGCCAGTTGCTCGGCGCGGGCGAAGAGAATGTGGCAGGTTTCGCCGTTTACCGCCTTGCGGCCGAACAGGCAGGGGAATTCGGGGTTGGCga  <  2:348497/149‑1 (MQ=255)
                                                                                                                                         ggTGCGATCGGTGCGACGGTGCGCACATAGTTGGCCAGGCCCTGCGCGATATCATCGGCCAGTTGCTCGGCGCGGGCGAAGAGAATGTGGCAGGTTTCGCCGTTTACCGCCTTGCTGCCGAACAGGCAGGGGAATTCGGGGTTGGCga  <  8:108098/148‑1 (MQ=255)
                                                                                                                                         ggTGCGATCGGTGCGACGGTGCGCACATAGTTGGCCAGGCCCTGCGCGATATCATCGGCCAGTTGCTCGGCGCGGGCGAAGAGAATGTGGCAGGTTTCGCCGTTTACCGCCTTGCGGCCGAACCGGCAGGGGAATTCGGGGTTGGCga  <  8:42161/148‑1 (MQ=255)
                                                                                                                                                tCGGTGCGACGGTGCGCACATAGTTGGCCAGGCCCTGCGCGATATCATCGGCCAGTTGCTCGGCGCGGGCGAAGAGAATGTGGCAGGTTTCGCCGTTTACCGCCTTGCgg                                 >  7:262483/1‑110 (MQ=255)
                                                                                                                                                    |                                                                                                                                        
TGCCTTGAGGCCAAGGGTGCGGGTCGCGCGCATGCACGCCGCGCAGTACCTTCCAGGCCAAGGCCTGCTGCTCGGCCAGGCTGCTGCCGGCTGCGGTTTCGAGAAACACCACCAAAGGGCTGCCGATGCGCTGCTTGGGTGCGATCGGTGCGACGGTGCGCACATAGTTGGCCAGGCCCTGCGCGATATCATCGGCCAGTTGCTCGGCGCGGGCGAAGAGAATGTGGCAGGTTTCGCCGTTTACCGCCTTGCGGCCGAACAGGCAGGGGAATTCGGGGTTGGCGA  >  NC_002947/2727713‑2727997

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