Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,010,677 A→T 7.4% L129Q (CTG→CAG)  PP_0871 ← glycine betaine/carnitine/choline ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,010,6770AT 7.4% 167.3 / 2.6 54L129Q (CTG→CAG) PP_0871glycine betaine/carnitine/choline ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (25/25);  new base T (2/2);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCCGCCGACGATGATCGGCACGGCGTTGGGCAGTTCCACTTGCCAGAGTTGCTGGCGCGGGGTCATGCCAATGCCGGTAGCGGCTTCCTTGAGCGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTC  >  NC_002947/1010533‑1010813
                                                                                                                                                |                                                                                                                                        
aCGCCGCCGACGATGATCGGCACGGCGTTGGGCAGTTCCACTTGCCAGAGTTGCTGGCGCGGGGTCATGCCAATGCCGGTAGCGGCTTCCTTGAGCGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCaggag                                                                                                                                      >  4:376393/1‑149 (MQ=255)
  gccgccGACGATGATCGGCACGGCGTTGGGCAGTTCCACTTGCCAGAGTTGCTGGCGCGGGGTCATGCCAATGCCGGTAGCGGCTTCCTTGAGCGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAggg                                                                                                                                    <  2:390308/149‑1 (MQ=255)
         acgaTGATCGGCACGGCGTTGGGCAGTTCCACTTGCCAGAGTTGCTGGCGCGGGGTCATGCCAATGCCGGTAGCGGCTTCCTTGAGCGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAg                                                                                                                             <  5:374011/149‑1 (MQ=255)
           gatgatCGGCACGGCGTTGGGCAGTTCCACTTGCCAGAGTTGCTGGCGCGGGGTCATGCCAATGCCGGTAGCGGCTTCCTTGAGCGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCTGGAGGGACGCGAgg                                                                                                                            >  5:122059/1‑148 (MQ=255)
           gatgatCGGCACGGCGTTGGGCAGTTCCACTTGCCAGAGTTGCTGGCGCGGGGTCATGCCAATGCCGGTAGCGGCTTCCTTGAGCGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCTGGAGGGACGCGAgg                                                                                                                            >  5:116249/1‑148 (MQ=255)
                tCGGCACGGCGTTGGGCAGTTCCACTTGCCAGAGTTGCTGGCGCGGGGTCATGCCAATGCCGGTAGCGGCTTCCTTGAGCGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCTGGAGGGACGCGAGGAACAGt                                                                                                                      <  6:122059/149‑1 (MQ=255)
                tCGGCACGGCGTTGGGCAGTTCCACTTGCCAGAGTTGCTGGCGCGGGGTCATGCCAATGCCGGTAGCGGCTTCCTTGAGCGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCTGGAGGGACGCGAGGAACAGt                                                                                                                      <  6:116249/149‑1 (MQ=255)
                    cACGGCGTTGGGCAGTTCCACTTGCCAGAGTTGCTGGCGCGGGGTCATGCCAATGCCGGTAGCGGCTTCCTTGAGCGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCg                                                                                                                   >  8:266917/1‑148 (MQ=255)
                     aCGGCGTTGGGCAGTTCCACTTGCCAGAGTTGCTGGCGCGGGGTCATGCCAATGCCGGTAGCGGCTTCCTTGAGCGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAgg                                                                                                                            >  8:26250/1‑138 (MQ=255)
                     aCGGCGTTGGGCAGTTCCACTTGCCAGAGTTGCTGGCGCGGGGTCATGCCAATGCCGGTAGCGGCTTCCTTGAGCGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAgg                                                                                                                            <  7:26250/138‑1 (MQ=255)
                               gcggTTCCACTTGCCAGAGTTGCTGGCGCGGGGTCATGCCAATGCCGGTAGCGGCTTCCTTGAGCGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGc                                                                                                       >  4:78783/4‑149 (MQ=255)
                                   ttCCACTTGCCAGAGTTGCTGGCGCGGGGTCATGCCAATGCCGGTAGCGGCTTCCTTGAGCGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCgg                                                                                                   >  8:131995/1‑149 (MQ=255)
                                        cTTGCCAGAGTTGCTGGCGCGGGGTCATGCCAATGCCGGTAGCGGCTTCCTTGAGCGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCg                                                                                              <  8:216588/149‑1 (MQ=255)
                                           gCCAGAGTTGCTGGCGCGGGGTCATGCCAATGCCGGTAGCGGCTTCCTTGAGCGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGa                                                                                             <  4:198212/147‑1 (MQ=255)
                                            ccAGAGTTGCTGGCGCGGGGTCATGCCAATGCCGGTAGCGGCTTCCTTGAGCGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGc                                                                                          >  1:381506/1‑149 (MQ=255)
                                             cAGAGTTGCTGGCGCGGGGTCATGCCAATGCCGGTAGCGGCTTCCTTGAGCGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGcc                                                                                         >  7:190501/1‑149 (MQ=255)
                                               gagTTGCTGGCGCGGGGTCATGCCAATGCCGGTAGCGGCTTCCTTGAGCGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGCCCa                                                                                       <  8:190501/149‑1 (MQ=255)
                                                agTTGCTGGCGCGGGGTCATGCCAATGCCGGTAGCGGCTTCCTTGAGCGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGCCCAg                                                                                      >  7:207423/1‑149 (MQ=255)
                                                 gTTGCTGGCGCGGGGTCATGCCAATGCCGGTAGCGGCTTCCTTGAGCGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGCCCAgg                                                                                     <  6:250216/149‑1 (MQ=255)
                                                              ggTCATGCCAATGCCGGTAGCGGCTTCCTTGAGCGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAAc                                                                                                                         >  8:227963/1‑100 (MQ=255)
                                                              ggTCATGCCAATGCCGGTAGCGGCTTCCTTGAGCGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAAc                                                                                                                         <  7:227963/100‑1 (MQ=255)
                                                                 cATGCCAATGCCGGTAGCGGCTTCCTTGAGCGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGCCCAGGATACTCAGGGCGATgg                                                                     >  4:178309/1‑149 (MQ=255)
                                                                  aTGCCATTGCCGGTAGCGGCTTCCTTGAGCGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCgg                                                                                                   <  1:164386/118‑1 (MQ=255)
                                                                       aaTGCCGGTAGCGGCTTCCTTGAGCGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGaa                                                               >  8:164236/1‑149 (MQ=255)
                                                                           ccGGTAGCGGCTTCCTTGAGCGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGc                                                           <  8:207423/149‑1 (MQ=255)
                                                                               tAGCGGCTTCCTTGAGCGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCa                                                         >  6:368905/1‑147 (MQ=255)
                                                                                   ggCTTCCTTGAGCGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAgggg                                                   >  5:89032/1‑149 (MQ=255)
                                                                                       tCCTTGAGCGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGAt                                               >  4:303690/1‑149 (MQ=255)
                                                                                         cTTGAGCGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATg                                              <  3:376393/148‑1 (MQ=255)
                                                                                          ttGAGCGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAgggg                                                   >  4:202986/1‑142 (MQ=255)
                                                                                          ttGAGCGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAgggg                                                   <  3:202986/142‑1 (MQ=255)
                                                                                          ttGAGCGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGgt                                            >  6:356849/1‑149 (MQ=255)
                                                                                            gAGCGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGgt                                            <  6:89032/147‑1 (MQ=255)
                                                                                              gCGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTg                                        >  7:247583/1‑149 (MQ=255)
                                                                                              gCGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTg                                        <  7:266917/149‑1 (MQ=255)
                                                                                               cGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGc                                       <  8:321940/149‑1 (MQ=255)
                                                                                                 aTGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCAAGAACGGCCAGAGGGGGGATGGTGTTGcc                                      >  6:112966/1‑148 (MQ=255)
                                                                                                       ggACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGtt                                >  6:218912/1‑148 (MQ=255)
                                                                                                       ggACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGtt                                >  2:254061/1‑148 (MQ=255)
                                                                                                       ggACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTg                               <  5:112966/149‑1 (MQ=255)
                                                                                                         aCGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTg                               >  2:108101/1‑147 (MQ=255)
                                                                                                         aCGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGACGGACGCGAGGAACAGTGCGAAGATCGCACGCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTg                               >  6:45351/1‑147 (MQ=255)
                                                                                                           gTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTgaag                            >  5:9642/1‑148 (MQ=255)
                                                                                                               ttCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACt                        <  3:178309/148‑1 (MQ=255)
                                                                                                                tCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTg                       <  7:131995/148‑1 (MQ=255)
                                                                                                                  aGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCAt                    <  3:78783/149‑1 (MQ=255)
                                                                                                                  aGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCAt                    <  8:247583/149‑1 (MQ=255)
                                                                                                                       ttCGTAGGTGTTGGGCAAGTTGGGGATGTGGGACGGGTTGAACAGTTCGTAGATCGCAGGCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTgaag                            <  2:86830/136‑1 (MQ=255)
                                                                                                                       ttCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGCCCAGGATAATCAGGGCTATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCAAACGTTgaag                            >  1:86830/1‑136 (MQ=255)
                                                                                                                             ggtggtGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCg          <  2:381506/148‑1 (MQ=255)
                                                                                                                                    cgcACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTc  <  5:217633/149‑1 (MQ=255)
                                                                                                                                     gcACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGtt                                <  5:123950/118‑1 (MQ=255)
                                                                                                                                     gcACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGtt                                >  6:123950/1‑118 (MQ=255)
                                                                                                                                     gcACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCt            <  7:351673/138‑1 (MQ=255)
                                                                                                                                     gcACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCt            >  8:351673/1‑138 (MQ=255)
                                                                                                                                                |                                                                                                                                        
ACGCCGCCGACGATGATCGGCACGGCGTTGGGCAGTTCCACTTGCCAGAGTTGCTGGCGCGGGGTCATGCCAATGCCGGTAGCGGCTTCCTTGAGCGATGCGGGGACGTTTTTCAGGCCTTCGTAGGTGTTGCGCACGATGGGCAGGAGGGACGCGAGGAACAGTGCGAAGATCGCAGGCCCGGCGCCGATGCCCAGGATACTCAGGGCGATGGCCAGAACGGCCAGAGGGGGGATGGTGTTGCCAACGTTGAAGAACTGCATGAAGCGCTCGGCTTTGTC  >  NC_002947/1010533‑1010813

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