Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,033,219 C→T 9.5% T286I (ACC→ATC)  PP_3555 → transmembrane sensor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,033,2190CT 9.5% 111.7 / 2.6 42T286I (ACC→ATC) PP_3555transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base C (20/18);  new base T (2/2);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.64e-01

GCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGGCGCAGCGGCAAGCTGGTGCTCGACAACCTGAGCCTGGAACAGGCGCTGCCGGTGATCAACCGCTACCTCGACGCCCCGCTGCTACTGGCCGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTGCCGGTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTGCTCAATCGCATCTCGCCGCCGCCCGGCAAGG  >  NC_002947/4033073‑4033365
                                                                                                                                                  |                                                                                                                                                  
gCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGGCGCAGCGGCAAGCTGGTGCTCGACAACCTGAGCCTGGAACAGGCGCTGCCGGTGATCAACCGCTACCTCGACGCCCCGCTGCTACTGGCCGATGCCAGCACCg                                                                                                                                                  <  3:175203/149‑1 (MQ=255)
   gCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGGCGCAGCGGCAAGCTGGTGCTCGACAACCTGAGCCTGGAACAGGCGCTGCCGGTGATCAACCGCTACCTCGACGCCCCGCTGCTACTGGCCGATGCCAGCACCGGcc                                                                                                                                               >  2:193100/1‑149 (MQ=255)
            gCAAAGCGACGACTATGGCAACAGCCTGGCATGGCGCAGCGGCAAGCTGGTGCTCGACAACCTGAGCCTGGAACAGGCGCTGCCGGTGATCAACCGCTACCTCGACGCCCCGCTGCTACTGGCCGATGCCAGCACCGGCCGCATTCGc                                                                                                                                       >  6:224334/1‑148 (MQ=255)
              aaaGCGACGACTATTGCATCAGCCTGGCATGGCGCAGCGGCAAGCTGGTGCTCGACAACCTGAGCCTGGAACAGGCGCTGCCGGTGATCAACCGCTACCTCGACGCCCCGCTGCTACTGGCCGATGCCAGCACCGGc                                                                                                                                                <  2:288995/137‑1 (MQ=255)
              aaaGCGACGACTATGGCAACAGCCTGGCATGGCGCACCGGCAAGCTGGTGCTCGACAACCTGAGCCTGGAACAGGCGCTTCCGCTGATCAACCGCTACCTCGACGCCCCGCTGCTACTGGCCGATGCCAGCACCgaa                                                                                                                                                >  1:288995/1‑135 (MQ=255)
                    acgacTATGGCAACAGCCTGGCATGGCGCAGCGGCAAGCTGGTGCTCGACAACCTGAGCCTGGAACAGGCGCTGCCGGTGATCAACCGCTACCTCGACGCCCCGCTGCTACTGGCCGATGCCAGCACCGGCCGCATTCGCATCAGCgg                                                                                                                               >  7:123309/1‑148 (MQ=255)
                      gacTATGGCAACAGCCTGGCATGGTGCAGCGGCAAGCTGGTGCTCGACAACCTGAGCCTGGAACAGGCGCTGCCGGTGATCAACCGCTACCTCGACGCCCCGCTGCTACTGGCCGATGCCAGCACCGGCCGCATTCGCATCAGCGGCa                                                                                                                             >  7:213474/1‑148 (MQ=255)
                      gacTATGGCAACAGCCTGGCATGGCGCAGCGGCAAGCTGGTGCTCGACAACCTGAGCCTGGAACAGGCGCTGCCGGTGATCAACCGCTACCTCGACGCCCCGCTGCTACTGGCCGATGCCAGCACCGGCCGCATTCGCATCAGCGGCa                                                                                                                             >  3:181970/1‑148 (MQ=255)
                           tGGCAACAGCCTGGCATGGCGCAGCGGCAAGCTGGTGCTCGACAACCTGAGCCTGGAACAGGCGCTGCCGGTGATCAACCGCTACCTCGACGCCCCGCTGCTACTGGCCGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTac                                                                                                                        <  8:42174/148‑1 (MQ=255)
                                 cAGCCTGGCATGGCGCAGCGGCAAGCTGGTGCTCGACAACCTGAGCCTGGAACAGGCGCTGCCGGTGATCAACCGCTACCTCGACGCCCCGCTGCTACTGGCCGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACAccc                                                                                                                 <  8:217368/149‑1 (MQ=255)
                                    ccTGGCATGGCGCAGCGGCAAGCTGGTGCTCGACAACCTGAGCCTGGAACAGGCGCTGCCGGTGATCAACCGCTACCTCGACGCCCCGCTGCTACTGGCCGATGCCAGCACCGGCCGCATTCGCATCAGCgg                                                                                                                               <  4:15935/132‑1 (MQ=255)
                                    ccTGGCATGGCGCAGCGGCAAGCTGGTGCTCGACAACCTGAGCCTGGAACAGGCGCTGCCGGTGATCAACCGCTACCTCGACGCCCCGCTGCTACTGGCCGATGCCAGCACCGGCCGCATTCGCATCAGCgg                                                                                                                               >  3:15935/1‑132 (MQ=255)
                                    ccTGGCATGGCGCAGCGGCAAGCTGGTGCTCGACAACCTGAGCCTGGAACAGGCGCTGCCGGTGATCAACCGCTACCTCGACGCCCCGCTGCTACTGGCCGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCgg                                                                                                               >  3:260358/1‑148 (MQ=255)
                                        gCATGGCGCAGCGGCAAGCTGGTGCTCGACAACCTGAGCCTGGAACAGGCGCTGCCGGTGATCAACCGCTACCTCGACGCCCCGCTGCTACTGGCCGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGaa                                                                                                            <  5:224334/147‑1 (MQ=255)
                                        gCATGGCGCAGCGGCAAGCTGGTGCTCGACAACCTGAGCCTGGAACAGGCGCTGCCGGTGATCAACCGCTACCTCGACGCCCCGCTGCTACTGGCCGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGt                                                                                                          >  8:168796/1‑149 (MQ=255)
                                         cATGGCGCAGCGGCAAGCTGGTGCTCGACAACCTGAGCCTGGAACAGGCGCTGCCGGTGATCAACCGCTACCTCGACGCCCCGCTGCTACTGGCCGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTg                                                                                                         <  7:168796/149‑1 (MQ=255)
                                             gcgcAGCGGCAAGCTGGTGCTCGACAACCTGAGCCTGGAACAGGCGCTGCCGGTGATCAACCGCTACCTCGACGCCCCGCTGCTACTGGCCGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGgc                                                                                                     >  2:21143/1‑149 (MQ=255)
                                               gcAGCGGCAAGCTGGTGCTCGACAACCTGAGCCTGGAACAGGCGCTGCCGGTGATCAACCGCTACCTCGACGCCCCGCTGCTACTGGCCGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGgcgc                                                                                                   <  4:260358/149‑1 (MQ=255)
                                                cAGCGGCAAGCTAGTGCTCGACAACCTGAGTCTGGAACAGGCGCTGCCGGTGATCAACCGCTACCTCGACGCCCCGCTGCTACTGGCCGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGcc                                                                                                  <  4:68085/149‑1 (MQ=255)
                                                             gTGCTCGACAACCTGAGCCTGGAACAGGCGCTGCCGGTGATCAACCGCTACCTCGACGCCCCGCTGCTACTGGCCGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAAcc                                                                                      <  1:193100/148‑1 (MQ=255)
                                                                cTCGACAACCTGAGCCTGGAACAGGCGCTGCCGGTGATCAACCGCTACCTCGACGCCCCGCTGCTACTGGCCGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGc                                                                                   <  2:53399/148‑1 (MQ=255)
                                                                 tCGACAACCTGAGCCTGGAACAGGCGCTGCCGGTGATCAACCGCTACCTCGACGCCCCGCTGCTACTGGCCGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGcc                                                                                  >  3:31567/1‑148 (MQ=255)
                                                                    aCAACCTGAGCCTGGAACAGGCGCTGCCGGTGATCAACCGCTACCTCGACGCCCCGCTGCTACTGGCCGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCaaa                                                                              <  3:56278/149‑1 (MQ=255)
                                                                        ccTGAGCCTGGAACAGGCGCTGCCGGTGATCAACCGCTACCTCGACGCCCCGCTGCTACTGGCCGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGaa                                                                                                            <  4:51625/115‑1 (MQ=255)
                                                                        ccTGAGCCTGGAACAGGCGCTGCCGGTGATCAACCGCTACCTCGACGCCCCGCTGCTACTGGCCGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGaa                                                                                                            >  3:51625/1‑115 (MQ=255)
                                                                              ccTGGAACAGGCGCTGCCGGTGATCAACCGCTACCTCGACGCCCCGCTGCTACTGGCCGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTGCCgg                                                                    >  2:191662/1‑149 (MQ=255)
                                                                                   aaCAGGCGCTGCCGGTGATCAACCGCTACCTCGCCGCCCCGCTGCTACTGGCCGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTGCCGGTCTAc                                                               <  8:123309/149‑1 (MQ=255)
                                                                                     cAGGCGCTGCCGGTGATCAACCGCTACCTCGACGCCCCGCTGCTACTGGCCGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCaacaac                                                                                       >  4:106022/1‑123 (MQ=255)
                                                                                     cAGGCGCTGCCGGTGATCAACCGCTACCTCGACGCCCCGCTGCTACTGGCCGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCaacaac                                                                                       <  3:106022/123‑1 (MQ=255)
                                                                                                 gTGATCAACCGCTACCTCGACGCCCCGCTGCTACTGGCCGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTGCCGGTCTACCTGACCCgcagc                                                   >  6:335296/1‑147 (MQ=255)
                                                                                                                 tCGACGCCCCGCTGCTACTGGCCGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTGCCGGTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTg                                 >  7:37193/1‑149 (MQ=255)
                                                                                                                   gACGCCCCGCTGCTACTGGCCGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTGCCGGTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTg                                 <  8:213474/147‑1 (MQ=255)
                                                                                                                       ccccGCTGCTACTGGCCGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTGCCGGTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTGCTCAAt                           >  5:192176/1‑149 (MQ=255)
                                                                                                                         ccGCTGCTACTGGCCGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACACGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTGCCGGTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTGCTCAATCg                         >  1:222323/1‑149 (MQ=255)
                                                                                                                          cGCTGCTACTGGCCGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTGCCGGTc                                                                  <  5:239925/107‑1 (MQ=255)
                                                                                                                          cGCTGCTACTGGCCGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTGCCGGTc                                                                  >  6:239925/1‑107 (MQ=255)
                                                                                                                              gaatgcggCCGGTGCTGGCATCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTGCCGGTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTGCTCAATCGCAtct                    >  7:9095/7‑149 (MQ=255)
                                                                                                                               aatgcggCCGGTGCTGGCATCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTGCCGGTCTACCTGACCCGCAGCAAGGAc                                             <  5:132726/118‑1 (MQ=255)
                                                                                                                               aatgcggCCGGTGCTGGCATCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTGCCGGTCTACCTGACCCGCAGCAAGGAc                                             >  6:132726/6‑123 (MQ=255)
                                                                                                                                     gCCGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTGCCGGTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTGCTCAATCGCATCTcgccg               <  5:335296/147‑1 (MQ=255)
                                                                                                                                            ctggCATCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTGCCGGTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTGCTCAATCGCATCTCGCCGCCGCCCGGc      <  8:9095/146‑1 (MQ=255)
                                                                                                                                                cACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTGCCGGTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTGCTCAATCGCATCTCGCCGCCGCCCGGCACgg  >  5:49290/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCCGCAGCTGGAGCAAAGCGACGACTATGGCAACAGCCTGGCATGGCGCAGCGGCAAGCTGGTGCTCGACAACCTGAGCCTGGAACAGGCGCTGCCGGTGATCAACCGCTACCTCGACGCCCCGCTGCTACTGGCCGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTGCCGGTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTGCTCAATCGCATCTCGCCGCCGCCCGGCAAGG  >  NC_002947/4033073‑4033365

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