Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,928,943 C→T 9.5% V160V (GTC→GTT)  PP_3464 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,928,9430CT 9.5% 115.5 / 3.0 43V160V (GTC→GTT) PP_3464hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (22/16);  new base T (2/2);  total (24/19)
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.95e-01

TCGAAGACCTGTACCTGGGCTGGCGCTCGGCCAACCTGTTCCCGGCGCTGGGTGAAGATGGCGTGGACATTTCCGGTGGCCGCCAGGTCGTCACCCTGGGTGATGGTTTCCTGATTCAGGGTGACCCGGTCAACCTGGGCAAGGTCGACCTGGGGGCCAACTTCGACCGCGGCGGTGCCTACTACCTGGCCGCGCGCAAGGCCTTCGACCGCACCGCCGTGCTGCGCCTGGGCGGCAACCAGGGCTGGCGCGGTGACCTGATGTGGTTGAAGTCCGACAACCATTACCAGGC  >  NC_002947/3928798‑3929089
                                                                                                                                                 |                                                                                                                                                  
tCGAAGACCTGTACCTGGGCTGGCGCTCGGCCAACCTGTTCCCGGCGCTGGGTGAAGATGGCGTGGACATTTCCGGTGGCCGCCAGGTCGTCACCCTGGGTGATGGTTTCCTGATTCAGGGTGACCCGGTCAACCTGGGCAAGGTCGAc                                                                                                                                                 <  6:215429/149‑1 (MQ=255)
      aCCTGTACCTGGGCTGGCGCTCGGCCAACCTGTTCCCGGCGCTGGGTGAAGATGGCGTGGACATTTCCGGTGGCCGCCAGGTCGTCACCCTGGGTGATGGTTTCCTGATTCAGGGTGACCCGGTCAACCTGGGCAAGGTCGACCTgggg                                                                                                                                           <  4:101652/149‑1 (MQ=255)
             ccTGGGCTGGCGCTCGGCCAACCTGTTCCCGGCGCTGGGTGAAGATGGCGTGGACATTTCCGGTGGCCGCCAGGTCGTCACCCTGGGTGATGGTTTCCTGATTCAGGGTGACCCGGTCAACCTGGGCAAGGTCGACCTGGGGGCCAACt                                                                                                                                    >  3:48323/1‑149 (MQ=255)
                        gcTCGGCCAACCTGTTCCCGGCGCTGGGTGAAGATGGCGTGGACATTTCCGGTGGCCGCCAGGTCGTCACCCTGGGTGATGGTTTCCTGATTCAGGGTGACCCGGTCAACCTGGGCAAGGTCGACCTGGGGGCCAACTTCGAccatgtg                                                                                                                         >  6:59557/1‑144 (MQ=255)
                        gcTCGGCCAACCTGTTCCCGGCGCTGGGTGAAGATGGCGTGGACATTTCCGGTGGCCGCCAGGTCGTCACCCTGGGTGATGGTTTCCTGATTCAGGGTGACCCGGTCAACCTGGGCAAGGTCGACCTGGGGGCCAACTTCGACCgcggc                                                                                                                         >  1:125963/1‑149 (MQ=255)
                        gcTCGGCCAACCTGTTCCCGGCGCTGGGTGAAGATGGCGTGGACATTTCCGGTGGCCGCCAGGTCGTCACCCTGGGTGATGGTTTCCTGATTCAGGGTGACCCGGTCAACCTGGGCAAGGTCGACCTGGGGGCCAACTTCGACCgcggc                                                                                                                         >  8:275350/1‑149 (MQ=255)
                            ggCCAACCTGTTCCCGGCGCTGGGTGAAGATGGCGTGGACATTTCCGGTGGACGCCAGGTCGTCACCCTGGGTGATGGTTTCCTGATTCAGGGTGACCCGGTCAACCTGGGCAAGGTCGACCTGGGGGCCAACTTCGACCGCGGCGGTg                                                                                                                     >  8:26471/1‑149 (MQ=255)
                             gCCAACCTGTTCCCGGCGCTGGGTGAAGATGGCGTGGACATTTCCGGTGGCCGCCAGGTCGTCACCCTGGGTGATGGTTTCCTGATTCAGGGTGACCCGGTCAACCTGGGCAAGGTCGACCTGGGGGCCAACTTCGACCGCGGCGGTGc                                                                                                                    <  5:49962/149‑1 (MQ=255)
                               cAACCTGTTCCCGGCGCTGGGTGAAGATGGCGTGGACATTTCCGGTGGCCGCCAGGTCGTCACCCTGGGTGATGGTTTCCTGATTCAGGGTGACCCGGTCAACCTGGGCAAGGTCGACCTGGGGGCCAACTTCGACCGCGGCGGTGCct                                                                                                                  <  4:138574/149‑1 (MQ=255)
                                   cTGTTCCCGGCGCTGGGTGAAGATGGCGTGGACATTTCCGGTGGCCGCCAGGTCGTCACCCTGGGTGATGGTTTCCTGATTCAGGGTGACCCGGTCAACCTGGGCAAGGTCGACCTGGGGGCCAACTTCGACCGCGGCGGTGCctacta                                                                                                              >  8:293234/1‑149 (MQ=255)
                                      ttCCCGGCGCTGGGTGAAGATGGCGTGGACATTTCCGGTGGCCGCCAGGTCGTCACCCTGGGTGATGGTTTCCTGATTCAGGGTGACCCGGTCAACCTGGGCAAGGTCGACCTGGGGGCCAACTTCGACCGCGGCGGTGCCtactac                                                                                                             <  4:48323/147‑1 (MQ=255)
                                            gcgcTGGGTGAAGATGGCGTGGACATTTCCGGTGGCCGCCAGGTCGTCACCCTGGGTGATGGTTTCCTGATTCAGGGTGACCCGGTCAACCTGGGCAAGGTCGACCTGGGGGCCAACTTCGACCGCGGCGGTGCCTACTACCTGGCcgc                                                                                                     <  3:24199/149‑1 (MQ=255)
                                                   gTGAAGATGGCGTGGACATTTCCGGTGGCCGCCAGGTCGTCACCCTGGGTGATGGTTTCCGGATTCAGGGTGACCCGGTCAACCTGGGCAAGGTCGACCTGGGGGCCAACTTCGACCGCGGCGGTGCCTACTACCTGGCCGCGCGCAAg                                                                                              >  8:123812/1‑149 (MQ=255)
                                                    tGAAGATGGCGTGGACATTTCCGGTGGCCGCCAGGTCGTCACCCTGGGTGATGGTTTCCTGATTCAGGGTGACCCGGTCAACCTGGGCAAGGTCGACCTGGGGGCCAACTTCGACCGCGGCGGTGCCTACTACCTGGCCGCGCGCAAgg                                                                                             <  8:166629/149‑1 (MQ=255)
                                                        gATGGCGTGGACATTTCCGGTGGCCGCCAGGTCGTCACCCTGGGTGATGGTTTCCTGATTCAGGGTGACCCGGTCAACCTGGGCAAGGTCGACCTGGGGGCCAACTTCGACCGCGGCGGTGCCTACTACCTGGCCGCGCGCAAGGCCtt                                                                                         >  6:2567/1‑149 (MQ=255)
                                                               tGGACATTTCCGGTGGCCGCCAGGTCGTCACCCTGGGTGATGGTTTCCTGATTCAGGGTGACCCGGTCAACCTGGGCAAGGTCGACCTGGGGGCCAACTTCGACCGCGGCGGTGCCTACTACCTGGCCGCGCGCAAGGCCTTCGaccgc                                                                                  <  7:293234/149‑1 (MQ=255)
                                                                  aCATTTCCGGTGGCCGCCAGGTCGTCACCCTGGGTGATGGTTTCCTGATTCAGGGTGACCCGGTCAACCTGGGCAAGGTCGACCTGGGGGCCAACTTCGACCGCGGCGGTGCCTACTACCt                                                                                                           >  1:213597/1‑121 (MQ=255)
                                                                  aCATTTCCGGTGGCCGCCAGGTCGTCACCCTGGGTGATGGTTTCCTGATTCAGGGTGACCCGGTCAACCTGGGCAAGGTCGACCTGGGGGCCAACTTCGACCGCGGCGGTGCCTACTACCt                                                                                                           <  2:213597/121‑1 (MQ=255)
                                                                         cGGTGGCCGCCAGGTCGTCACCCTGGGTGATGGTTTCCTGATTCAGGGTGACCCGGTCAACCTGGGCAAGGTCGACCTGGGGGCCAACTTCGACCGCGGCGGTGCCTACTACCTGGCCGCGCGCAAGGCCTTCGACCGCACCGCCGtgc                                                                        >  8:253542/1‑149 (MQ=255)
                                                                              gccgccAGGTCGTCACCCTGGGTGATGGTTTCCTGATTCAGGGTGACCCGGTCAACCTGGGCAAGGTCGACCTGGGGGCCAACTTCGACCGCGGCGGTGCCTACTACCTGGCCGCGCGCAAGGCCTTCGACCGCACCGCCGTGCTgcgc                                                                   >  5:299386/1‑149 (MQ=255)
                                                                                      gtcgtcACCCTGGGTGATGGTTTCCTGATTCAGGGTGACCCGGTCAACCTGGGCAAGGTCGACCTGGGGGCCAACTTCGACCGCGGCGGTGCCTACTACCTGGCCGCGCGCAAGGCCTTCGACCGCACCGCCGTGCTGCGCCTGggcgg                                                           <  7:26471/149‑1 (MQ=255)
                                                                                         gtcACCCTGGGTGATGGTTTCCTGATTCAGGGTGACCCGGTCAACCTGGGCAAGGTCGACCTGGGGGCCAACTTCGACCGCGGCGGTGCCTACTACCTgg                                                                                                         >  8:78774/1‑100 (MQ=255)
                                                                                         gtcACCCTGGGTGATGGTTTCCTGATTCAGGGTGACCCGGTCAACCTGGGCAAGGTCGACCTGGGGGCCAACTTCGACCGCGGCGGTGCCTACTACCTgg                                                                                                         <  7:78774/100‑1 (MQ=255)
                                                                                                 gggTGATGGTTTCCTGATTCAGGGTGACCCGGTCAACCTGGGCAAGGTCGACCTGGGGGCCAACTTCGACCGCGGCGGTGCCTACTACCTGGCCGCGCGCAAg                                                                                              <  8:183019/103‑1 (MQ=255)
                                                                                                 gggTGATGGTTTCCTGATTCAGGGTGACCCGGTCAACCTGGGCAAGGTCGACCTGGGGGCCAACTTCGACCGCGGCGGTGCCTACTACCTGGCCGCGCGCAAg                                                                                              >  7:183019/1‑103 (MQ=255)
                                                                                                  ggTGATGGTTTCCTGATTCAGGGTGACCCGGTCAACCTGGGCAAGGTCGACCTGGGGGCCAACTTCGACCGCGGCGGTGCCTACTACCTGGCCGCGCGCAAGGCCTTCGACCGCACCGCCGTGCTGCGCCTGGGCGGCAACCAGGGCTg                                               <  7:253542/149‑1 (MQ=255)
                                                                                                    tGATGGTTTCCTGATTCAGGGTGACCCGGTCAACCTGGGCAAGGTCGACCTGGGGGCCAACTTCGACCGCGGCGGTGCCTACTACCTGGCCGCGCGCAAGGCCTTCGACCGCACCGCCGTGCTGCGCCTGGGCGGCAACCAGGGCTGgc                                             >  4:43669/1‑149 (MQ=255)
                                                                                                             ccTGATTCAGGGTGACCCGGTCAACCTGGGCAAGGTCGACCTGGGGGCCAACTTCGACCGCGGCGGTGCCTACTACCTGGCCGCGCGCAAGGCCTTCGACCGCACCGCCGTGCTGCGCCTGGGCGGCAACCAGGGCTGGCGCGGTGAc                                     <  3:43669/148‑1 (MQ=255)
                                                                                                              cTGATTCAGGGTGACCCGGTCAACCTGGGCAAGGTGGACCTGGGGGCCAACTTGGACCTCGGCGGTGCCTACTACCTGGCCGCGCGCTAGGCCTTCGACCGCACCGCCGTGCTGCGCCTGGGCGTCAACCAGGGCTTGCGCGGTGAcc                                    <  4:274254/148‑1 (MQ=255)
                                                                                                              cTGATTCAGGGTGACCCGGTCAACCTGGGCAAGGTCGACCTGGGGGCCAACTTCGACCGCGGCGGTGCCTACTACCTGGCCGCGCGCAAGGCCTTCGACCGCACCGCCGTGCTGCGCCTGGGCGGCAACCAGGGCTGGCGCGGTGACCt                                   >  4:263241/1‑149 (MQ=255)
                                                                                                                         cccccaggTCGACCTTGCCCAGGTTGACCGGGGGGCCAACTTCGACCGCGGCGGTGCCTACTACCTGGCCGCGCGCAAGGCCTTCGACCGCACCGCCGTGCTGCGCCTGGGCGGCAACCAGGGCTGGCGCGGTGACCTGATGTGGTTg                         >  6:313648/7‑148 (MQ=255)
                                                                                                                            ccaggTCGACCTTGCCCAGGTTGACCTGGGGGCCAACTTCGACCGCGGCGGTGCCTACTACCTGGCCGCGCGCAAGGCCTTCGACCGCACCGCCGTGCTGCGCCTGGGCGGCAACCAGGGCTGGCGCGGTGACCTGATGTGGTTGaa                       <  7:238561/144‑1 (MQ=255)
                                                                                                                            ccaggTCGACCTTGCCCAGGTTGACCTGGGGGCCAACTTCGACCGCGGCGGTGCCTACTACCTGGCCGCGCGCAAGGCCTTCGACCGCACCGCCGTGCTGCGCCTGGGCGGCAACCAGGGCTGGCGCGGTGACCTGATGTGGTTGaa                       >  8:238561/4‑147 (MQ=255)
                                                                                                                              cGGTCAACCTGGGCAAGGTCGACCTGGGGGCCAACTTCGACCGCGGCGGTGCCTACTACCTGGCCGCGCGCAAGGCCTTCGACCGCACCGCCGTGCTGCGCCTGGGCGGCAACCAGGGCTGGCGCGGTGACCTGATGTGGTTGAAGTcc                   >  4:213091/1‑149 (MQ=255)
                                                                                                                              aggTCGACCTTGCCCAGGTTGACCGGGGGGCCAACTTCGACCGCGGCGGTGCCTACTACCTGGCCGCGCGCAAGGCCTTCGACCGCACCGCCGTGCTGCGCCTGGGCGGCAACCAGGGCTGGCGCGGTGACCTGATGTGGTTGAAGTcc                   <  5:313648/148‑1 (MQ=255)
                                                                                                                                     ccTGGGCAAGGTCGACCTGGGGGCCAACTTCGACCGCGGCGGTGCCTACTACCTGGCCGCGCGCAAGGCCTTCGACCGCACCGCCGTGCTGCGCCTGGGCGGCAACCAGGGCTGGCACGGTGACCTGATGTGGTTGAAGTCCGACAAc             >  4:215685/1‑148 (MQ=255)
                                                                                                                                          gCAAGGTCGACCTGGGGGCCAACTTCGACCGCGGCGGTGCCTACTACCTGGCCGCGCGCAAGGCCTTCGACCGCACCGCCGTGCTGCGCCTGGGCGGCAACCAGGGCTGGCGCGGTGACCTGATGTGGTTGAAGTCCGACAACCATTAc       >  6:199521/1‑149 (MQ=255)
                                                                                                                                             aGGTCGACCTGGGGGCCAACTTCGACCGCGGCGGTGCCTACTACCTGGCCGCGCGCAAGGCCTTCGACCGCACCGCCGTGCTGCGCCTGGGCGGCAACCAGGGCTGGCGCGGTGACCTGATGTGGTTGAAGTCCGACAACCATTACCAg    >  2:315049/1‑149 (MQ=255)
                                                                                                                                               gTCGACCTGGGGGCCAACTTCGACCGCGGCGGTGCCTACTACCTGGCCGCGCGCAAGGCCTTCGACCGCACCGCCGTGCTGCGCCTGGGCGGCAACCAGGGCTGGCGCGGTGACCTGATGTGGt                           <  2:88043/124‑1 (MQ=255)
                                                                                                                                               gTCGACCTGGGGGCCAACTTCGACCGCGGCGGTGCCTACTACCTGGCCGCGCGCAAGGCCTTCGACCGCACCGCCGTGCTGCGCCTGGGCGGCAACCAGGGCTGGCGCGGTGACCTGATGTGGt                           >  1:88043/1‑124 (MQ=255)
                                                                                                                                               gTCGACCTGGGGGCCAACTTCGACCGCGGCGGTGCCTACTACCTGGCCGCGCGCAAGGCCTTCGACCGCACCGCCGTGCTGCGCCTGGGCGGCAACCAGGGCTGGCGCGGTGACCTGATGTGGTTGAAGTCCGACAACCATTACCAGGc  >  7:145590/1‑149 (MQ=255)
                                                                                                                                                tCGACCTGGGGGCCAACTTCGACCGCGGCGGTGCCTACTACCTGGCCGCGCGCAAGGCCTTCGACCGCACCGCCGTGCTGCGCCTGGGCGGCAACCAGGGCTGGCGCGGTGACCTGATGTGGTTGAAGTCCGACAAcc            <  7:5740/138‑1 (MQ=255)
                                                                                                                                                tCGACCTGGGGGCCAACTTCGACCGCGGCGGTGCCTACTACCTGGCCGCGCGCAAGGCCTTCGACCGCACCGCCGTGCTGCGCCTGGGCGGCAACCAGGGCTGGCGCGGTGACCTGATGTGGTTGAAGTCCGACAAcc            >  8:5740/1‑138 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TCGAAGACCTGTACCTGGGCTGGCGCTCGGCCAACCTGTTCCCGGCGCTGGGTGAAGATGGCGTGGACATTTCCGGTGGCCGCCAGGTCGTCACCCTGGGTGATGGTTTCCTGATTCAGGGTGACCCGGTCAACCTGGGCAAGGTCGACCTGGGGGCCAACTTCGACCGCGGCGGTGCCTACTACCTGGCCGCGCGCAAGGCCTTCGACCGCACCGCCGTGCTGCGCCTGGGCGGCAACCAGGGCTGGCGCGGTGACCTGATGTGGTTGAAGTCCGACAACCATTACCAGGC  >  NC_002947/3928798‑3929089

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