Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,750,737 T→A 8.5% K88M (AAG→ATG) ‡ PP_3315 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,750,7370TA 8.5% 133.4 / 2.7 47K88M (AAG→ATG) ‡PP_3315membrane protein
Reads supporting (aligned to +/- strand):  ref base T (18/24);  new base A (2/2);  total (21/26)
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.36e-01

GGCGCCAGCGAGCGCCGCGACTTTGGCCGGGCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAGCTGTACCAGCGACACCACGATAACCA  >  NC_002947/3750589‑3750878
                                                                                                                                                    |                                                                                                                                             
ggCGCCAGCGAGCGCCGCGACTTTGGCCGGGCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCt                                                                                                                                               <  4:177593/149‑1 (MQ=255)
 gcgcCAGCGAGCGCCGCGACTTTGGCCGGGCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCtt                                                                                                                                              <  3:19153/149‑1 (MQ=255)
   gcCAGCGAGCGCCGCGACTTTGGCCGGGCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGa                                                                                                                                            <  5:29352/149‑1 (MQ=255)
               cgcgACTTTGGCCGGGCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAAt                                                                                                                                          <  4:34539/139‑1 (MQ=255)
               cgcgACTTTGGCCGGGCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAAt                                                                                                                                          >  3:34539/1‑139 (MQ=255)
                     tttGGCCGGGCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATTCAAGCATTCCATc                                                                                                                                    >  1:346456/1‑139 (MQ=255)
                     tttGGCCGGGCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATTCAAGCATTCCATc                                                                                                                                    <  2:346456/139‑1 (MQ=255)
                       tGGCCGGGCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATcagca                                                                                                                          <  7:8359/147‑1 (MQ=255)
                        ggCCGGGCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGAGCCACCCGATCAAGAGGACCGGCTGGCCCTGAATGAGGCGGCCATGGCGCACGCCAACGCGCCCGATGGACTGCGTGACTGCGATCGCGGTgagcctcg                                                                                                                       >  4:293658/1‑145 (MQ=255)
                                    aTCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATTCAAGCATTCCATcgt                                                                                                                                  <  7:337198/126‑2 (MQ=255)
                                    aTCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATTCAAGCATTCCATcgt                                                                                                                                  >  8:337198/1‑125 (MQ=255)
                                    aTCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATAc                                                                                                           <  3:55632/149‑1 (MQ=255)
                                    aTCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATAc                                                                                                           >  7:60698/1‑149 (MQ=255)
                                        ggATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGaa                                                                                                       <  8:321371/149‑1 (MQ=255)
                                                  tCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTgg                                                                                             <  6:338556/149‑1 (MQ=255)
                                                    aagATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGTCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGaa                                                                                           <  8:18043/147‑1 (MQ=255)
                                                    aCGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGaa                                                                                           <  6:67301/149‑1 (MQ=255)
                                                    aCGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGa                                                                                            <  3:357723/148‑1 (MQ=255)
                                                    aCGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGa                                                                                            <  7:242021/148‑1 (MQ=255)
                                                     cGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGCCCATACCGAACACATTCTgg                                                                                             >  8:317592/1‑146 (MQ=255)
                                                      gATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATTAGCAGCGTCATGACCATACCGAACACATTCTgg                                                                                             >  6:272548/1‑145 (MQ=255)
                                                                        gATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCt                                                                                               >  7:284077/1‑125 (MQ=255)
                                                                            aaCCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGCTGGACCACCGGAtcttc                                                                   <  4:122063/149‑5 (MQ=255)
                                                                                     aaGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCgat                                                          <  6:201337/149‑1 (MQ=255)
                                                                                       gagGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCgatga                                                        >  3:278224/1‑149 (MQ=255)
                                                                                         ggACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATTAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATc                                                      <  5:272548/149‑1 (MQ=255)
                                                                                           aCCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATcag                                                    <  7:198804/149‑1 (MQ=255)
                                                                                             cGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATcagca                                                  <  8:60698/149‑1 (MQ=255)
                                                                                                 ttGACCTGAATGATACTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCt                                              >  6:127319/1‑149 (MQ=255)
                                                                                                   gACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTg                                             >  6:272537/1‑148 (MQ=255)
                                                                                                    aCCTGAATGATGCTGCCATGGCGCAAGACAACGCGCAGGATGGAATGCTTGAATTCGATCGCGATCAGCAGCGTCATGACGATACCGAACACTTTCTGGAAGATGCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCATGCTGaa                                           <  1:386373/149‑1 (MQ=255)
                                                                                                    aCCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTg                                             >  5:335616/1‑147 (MQ=255)
                                                                                                     ccTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGaa                                           <  7:45996/148‑1 (MQ=255)
                                                                                                      cTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGa                                         >  1:83602/1‑149 (MQ=255)
                                                                                                         aaTGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCg                                                            <  7:102429/127‑1 (MQ=255)
                                                                                                         aaTGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCg                                                            >  8:102429/1‑127 (MQ=255)
                                                                                                                ctgcCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCTGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGAt                               <  1:271688/149‑1 (MQ=255)
                                                                                                                  gcCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGAt                               >  6:264232/1‑147 (MQ=255)
                                                                                                                  gcCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCAGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGAt                               >  4:375534/1‑147 (MQ=255)
                                                                                                                     aTGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATc                                    >  3:273696/1‑139 (MQ=255)
                                                                                                                     aTGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATc                                    <  4:273696/139‑1 (MQ=255)
                                                                                                                          gcAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAGCTGTAc                      >  5:141398/1‑148 (MQ=255)
                                                                                                                           cAAGACAACGCGCCCGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAGCTGTACCa                    <  3:293658/149‑1 (MQ=255)
                                                                                                                             aGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAGCTGTACCAGc                  <  5:272537/149‑1 (MQ=255)
                                                                                                                                      gcACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATcagc                                                   >  6:267099/1‑107 (MQ=255)
                                                                                                                                      gcACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATcagc                                                   <  5:267099/107‑1 (MQ=255)
                                                                                                                                            tGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAGCTGTACCAGCGACACCACGATAacc   >  6:272482/1‑149 (MQ=255)
                                                                                                                                                aTGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAGCTGTACCAGCGACACCACGATAacca  >  5:197098/1‑146 (MQ=255)
                                                                                                                                                aTGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAGCTGTACCAGCGACACCACGATAacca  <  6:197098/146‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
GGCGCCAGCGAGCGCCGCGACTTTGGCCGGGCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAGCTGTACCAGCGACACCACGATAACCA  >  NC_002947/3750589‑3750878

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