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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,750,7380TA 8.5% 134.6 / 3.2 47K88* (AAG→TAG) ‡PP_3315membrane protein
Reads supporting (aligned to +/- strand):  ref base T (19/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 = 7.76e-01

GCGCCAGCGAGCGCCGCGACTTTGGCCGGGCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAGCTGTACCAGCGACACCACGATAACCA  >  NC_002947/3750590‑3750878
                                                                                                                                                    |                                                                                                                                            
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                                                                                                                                  >  8:337198/1‑125 (MQ=255)
                                   aTCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATTCAAGCATTCCATcgt                                                                                                                                  <  7:337198/126‑2 (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)
                                                                                                                                                    |                                                                                                                                            
GCGCCAGCGAGCGCCGCGACTTTGGCCGGGCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAGCTGTACCAGCGACACCACGATAACCA  >  NC_002947/3750590‑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: ‑