Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,731,849 C→T 9.4% A464A (GCG→GCA PP_3299 ← lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,731,8490CT 9.4% 145.7 / 3.4 53A464A (GCG→GCAPP_3299lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (25/23);  new base T (2/3);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.07e-01

TTACGCGGTATGAGCCCAATCCTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTC  >  NC_002947/3731706‑3731989
                                                                                                                                               |                                                                                                                                            
ttACGCGGTATGAGCCCAATCCTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCgctgc                                                                                                                                         >  5:165554/1‑149 (MQ=255)
     cgGTATTAGCCAAATCCTCAAAGAATGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTTTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCGGGTgctgctattg                                                                                                                                    <  6:391121/149‑5 (MQ=255)
          tGAGCCCAATCCTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTg                                                                                                                                <  3:382507/148‑1 (MQ=255)
          tGAGCCCAATCCTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTg                                                                                                                                <  3:376441/148‑1 (MQ=255)
              cccAATCCTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCaa                                                                                                                           >  6:125978/1‑149 (MQ=255)
              cccAATCCTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCaa                                                                                                                           <  5:7277/149‑1 (MQ=255)
                cAATCCTCAACGAAGGTGAGCACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACACAGGGAGAGGTAGGCAGAAATAGCAGCACGCGCTGCTTCAGTTTGCGCaaat                                                                                                                         >  7:13699/1‑147 (MQ=255)
                     cTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCGCGTgctgctat                                                                                                                                      <  1:215845/131‑3 (MQ=255)
                     cTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCGCGTgctgctat                                                                                                                                      >  2:215845/1‑129 (MQ=255)
                     cTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCg                                                                                                                    <  3:223885/149‑1 (MQ=255)
                      tCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCGCGTgctgctattg                                                                                                                                    <  1:207556/132‑5 (MQ=255)
                      tCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCGCGTgctgctattg                                                                                                                                    >  2:207556/1‑128 (MQ=255)
                       cAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCgg                                                                                                                   >  1:47331/1‑148 (MQ=255)
                             aGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCg                                                                                                            <  8:49132/149‑1 (MQ=255)
                                  aGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTc                                                                                                       <  6:121514/149‑1 (MQ=255)
                                          ttctGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAg                                                                                               >  5:271977/1‑149 (MQ=255)
                                                aCCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAgctgct                                                                                          >  7:119325/1‑148 (MQ=255)
                                                aCCAAGGAATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACAACCCGCCACCCGGGTCTAGGTAGGCAGCAAAAGCAGCACGCGCGCCTTCATTTTGCGAAACAGCTCGCAAATCCGATGTCGTGAGAAgctgat                                                                                          >  6:156130/1‑146 (MQ=255)
                                                 ccAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGa                                                                                        >  1:391644/1‑149 (MQ=255)
                                                    aGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCt                                                                                     >  2:123583/1‑149 (MQ=255)
                                                         tttACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGAACTGCGtt                                                                                <  3:235510/149‑1 (MQ=255)
                                                                                tgCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCACTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCg                                                                                  >  5:124622/1‑124 (MQ=255)
                                                                                tgCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCACTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCg                                                                                  <  6:124622/124‑1 (MQ=255)
                                                                                  ccACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTaca                                                       >  5:210923/1‑149 (MQ=255)
                                                                                    acaTCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACgcgc                                                                                                                                            >  3:236638/1‑62 (MQ=255)
                                                                                    acaTCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACgcgc                                                                                                                                            <  4:236638/62‑1 (MQ=255)
                                                                                    acaTCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACAcc                                                     >  8:195638/1‑149 (MQ=255)
                                                                                           aGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTAGACCGTTCTg                                               <  3:154782/148‑1 (MQ=255)
                                                                                           aGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTgg                                              <  8:119325/149‑1 (MQ=255)
                                                                                           aGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTg                                               <  6:165554/148‑1 (MQ=255)
                                                                                           aGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTg                                               <  2:47331/148‑1 (MQ=255)
                                                                                               aGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATg                                          >  8:140764/1‑149 (MQ=255)
                                                                                               aGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATg                                          >  3:241346/1‑149 (MQ=255)
                                                                                                 ttGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCt                                        <  5:125978/149‑1 (MQ=255)
                                                                                                   gACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTAc                                      >  8:119927/1‑149 (MQ=255)
                                                                                                      caccCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACg                                     >  6:202718/1‑147 (MQ=255)
                                                                                                       accCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAg                                  >  8:355756/1‑149 (MQ=255)
                                                                                                         ccGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGa                                <  8:13699/149‑1 (MQ=255)
                                                                                                         ccGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGa                                >  1:49247/1‑149 (MQ=255)
                                                                                                          cGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGAc                               >  3:262613/1‑149 (MQ=255)
                                                                                                           gCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGAc                               >  8:22467/1‑148 (MQ=255)
                                                                                                            ccACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACg                              >  2:16171/1‑148 (MQ=255)
                                                                                                                ccAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAAcgt                         >  3:191520/1‑149 (MQ=255)
                                                                                                                 cAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAAcgtc                        <  2:391644/149‑1 (MQ=255)
                                                                                                                      cGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTgct                   <  1:16171/149‑1 (MQ=255)
                                                                                                                       gAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTgctc                  <  4:241346/149‑1 (MQ=255)
                                                                                                                                agcaATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAgg          >  2:351967/1‑148 (MQ=255)
                                                                                                                                 gcaATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGCATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGa        >  4:48228/1‑149 (MQ=255)
                                                                                                                                 gcaATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGa        <  5:202718/149‑1 (MQ=255)
                                                                                                                                   aaTAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATa      >  4:123043/1‑149 (MQ=255)
                                                                                                                                    aTAGCAGCACGCGCTTCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCTCTCCAGGTATAc     <  7:296682/149‑1 (MQ=255)
                                                                                                                                    aTAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATAc     <  7:195638/149‑1 (MQ=255)
                                                                                                                                    aTAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATAc     <  4:191520/149‑1 (MQ=255)
                                                                                                                                    aTAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATAc     <  1:123583/149‑1 (MQ=255)
                                                                                                                                       gcagcaCGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTc  <  7:140764/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                            
TTACGCGGTATGAGCCCAATCCTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTC  >  NC_002947/3731706‑3731989

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: