Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,189,279 G→C 10.9% H457Q (CAC→CAG)  PP_2799 ← class III aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,189,2790GC10.9% 152.5 / 6.4 52H457Q (CAC→CAG) PP_2799class III aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (23/23);  new base C (3/2);  total (26/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 = 6.64e-01

TTGGGCCATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAG‑CGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCTTGTGA  >  NC_002947/3189133‑3189426
                                                                                                                                                  |                                                                                                                                                    
ttGGGCCATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTg                                                                                                                                                    >  2:188822/1‑149 (MQ=255)
   ggCCATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAGCCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCtg                                                                                                                                                 >  7:31598/1‑149 (MQ=255)
      cATGATGGTTACCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTCGCTCACAGCACATGCTGTGa                                                                                                                                              <  2:194076/149‑1 (MQ=255)
        tgatgGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAgc                                                                                                                                            <  1:188822/149‑1 (MQ=255)
          atgGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTg                                                                                                                                                    >  2:327163/1‑139 (MQ=255)
          atgGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTg                                                                                                                                                    <  1:327163/139‑1 (MQ=255)
          atgGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAgcgc                                                                                                                                          >  5:153686/1‑149 (MQ=255)
             gTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTCGCTCACAGCACCTGCTGTGAgcaaag                                                                                                                                        <  2:247362/148‑5 (MQ=255)
              ttCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAGCCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGcc                                                                                                                                       >  6:86292/1‑148 (MQ=255)
                 ccATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAAt                                                                                                                                   >  3:78851/1‑149 (MQ=255)
                        aCCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTcgc                                                                                                                            >  6:96802/1‑149 (MQ=255)
                          cGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGt                                                                                                                          >  2:274564/1‑149 (MQ=255)
                            aGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAg                                                                                                                                             <  1:328154/128‑1 (MQ=255)
                            aGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAg                                                                                                                                             >  2:328154/1‑128 (MQ=255)
                            aGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTca                                                                                                                        <  4:78851/149‑1 (MQ=255)
                                cGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTcagcca                                                                                                                    <  3:249637/149‑1 (MQ=255)
                                  aTCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTcagccag‑c                                                                                                                 <  5:346809/149‑1 (MQ=255)
                                             gTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAG‑CGTCCGCTCGaa                                                                                                      <  5:314864/149‑1 (MQ=255)
                                               cGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAG‑CGTCCGCTCGAAgg                                                                                                    >  3:238421/1‑149 (MQ=255)
                                                tctcTCGAAGCTTTCATCACTGTGTTCACGCACGTTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAG‑CGTCCGCTCGAAg                                                                                                     <  7:156027/144‑1 (MQ=255)
                                                    tCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATcc                                                                                                                                 >  3:205795/1‑116 (MQ=255)
                                                    tCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATcc                                                                                                                                 <  4:205795/116‑1 (MQ=255)
                                                    tCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCTTCGCTTCAGCCAG‑CTTCCGCTCGAAg                                                                                                     >  8:156027/1‑143 (MQ=255)
                                                      gAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAG‑CGTCCGCTCGAAGGCGCTg                                                                                               >  8:19874/1‑147 (MQ=255)
                                                                  aCTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAG‑CGTCCGCTCGAAGGCGCTGAAGATCGCCTCGat                                                                                 <  6:153686/149‑1 (MQ=255)
                                                                   cTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAG‑CGTCCGCTCGAAGGCGCTGAAGATCGCCTCGata                                                                                <  5:96802/149‑1 (MQ=255)
                                                                           acgcacgGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAG‑CGTCCGCTCGAAGCCGCTGAAGATCGCCTCGATATCGTCCTg                                                                        >  2:92877/1‑149 (MQ=255)
                                                                            cgcacgGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAG‑CGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTgg                                                                       <  7:19874/149‑1 (MQ=255)
                                                                               acgGTCAAGCCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAG‑CGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGggtg                                                                    <  5:86292/149‑1 (MQ=255)
                                                                                      aaCCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAG‑CGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATc                                                              >  3:211566/1‑148 (MQ=255)
                                                                                        cccTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAG‑CGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAAc                                                           <  4:211566/149‑1 (MQ=255)
                                                                                         ccTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAG‑CGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTgg                                                                       <  1:191462/136‑1 (MQ=255)
                                                                                         ccTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAG‑CGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTgg                                                                       >  2:191462/1‑136 (MQ=255)
                                                                                             aTGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTTCTGTGAGCGCGCCCAATCGGTCGCGTCAGCGAG‑CGTCCGCTCGAAGGCGCTGAAGATCGCCTCTATATCTTCCTGGGTGGTAATCAACgga                                                        <  1:92877/147‑1 (MQ=255)
                                                                                             aTGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCCGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACggag                                                       <  1:274564/149‑1 (MQ=255)
                                                                                              tGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAG‑CGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACggagga                                                     <  4:238421/149‑1 (MQ=255)
                                                                                                      gttTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAG‑CGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACggagga                                                     <  1:203908/141‑1 (MQ=255)
                                                                                                      gttTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAG‑CGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACggagga                                                     >  2:203908/1‑141 (MQ=255)
                                                                                                                     ccATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAG‑CGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGAt                              >  3:197208/1‑149 (MQ=255)
                                                                                                                      cATCAGATGCCGCCTTTTCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAG‑CGTCCGCTCGAAGGCGCTGAAGATAGCCTCGATATCGTCCTGGGTGGTAATAAACGGAGGAAAGAAGGCGATAGTGTCGTTGATg                             >  8:94080/1‑149 (MQ=255)
                                                                                                                       aTCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAG‑CGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATgg                            >  4:186269/1‑149 (MQ=255)
                                                                                                                         cAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAG‑CGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGACTGCGATCGTGTCGTAGATGgcg                          >  8:223471/1‑149 (MQ=255)
                                                                                                                          aGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAG‑CGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGgcg                          <  3:186269/148‑1 (MQ=255)
                                                                                                                                  ccTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAG‑CGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATc                 <  4:197208/149‑1 (MQ=255)
                                                                                                                                       gCTCACAGCACCTGCTGTGATCGCGCACAATCCGTCGAGTCAGCCCT‑CGTCCGCTCGAAGGCGCTGAAGATCGCCGCGATATCGTCCTGTGTGATCCTCAACGTCGGACAGAAGGCGCTCGGGGCGTAGAGGGCGCGAATGCTCAGcg             >  8:79303/1‑147 (MQ=255)
                                                                                                                                       gCTCACAGCACCTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAG‑CGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACg                                                          <  7:87047/103‑1 (MQ=255)
                                                                                                                                       gCTCACAGCACCTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAG‑CGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACg                                                          >  8:87047/1‑103 (MQ=255)
                                                                                                                                         tCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAG‑CGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACg                                                          <  4:239674/101‑1 (MQ=255)
                                                                                                                                         tCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAG‑CGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCAt          >  2:334233/1‑149 (MQ=255)
                                                                                                                                         tCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAG‑CGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCAt          >  3:43655/1‑149 (MQ=255)
                                                                                                                                         tCACAGCACCTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAG‑CGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCAt          >  7:78492/1‑149 (MQ=255)
                                                                                                                                               cagGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAG‑CGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCTtgt    <  1:334233/149‑1 (MQ=255)
                                                                                                                                                 ggTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAG‑CGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCTtgt    >  6:173745/1‑147 (MQ=255)
                                                                                                                                                 ggTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAG‑CGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCTTGTGa  <  6:94359/149‑1 (MQ=255)
                                                                                                                                                 ggTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAG‑CGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCTTGTGa  <  7:94080/149‑1 (MQ=255)
                                                                                                                                                 ggTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAG‑CGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCTTGTGa  <  4:43655/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TTGGGCCATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAG‑CGTCCGCTCGAAGGCGCTGAAGATCGCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCTTGTGA  >  NC_002947/3189133‑3189426

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