Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,327,253 A→T 8.9% 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,2530AT 8.9% 121.7 / 2.8 45N74I (AAC→ATC) ‡PP_3797membrane protein
Reads supporting (aligned to +/- strand):  ref base A (20/21);  new base T (2/2);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

GGTTCTTTGTCCTGATCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCGGAACCATTCAGGGCAGGACCCCA  >  NC_002947/4327106‑4327401
                                                                                                                                                   |                                                                                                                                                    
ggTTCTTTGTCCTGATCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGaa                                                                                                                                                      >  5:263971/1‑148 (MQ=255)
 gaacTTTCTCCAGAACGTCTCTTTCACCGTCAAGCTCGCCGTGTTCAGTCAGGAAATCGTCTGGCTCGACGATCCCGCCGCACGCTCCAACAAGCCACACGCGGACACCGACCGCATGAGCTATCAGCAGGTACTGGATCACTTGAAc                                                                                                                                                     >  4:120355/4‑148 (MQ=255)
    cTTTGTCCTGATCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATTAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAg                                                                                                                                                  <  3:120355/148‑1 (MQ=255)
      ttGTCCTGATCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAgg                                                                                                                                                 <  4:292119/147‑1 (MQ=255)
      ttGTCCTGATCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAgg                                                                                                                                                 >  3:292119/1‑147 (MQ=255)
       tGTCCTGATCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGTTCAAGGc                                                                                                                                                <  7:280788/147‑1 (MQ=255)
         tCCTGATCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAg                                                                                                                                            <  3:38113/149‑1 (MQ=255)
          ccTGATCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAg                                                                                                                                            <  3:42157/148‑1 (MQ=255)
          ccTGATCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGc                                                                                                                                           >  2:56040/1‑149 (MQ=255)
              aTCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGACAACAGGCCAGACCCCGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAATGCCGAGCCGGa                                                                                                                                       >  1:282682/1‑149 (MQ=255)
                  tCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACAtg                                                                                                                                   <  4:153575/149‑1 (MQ=255)
                  tCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACAtg                                                                                                                                   >  3:153575/1‑149 (MQ=255)
                                   cTCGCCGTGGTCAGTCAGGAAATCGTNTGGCTCGCCNATCCNGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGTGCGCCTGag                                                                                                                  <  4:127544/149‑1 (MQ=255)
                                   cTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGTGCGCCTGag                                                                                                                  <  4:122847/149‑1 (MQ=255)
                                   cTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGag                                                                                                                  <  1:56040/149‑1 (MQ=255)
                                                gtcagGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGa                                                                                                                   >  8:207674/1‑135 (MQ=255)
                                                gtcagGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGa                                                                                                                   <  7:207674/135‑1 (MQ=255)
                                                    ggAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGa                                                                                                                                       >  8:270612/1‑111 (MQ=255)
                                                    ggAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGa                                                                                                                                       <  7:270612/111‑1 (MQ=255)
                                                    ggAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGtt                                                                                                    <  1:93565/146‑1 (MQ=255)
                                                    ggAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGtt                                                                                                    >  2:93565/1‑146 (MQ=255)
                                                           gTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTcg                                                                                           <  7:80451/148‑1 (MQ=255)
                                                                    gCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACAtggtggtg                                                                                                                             >  5:291499/1‑105 (MQ=255)
                                                                    gCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACAtggtggtg                                                                                                                             <  6:291499/105‑1 (MQ=255)
                                                                    gCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGc                                                                                 >  7:154430/1‑149 (MQ=255)
                                                                                  aCGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACcctgc                                                                   <  4:26888/149‑1 (MQ=255)
                                                                                                 aGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGcc                                                    >  5:74923/1‑149 (MQ=255)
                                                                                                 aGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGcc                                                    >  6:99562/1‑149 (MQ=255)
                                                                                                          cacCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCtgt                                           <  8:154430/149‑1 (MQ=255)
                                                                                                                     tGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCAcc                                                        >  5:229425/1‑125 (MQ=255)
                                                                                                                     tGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCAcc                                                        <  6:229425/125‑1 (MQ=255)
                                                                                                                        gCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTa                              <  6:74923/148‑1 (MQ=255)
                                                                                                                                cagGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCGGAAc                     >  1:4288/1‑149 (MQ=255)
                                                                                                                                       tGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCgg                        >  6:184309/1‑139 (MQ=255)
                                                                                                                                       tGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCgg                        <  5:184309/139‑1 (MQ=255)
                                                                                                                                       gctcggCCTTGTTCAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCGGAACCATTCAg              >  8:104480/5‑149 (MQ=255)
                                                                                                                                        ggAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCGGa                       <  7:81688/139‑1 (MQ=255)
                                                                                                                                        ggAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCGGa                       >  8:81688/1‑139 (MQ=255)
                                                                                                                                        ctcggCCTTGTTCAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACAcc                          >  7:27727/4‑136 (MQ=255)
                                                                                                                                        ctcggCCTTGTTCAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACAcc                          <  8:27727/133‑1 (MQ=255)
                                                                                                                                         gAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCCCTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTACACCCCTACACCGGAACCATTCAggg            >  6:229417/1‑149 (MQ=255)
                                                                                                                                         gAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCGGAACCATTCAg              <  6:263971/147‑1 (MQ=255)
                                                                                                                                              cTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCGGAACCATTCAGGGCAGGa       <  7:23575/149‑1 (MQ=255)
                                                                                                                                              cTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCGGAACCATTCAGGGCAGGa       >  8:169265/1‑149 (MQ=255)
                                                                                                                                              cTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCGGAACCATTCAGGGCAGGa       >  8:23575/1‑149 (MQ=255)
                                                                                                                                              cTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGAGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACACCTACACCGGAACCATTCAGGGCAGGa       >  3:86607/1‑149 (MQ=255)
                                                                                                                                                   aCAAGGACGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACTAGCCCGACGCTGTATGTAAACCCCTACACCGGAACCATTCAGGGCAGGACCCCa  <  2:282682/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGTTCTTTGTCCTGATCGTCTGTTTCACCGGCATGCTCGCCGTGGTCAGTCAGGAAATCGTCTGGCTCGCCGATCCGGCCGCACGCGCCAACAAGCCAGACGCGGACACCGAGCGCATGAGCTTTCAGCAGGTACTGGAGGCCTTGAACAAGGCCGAGCCGGACATGGTGGTGGAGCGCCTGAGCCAGCCGGATGGTTCGCACTTCGCGGTAAAAGCCAATGTCACCCTGCCTGACGGCACCAGCCCGACGCTGTATGTAAACCCCTACACCGGAACCATTCAGGGCAGGACCCCA  >  NC_002947/4327106‑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: ‑