Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,597,755 T→A 9.0% K685M (AAG→ATG) ‡ PP_4071 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,597,7550TA 9.0% 150.8 / 3.5 53K685M (AAG→ATG) ‡PP_4071membrane protein
Reads supporting (aligned to +/- strand):  ref base T (20/28);  new base A (3/2);  total (23/30)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTTTGGCAGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGTCTCGTTGCAGCTGCCC  >  NC_002947/4597611‑4597903
                                                                                                                                                |                                                                                                                                                    
ccTTTGGCAGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAg                                                                                                                                                  <  6:22510/149‑1 (MQ=255)
  tttGGCAGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGc                                                                                                                                                >  8:116639/1‑149 (MQ=255)
    tGGCAGTGACCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGcc                                                                                                                                               <  8:285149/148‑1 (MQ=255)
     ggCAGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGccc                                                                                                                                              <  3:295688/148‑1 (MQ=255)
      gCAGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCg                                                                                                                                             <  8:203365/148‑1 (MQ=255)
      gCAGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGc                                                                                                                                            <  8:151701/149‑1 (MQ=255)
        aGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGccc                                                                                                                                              <  8:61395/145‑1 (MQ=255)
        aGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGccc                                                                                                                                              >  7:61395/1‑145 (MQ=255)
         gTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGAAGCAGCAGGcccagga                                                                                                                                          <  4:125972/148‑5 (MQ=255)
         gTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGAAGCAGCAGGcccag                                                                                                                                            >  3:125972/1‑144 (MQ=255)
            aGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGAAGCAGCAGGcc                                                                                                                                               <  6:17129/140‑1 (MQ=255)
            aGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGAAGCAGCAGGcc                                                                                                                                               >  5:17129/1‑140 (MQ=255)
              cATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCa                                                                                                                                                   <  8:330597/134‑1 (MQ=255)
              cATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCa                                                                                                                                                   >  7:330597/1‑134 (MQ=255)
                   gcgcCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTtgt                                                                                                                               <  7:21415/149‑1 (MQ=255)
                    cgcCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGGTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTtgtg                                                                                                                              <  5:332727/149‑1 (MQ=255)
                    cgcCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCCGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATCCAACTCCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTtgtg                                                                                                                              >  8:195939/1‑149 (MQ=255)
                     gcCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACCCGTTGTgg                                                                                                                             >  3:134023/1‑149 (MQ=255)
                          cAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACg                                                                                                                        >  4:337281/1‑149 (MQ=255)
                           aaaaCACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACgg                                                                                                                       <  8:39776/149‑1 (MQ=255)
                              acacCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAg                                                                                                                                                  <  5:147954/119‑1 (MQ=255)
                              acacCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAg                                                                                                                                                  >  6:147954/1‑119 (MQ=255)
                                acCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCt                                                                                                                  <  3:74829/149‑1 (MQ=255)
                                          cccTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGt                                                                                                        >  8:249271/1‑149 (MQ=255)
                                              ggCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTtg                                                                                                                                >  6:210473/1‑121 (MQ=255)
                                              ggCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTtg                                                                                                                                <  5:210473/121‑1 (MQ=255)
                                               gCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTctgc                                                                                                    <  2:30927/148‑1 (MQ=255)
                                                       cGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTgg                                                                                           <  7:249271/149‑1 (MQ=255)
                                                                    gTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGc                                                                              <  7:195939/149‑1 (MQ=255)
                                                                         cGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTa                                                                                                                <  8:55634/110‑1 (MQ=255)
                                                                         cGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTa                                                                                                                >  7:55634/1‑110 (MQ=255)
                                                                             aTCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCg                                                                     >  6:8273/1‑149 (MQ=255)
                                                                             aTCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTAAGCAATGTCTGCTGCACGTGGCTGACACTCTGGCGATACTGCg                                                                     >  8:144298/1‑149 (MQ=255)
                                                                              tCGAACGCAGGCCGCACCTGATGACGCAGGGCCAGACCGGGATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCTCAAACACGAGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGc                                                                      <  4:273312/147‑1 (MQ=255)
                                                                              tCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCAGGGCCTGCTGCTTCAGGCCCGCATACACGCCTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCg                                                                     >  3:273312/1‑148 (MQ=255)
                                                                               cGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCACCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTATGCTGCACGTGGCTGACACGCTGGCGATACTGCgg                                                                    >  2:19495/1‑148 (MQ=255)
                                                                                  aCGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGCTACTGCGGCAgg                                                                >  7:289364/1‑149 (MQ=255)
                                                                                          cccACCTGATGACGCAGGTCGAGACCGTTATGCAACTGCTCCTGGTCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCg                                                                     <  3:258565/136‑1 (MQ=255)
                                                                                               cTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGt                                                   <  3:337281/149‑1 (MQ=255)
                                                                                                  atgaCGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATTTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAg                                                 <  4:134023/148‑1 (MQ=255)
                                                                                                     aCGCAGGTCCAGACCGGTGTGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTctgc                                                                                                    >  7:48878/1‑94 (MQ=255)
                                                                                                     aCGCAGGTCCAGACCGGTGTGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTctgc                                                                                                    <  8:48878/94‑1 (MQ=255)
                                                                                                       gCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAgc                                           <  7:116639/149‑1 (MQ=255)
                                                                                                        cAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAgcg                                          <  6:304789/149‑1 (MQ=255)
                                                                                                                 aCCGGTATGCAACTGCTCCTGGGCCTGCTGCTGCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACtt                                   >  5:88466/1‑147 (MQ=255)
                                                                                                                  ccGGTTTGCTACTTCTCCTTGGCCGGCTTCTTCATGCCCGCATACACCCGTTGTGGCAACGGCACCTTACGCAATTTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATATGCGTAGGTCAGCGCATACTTCa                                 <  7:55097/148‑1 (MQ=255)
                                                                                                                  ccGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAg                                <  5:8273/149‑1 (MQ=255)
                                                                                                                              tGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAg                    >  8:31218/1‑149 (MQ=255)
                                                                                                                              tGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAg                    >  7:192951/1‑149 (MQ=255)
                                                                                                                                cTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAg                    <  7:144298/147‑1 (MQ=255)
                                                                                                                                   cTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGtctc               >  3:255013/1‑149 (MQ=255)
                                                                                                                                     gggCCTGAAGCAGCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGTCTCGt             >  2:100303/1‑149 (MQ=255)
                                                                                                                                      ggCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGAACCTTACGGATTGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGTCTCg              <  1:19495/147‑1 (MQ=255)
                                                                                                                                            ctgctTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGTCTCGTTGCAg        <  3:131024/147‑1 (MQ=255)
                                                                                                                                                ttCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAg                    <  3:119776/131‑1 (MQ=255)
                                                                                                                                                ttCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAg                    >  4:119776/1‑131 (MQ=255)
                                                                                                                                                ttCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGTCTCGTTGCAGCTGccc  >  2:11725/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CCTTTGGCAGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGTCTCGTTGCAGCTGCCC  >  NC_002947/4597611‑4597903

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