Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,327,253 A→T 17.5% N74I (AAC→ATC) ‡ PP_3797 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,327,2530AT17.5% 87.4 / 12.8 40N74I (AAC→ATC) ‡PP_3797membrane protein
Reads supporting (aligned to +/- strand):  ref base A (19/14);  new base T (2/5);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 2.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

GTTCTTTGTCCTGATCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCGGAACCATTCAGGGCAGGACCCCA  >  NC_002947/4327107‑4327401
                                                                                                                                                  |                                                                                                                                                     
gTTCTTTGTCCTGATCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAAc                                                                                                                                                      >  7:24781/1‑148 (MQ=255)
   cTTTGTCCTGATCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGTTCAAg                                                                                                                                                   <  3:249054/148‑1 (MQ=255)
      tGTCCTGATCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGcc                                                                                                                                                <  8:24781/148‑1 (MQ=255)
        tCCTGATCGTGTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGTTCAAGGcctc                                                                                                                                              <  1:109259/148‑3 (MQ=255)
        tCCTGATCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGTTCAAGGcctc                                                                                                                                              <  6:161983/148‑3 (MQ=255)
             aTCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGTTCAAGGcctc                                                                                                                                              <  8:101568/143‑3 (MQ=255)
             aTCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGTTCAAGGcctc                                                                                                                                              >  7:101568/1‑141 (MQ=255)
             aTCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCATCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGTTCAAGGcctc                                                                                                                                              <  3:82588/143‑3 (MQ=255)
             aTCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCATCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGTTCAAGGcctc                                                                                                                                              >  4:82588/1‑141 (MQ=255)
                     tttCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑tgg                                                                                                                                  <  1:200330/146‑1 (MQ=255)
                       tCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑tggtggtg                                                                                                                             >  5:219453/1‑149 (MQ=255)
                                  cTCGCCGTGGGCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGa                                                                                                                   <  7:288444/148‑1 (MQ=255)
                                   tCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGagc                                                                                                                 >  5:167726/1‑149 (MQ=255)
                                                 cagGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGtt                                                                                                    <  6:219453/148‑1 (MQ=255)
                                                 cagGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTc                                                                                                   >  3:62674/1‑149 (MQ=255)
                                                               gCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTa                                                                                      >  3:64297/1‑148 (MQ=255)
                                                               gCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTa                                                                                      >  7:145074/1‑148 (MQ=255)
                                                                    ccGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATAAGCTGTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGc                                                                                 >  8:90872/1‑148 (MQ=255)
                                                                     cGATCCGGCCGCACGCGCCCACAAGCCAGACGCGGACACCGAGAGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGCGCCGGCCA‑TGGTGGTGGAGGGCATGAGCCAGCCGGATGGTTCGCCCTTCGCGGGaactgcc                                                                                >  2:295204/1‑143 (MQ=255)
                                                                     cGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCa                                                                               <  4:62674/149‑1 (MQ=255)
                                                                                       caacaaGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAAt                                                                             >  8:235076/1‑133 (MQ=255)
                                                                                       caacaaGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGAct                                                             >  1:112682/1‑148 (MQ=255)
                                                                                           gaGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCAc                                                         <  5:68448/148‑1 (MQ=255)
                                                                                             gCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCAc                                                         >  6:102547/1‑147 (MQ=255)
                                                                                              ccAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAg                                                      <  8:145074/149‑1 (MQ=255)
                                                                                                      ggACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGAt                                                                                                        >  7:187254/1‑91 (MQ=255)
                                                                                                      ggACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGAt                                                                                                        <  8:187254/91‑1 (MQ=255)
                                                                                                      ggACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGc                                              <  6:167726/149‑1 (MQ=255)
                                                                                                      ggACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGc                                              <  8:6270/149‑1 (MQ=255)
                                                                                                          acCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCtgta                                          <  2:112682/149‑1 (MQ=255)
                                                                                                                   aTGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAAacc                                 >  7:68273/1‑149 (MQ=255)
                                                                                                                         tttCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTacac                           >  2:289850/1‑149 (MQ=255)
                                                                                                                           tCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGAACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCg                         >  8:296123/1‑149 (MQ=255)
                                                                                                                             agcagGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCGGa                       <  7:90872/149‑1 (MQ=255)
                                                                                                                                       ggAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCGGAACCATTCAgg             >  7:125988/1‑149 (MQ=255)
                                                                                                                                        gAGGCCTTGAACAAGGCCGAGCCGGACAGTGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCGGAACCATTCAgg             >  1:179362/1‑149 (MQ=255)
                                                                                                                                          ggCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCGGAACCATTCAGGGCa          <  2:179362/149‑1 (MQ=255)
                                                                                                                                              ttGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCGGAACCATTCAGGGCAGGAc      >  2:186990/1‑149 (MQ=255)
                                                                                                                                                gAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCGGAACCATTCa               <  5:212946/138‑1 (MQ=255)
                                                                                                                                                gAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCGGAACCATTCa               >  6:212946/1‑138 (MQ=255)
                                                                                                                                                  aCAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCGGAACCATTCAGGGCAGGACCCCa  <  4:64297/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GTTCTTTGTCCTGATCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACA‑TGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCGGAACCATTCAGGGCAGGACCCCA  >  NC_002947/4327107‑4327401

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