Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,591,576 T→A 10.4% L258M (TTG→ATG) ‡ PP_3169 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,591,5760TA10.4% 126.9 / 4.4 48L258M (TTG→ATG) ‡PP_3169membrane protein
Reads supporting (aligned to +/- strand):  ref base T (22/21);  new base A (2/3);  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.57e-01

GTCTGGCGGGTGCTGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATGGTC  >  NC_002947/3591429‑3591719
                                                                                                                                                   |                                                                                                                                               
gTCTGGCGGGTGCTGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCtt                                                                                                                                                <  8:344250/149‑1 (MQ=255)
      cGGGTGCTGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCAAGCGCTg                                                                                                                                          <  5:191502/149‑1 (MQ=255)
       gggTGCTGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGGACGATGACCCGTTGCACCAGCGCTTGCGCTgg                                                                                                                                         <  8:80322/149‑1 (MQ=255)
       gggTGCTGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTgg                                                                                                                                         <  3:325261/149‑1 (MQ=255)
        ggTGCGGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGt                                                                                                                                        <  4:109966/149‑1 (MQ=255)
            cTGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCAAGCGCTGGTGCaa                                                                                                                                    <  5:26985/149‑2 (MQ=255)
             tGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCAcc                                                                                                                                   >  4:265657/1‑149 (MQ=255)
             tGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCAc                                                                                                                                    >  4:43696/1‑148 (MQ=255)
               gCCAAGACCCCGAATGGAGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCAAGCGCTGGTGCaac                                                                                                                                   <  4:295936/147‑3 (MQ=255)
                  aaGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGa                                                                                                                              <  8:204168/149‑1 (MQ=255)
                               gCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTAt                                                                                                                 <  7:114603/149‑1 (MQ=255)
                                cGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTAt                                                                                                                 >  5:193038/1‑148 (MQ=255)
                                 gATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGa                                                                                                               >  3:45858/1‑149 (MQ=255)
                                       tacGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGAt                                                                                                         >  8:30232/1‑149 (MQ=255)
                                       tacGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGAt                                                                                                         >  3:88769/1‑149 (MQ=255)
                                       tacGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGAt                                                                                                         >  4:15411/1‑149 (MQ=255)
                                                    gcagcGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACtg                                                                                            <  6:36466/149‑1 (MQ=255)
                                                               gACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCAcc                                                                                   <  3:43696/147‑1 (MQ=255)
                                                               gACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGa                                                                                 <  5:330315/149‑1 (MQ=255)
                                                                  cAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCt                                                                              >  1:39207/1‑149 (MQ=255)
                                                                      aaCCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTgcg                                                                           >  2:352155/1‑148 (MQ=255)
                                                                        ccaccCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCAt                                                                        <  5:210748/149‑1 (MQ=255)
                                                                         caccCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACtggtg                                                                                         >  6:97397/1‑131 (MQ=255)
                                                                         caccCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACtggtg                                                                                         <  5:97397/131‑1 (MQ=255)
                                                                          accCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATgg                                                                      >  4:341316/1‑149 (MQ=255)
                                                                               aaaTGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATc                                                                 <  3:15411/149‑1 (MQ=255)
                                                                                    cAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCgg                                                             <  7:30232/148‑1 (MQ=255)
                                                                                            gTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCTCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGTCGACGCAGTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTac                                                     <  4:45858/148‑1 (MQ=255)
                                                                                                  cAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTc                                              >  7:3799/1‑149 (MQ=255)
                                                                                                   aaCCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCg                                             <  3:280083/149‑1 (MQ=255)
                                                                                                     ccTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTc                                              >  1:146696/1‑146 (MQ=255)
                                                                                                       tGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGGTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCCCCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCtt                                          >  2:242827/1‑148 (MQ=255)
                                                                                                         gAACTGGCGCAAGCACTTAACTATGACCCGTTGCACCAGCGCTTTCGGGGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCaa                                       <  7:151260/149‑1 (MQ=255)
                                                                                                               gcgcAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGc                                 >  3:313623/1‑149 (MQ=255)
                                                                                                                 gcAAGCACTGAACGATGACCCGTTGCACCAGCGCTTgc                                                                                                                                              <  4:228906/38‑1 (MQ=255)
                                                                                                                  cAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCa                              <  4:88769/149‑1 (MQ=255)
                                                                                                                       aCTGAAAGAAGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATCACTGCCTGCGTTATGACCAGAATCCCGACGCAATAGTGGTCACCGACCTGAGAATCAGGATCCCGGCGAATAACACCATTCGCTTCAACAAGGCACAgcgc                          >  8:148856/1‑148 (MQ=255)
                                                                                                                         tGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTAACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGTGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGat                       >  1:253648/1‑149 (MQ=255)
                                                                                                                           aaCGATGACGCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTTGCTGCGTTATGACCAGATTGGCGACGCACTGTTGGTCACCGACATGCGCATGTGGATCCCGGGGAATTAGACCTTTCGCTTCAACATGGCTCAGCTCGat                       <  7:148856/147‑1 (MQ=255)
                                                                                                                             cGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGatcat                    >  6:152969/1‑148 (MQ=255)
                                                                                                                                  aCCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAggg              <  3:265657/149‑1 (MQ=255)
                                                                                                                                      gTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAg                <  6:267877/143‑1 (MQ=255)
                                                                                                                                      gTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAg                >  5:267877/1‑143 (MQ=255)
                                                                                                                                       gtgaaCCAGCGCAAGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTg         >  3:77996/5‑149 (MQ=255)
                                                                                                                                       gtgaaCCAGCGCAAGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTg         >  3:72915/5‑149 (MQ=255)
                                                                                                                                           aCCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATg     >  6:140076/1‑149 (MQ=255)
                                                                                                                                           aCCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATg     >  4:202918/1‑149 (MQ=255)
                                                                                                                                           aCCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATg     >  3:233524/1‑149 (MQ=255)
                                                                                                                                              aGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATGGTc  <  3:341316/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
GTCTGGCGGGTGCTGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATGGTC  >  NC_002947/3591429‑3591719

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