Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,070,248 +GA intergenic (+48/+94) aroP‑I → / ← PP_0928 aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,070,2461.G91.3% 139.1 / 7.4 47intergenic (+46/+96)aroP‑I/PP_0928aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein
Reads supporting (aligned to +/- strand):  ref base . (3/1);  new base G (17/25);  total (21/26)
Fisher's exact test for biased strand distribution p-value = 3.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0029471,070,2462.A91.3% 137.2 / 7.0 47intergenic (+46/+96)aroP‑I/PP_0928aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein
Reads supporting (aligned to +/- strand):  ref base . (3/1);  new base A (18/24);  total (21/26)
Fisher's exact test for biased strand distribution p-value = 3.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGT‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAA‑GCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAA  >  NC_002947/1070098‑1070383
                                                                                                                                                       ||                                                                                                                                          
ccTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑tgt                                                                                                                                              <  7:144698/149‑1 (MQ=255)
  tGACACCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑tgt                                                                                                                                              >  6:311070/1‑147 (MQ=255)
    aCGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGT‑‑Ga                                                                                                                                          <  1:291337/147‑1 (MQ=255)
       ccgccgATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGT‑‑Gagatc                                                                                                                                      >  5:291040/1‑144 (MQ=255)
         gccgATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGT‑‑Gagatcgag                                                                                                                                   >  4:56977/1‑142 (MQ=255)
                  cATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTACGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCcgcgcg                                                                                                                            <  3:62310/149‑1 (MQ=255)
                   aTCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCgcgg                                                                                                                           >  8:258822/1‑149 (MQ=255)
                      tcGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCgcgggcg                                                                                                                        >  8:32649/1‑149 (MQ=255)
                            aTGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                  <  3:316263/149‑1 (MQ=255)
                                 gATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                             >  2:235984/1‑149 (MQ=255)
                                  aTCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGCGTATGCCGTTAGGTTATCATCCCC‑‑TGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGAtctc                                                                                                            <  1:50479/149‑1 (MQ=37)
                                         cATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGACGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCg                                                                                                     <  5:311070/149‑1 (MQ=255)
                                                  ctgtgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAA‑Gcc                                                                                           >  8:19623/1‑149 (MQ=255)
                                                     tgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCCAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAA‑GCCTaa                                                                                        <  7:258822/149‑1 (MQ=255)
                                                       ctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAA‑GCCTAAAg                                                                                      <  7:32649/149‑1 (MQ=255)
                                                           gggtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAA‑GCCTAAAGGCAt                                                                                  <  7:19623/149‑1 (MQ=255)
                                                           gggtggCTTTCAAGCGGAAGAAGGCTCGCCAGGCCAAGGGCTGAGCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAA‑GCCTAAAGGCAt                                                                                  <  5:133180/149‑1 (MQ=38)
                                                              tggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAA‑GCCTAAAGGCATGCg                                                                               >  3:260490/1‑149 (MQ=255)
                                                                gCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAA‑GCCTAAAGGCATGCGTc                                                                             >  5:38286/1‑149 (MQ=255)
                                                                  tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGT‑‑GATAGAGCGCCGCGCGGGGGGCGCGCGTGGTGAAAGGCGAGGAAA‑GCCGAAAGGCATGCGGCTg                                                                           >  3:136/1‑149 (MQ=255)
                                                                  tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAA‑GCCTAAAGGCATGCGTCTg                                                                           >  3:197829/1‑149 (MQ=255)
                                                                   ttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAA‑GCCTAAAGGCATGCGTCAgg                                                                          >  2:63429/1‑149 (MQ=255)
                                                                     cAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGATCGTACGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAA‑GCCTAAAGGCATGCGTCTggc                                                                         <  3:239967/148‑1 (MQ=38)
                                                                        gCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAA‑GCCTAAAGGCATGCGt                                                                              <  4:159906/140‑1 (MQ=255)
                                                                        gCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAA‑GCCTAAAGGCATGCGt                                                                              >  3:159906/1‑140 (MQ=255)
                                                                          tGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAA‑GCCTAAAGGCATGCGTCTGGCGGCCCt                                                                   >  5:306253/1‑149 (MQ=255)
                                                                             agaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAA‑GCCTAAAGGCATGCGTCTGGCGGCCCTGAc                                                                <  4:260490/149‑1 (MQ=255)
                                                                              gaagGCTAGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAACGCCTAAAGGCATGCGTCTGGCGGCCCTGAc                                                                <  1:195707/149‑1 (MQ=38)
                                                                                   cTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAA‑GCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTc                                                          <  1:327476/149‑1 (MQ=255)
                                                                                              aaGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAA‑GCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGc                                               <  4:197829/149‑1 (MQ=255)
                                                                                                 ggCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAA‑GCCTAAAg                                                                                      <  6:174251/107‑1 (MQ=255)
                                                                                                 ggCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAA‑GCCTAAAg                                                                                      >  5:174251/1‑107 (MQ=255)
                                                                                                 ggCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAA‑GCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGcccc                                            <  3:56977/149‑1 (MQ=255)
                                                                                                   cTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAA‑GCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCt                                           >  8:245264/1‑148 (MQ=255)
                                                                                                    tGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTTAGATAGAGCCCCCCGCGGGCGGCCCTCCATCTCACAAGCGATGAAA‑CCCAAAAAGCATTCGTATTGCGGCCCTGACGCGTTCCCCTTACATGCCCCta                                          >  6:210236/1‑147 (MQ=255)
                                                                                                     gATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAA‑GCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCa                                         >  6:1323/1‑148 (MQ=255)
                                                                                                       tCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAA‑GCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtg                                      <  7:245264/149‑1 (MQ=255)
                                                                                                         gggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAA‑GCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgct                                    <  6:291040/149‑1 (MQ=255)
                                                                                                            tttGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAA‑GCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGcc                                 >  6:318583/1‑149 (MQ=255)
                                                                                                               gTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCTGGCGGCGCTCGATCTCACAGGCGATGAAA‑GCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTcc                              >  8:241143/1‑149 (MQ=255)
                                                                                                                    tAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAA‑GCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGcc                           >  1:285550/1‑147 (MQ=255)
                                                                                                                           tATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAA‑GCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgca                  <  6:38286/149‑1 (MQ=255)
                                                                                                                           tATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAA‑GCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgca                  >  2:18427/1‑149 (MQ=255)
                                                                                                                                cGTTAGGTTTTCGGCGCC‑‑TGTGCGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAA‑GCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtg             <  4:136/149‑1 (MQ=255)
                                                                                                                                cGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAA‑GCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtg             <  1:18427/149‑1 (MQ=255)
                                                                                                                                  ttAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAA‑GCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtggt           <  5:210236/149‑1 (MQ=255)
                                                                                                                                    aGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAA‑GCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtggtg          <  2:94321/148‑1 (MQ=255)
                                                                                                                                       ttttCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAA‑GCCTAAAGGCATGCGTCTGGCGTCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACa      <  5:1323/149‑1 (MQ=255)
                                                                                                                                            ggCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAA‑GCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCaa  <  5:318583/148‑1 (MQ=255)
                                                                                                                                                       ||                                                                                                                                          
CCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGT‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAA‑GCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAA  >  NC_002947/1070098‑1070383

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