Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,597,751 C→A 9.4% Q686H (CAG→CAT) ‡ PP_4071 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,597,7510CA 9.4% 162.8 / 6.9 53Q686H (CAG→CAT) ‡PP_4071membrane protein
Reads supporting (aligned to +/- strand):  ref base C (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 = 8.51e-01

TTGAAGCCTTTGGCAGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGTCTCGTTGCAG  >  NC_002947/4597605‑4597897
                                                                                                                                                  |                                                                                                                                                  
ttGAAGCCTTTGGCAGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCctgctg                                                                                                                                                  >  7:39776/1‑149 (MQ=255)
      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                                                                                                                                             >  7:330597/1‑134 (MQ=255)
                    cATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCa                                                                                                                                             <  8:330597/134‑1 (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                                                                                                                                            >  6:147954/1‑119 (MQ=255)
                                    acacCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAg                                                                                                                                            <  5:147954/119‑1 (MQ=255)
                                      acCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCt                                                                                                            <  3:74829/149‑1 (MQ=255)
                                                cccTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGt                                                                                                  >  8:249271/1‑149 (MQ=255)
                                                    ggCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTtg                                                                                                                          <  5:210473/121‑1 (MQ=255)
                                                    ggCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTtg                                                                                                                          >  6:210473/1‑121 (MQ=255)
                                                     gCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTctgc                                                                                              <  2:30927/148‑1 (MQ=255)
                                                             cGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTgg                                                                                     <  7:249271/149‑1 (MQ=255)
                                                                          gTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGc                                                                        <  7:195939/149‑1 (MQ=255)
                                                                               cGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTa                                                                                                          >  7:55634/1‑110 (MQ=255)
                                                                               cGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTa                                                                                                          <  8:55634/110‑1 (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              >  7:192951/1‑149 (MQ=255)
                                                                                                                                    tGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAg              >  8:31218/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)
                                                                                                                                                  |                                                                                                                                                  
TTGAAGCCTTTGGCAGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGAGCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGTCTCGTTGCAG  >  NC_002947/4597605‑4597897

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