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.G92.9% 129.0 / 4.3 43intergenic (+46/+96)aroP‑I/PP_0928aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base G (17/22);  total (18/25)
Fisher's exact test for biased strand distribution p-value = 2.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.70e-01
*NC_0029471,070,2462.A93.0% 123.2 / 4.7 43intergenic (+46/+96)aroP‑I/PP_0928aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base A (18/22);  total (18/25)
Fisher's exact test for biased strand distribution p-value = 2.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.85e-01

GCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GT‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTAC  >  NC_002947/1070104‑1070378
                                                                                                                                                 ||                                                                                                                                    
gcggccgATTGGCATCTCTGTAATGCTTATCCCGGCATGGATAGCTGTGCTGTGGGTGGCTTTCAAGCTGAATAAGGCTCGCCAGGGCAAGGGGTGATCGGGTTTGTTCTTAAGGGGTATGCCTTTAGGTTTTCGGCGCCT‑‑GT‑‑Gagat                                                                                                                                 <  1:300644/145‑4 (MQ=255)
  cgccgATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GT‑‑Gagatcg                                                                                                                               <  2:255201/148‑6 (MQ=255)
            cATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCcgcgc                                                                                                                       <  1:70984/148‑1 (MQ=255)
             aTCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCcgcgcg                                                                                                                      >  3:105390/1‑148 (MQ=255)
              tctcGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCgcgg                                                                                                                     >  7:7140/1‑148 (MQ=255)
                  ggTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCgagggcgg                                                                                                                 >  8:130717/1‑148 (MQ=255)
                         cTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                          <  7:15886/148‑1 (MQ=255)
                                 tgtaTGGCTAGCTGTGCGGTGTGTTGCTTTCTAGCTGAAGATGGCTCGCATGGCCAAGGGCTGATCGGGTTAGTTGGGTATGGGTATTCGGTTAGGTTTTGGGCGCCTGGGA‑‑GATCGAGCGCGGCGCGGGCGGCGCTCGTTCTcaca                                                                                                   <  6:245746/144‑1 (MQ=255)
                                 ggCATGGATCGCTGTGCTGGGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAg                                                                                                  <  1:220624/148‑1 (MQ=255)
                                       gATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGAt                                                                                             >  8:181413/1‑147 (MQ=255)
                                       gATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCCCGGGCGGCGCGCCATCTCACAcgcgct                                                                                             >  1:162178/1‑145 (MQ=38)
                                            ctgtgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTCTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAg                                                                                        <  5:127065/147‑1 (MQ=255)
                                              gtgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGc                                                                                       >  6:298526/1‑146 (MQ=255)
                                                 ctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTaaa                                                                                  <  3:304119/148‑1 (MQ=255)
                                                  tgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTaa                                                                                   <  6:139336/146‑1 (MQ=255)
                                                      ggtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCa                                                                              <  7:258754/147‑1 (MQ=255)
                                                      gcgggCGTTGAAGCTGATGAAGGCTTGCGAGGCCAAGGGCTGTTCGGGTGTGTCCGGAAGGCGTATGCCGTTAGGTTTTCGGCGCAT‑‑GTGAGAGCGAGCGTCGGGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCa                                                                              <  8:248103/144‑1 (MQ=21)
                                                       gtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCa                                                                              <  5:146447/146‑1 (MQ=255)
                                                       gtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATg                                                                            >  4:180991/1‑148 (MQ=255)
                                                              tCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTg                                                                      <  3:180991/147‑1 (MQ=255)
                                                                  gCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTggcggc                                                                 >  4:69199/1‑148 (MQ=255)
                                                                      aagaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGccc                                                               >  5:312706/1‑146 (MQ=255)
                                                                       agaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAATGGTTAAGCCGTTAGGTTTTCCGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTg                                                             <  5:92446/147‑1 (MQ=37)
                                                                       agaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCg                                                                          <  8:1382/134‑1 (MQ=255)
                                                                       agaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCg                                                                          >  7:1382/1‑134 (MQ=255)
                                                                       agaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTg                                                             >  5:322865/1‑147 (MQ=255)
                                                                         aagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCGCAGGCGATGAAAGCCTAAAGGCATGCGTCTGGAGGCCCTGAcg                                                          >  7:137653/1‑148 (MQ=38)
                                                                             cTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGtt                                                      <  1:138758/148‑1 (MQ=255)
                                                                               cGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTATCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTcc                                                    >  2:263204/1‑148 (MQ=255)
                                                                                 ccAGGCCAAGGGCTGATCGGGGTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCGCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTcc                                                    <  5:325779/146‑1 (MQ=38)
                                                                                   aGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCtt                                                <  2:80303/148‑1 (MQ=255)
                                                                                   aGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGGTCCCCtt                                                >  4:49770/1‑148 (MQ=255)
                                                                                       cAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCAt                                            <  8:98613/148‑1 (MQ=255)
                                                                                        aaGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATg                                           <  4:105390/148‑1 (MQ=255)
                                                                                             cTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGccc                                        >  6:179561/1‑146 (MQ=255)
                                                                                                 tCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAt                                  <  8:7140/148‑1 (MQ=255)
                                                                                                                aGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTTAGATCGAGCGCAGCGCGGGCGGCGCTAGATCTCACCGGCGCTGCAAGCCTAAAGGCATGCGTCTTGCGGCCCGGCCGAGGTCCCCGTCCAGGCCCCGAAATGATTCCTCCGGCaaa                   >  4:179926/1‑147 (MQ=255)
                                                                                                                     tATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgc              >  1:96625/1‑148 (MQ=255)
                                                                                                                            ttAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCCATGCTGCCTCCGCCAAGCACGCACAgtgg       >  8:236123/1‑148 (MQ=255)
                                                                                                                            ttAGGTTTTCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtgg       <  2:162178/148‑1 (MQ=255)
                                                                                                                             tAGGTTTTCGGCGCCT‑‑GTGAGATGGAGCGGCGCGCGGGCGGCGCTCGATCTCACATGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTGCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtgg       <  5:298526/147‑1 (MQ=38)
                                                                                                                                 ttttCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTAc  <  7:181413/148‑1 (MQ=255)
                                                                                                                                  tttCGGCGCCT‑‑GTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTAc  <  5:291733/147‑1 (MQ=255)
                                                                                                                                                 ||                                                                                                                                    
GCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GT‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTAC  >  NC_002947/1070104‑1070378

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: