Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,054,778 T→A 14.6% intergenic (‑253/+55) PP_2666 ← / ← PP_2667 hypothetical protein/ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,054,7780TA14.6% 75.0 / 4.9 33intergenic (‑253/+55)PP_2666/PP_2667hypothetical protein/ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (10/17);  new base A (2/3);  total (13/20)
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

GAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGAAGGTGAGGATCGCCAGCAGGGTGAACAATGCGGTCAAGCCCA  >  NC_002947/3054634‑3054912
                                                                                                                                                |                                                                                                                                      
gAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATGGCGGCAGCATGCATAGATTGCTGGGGCCGCGTTGc                                                                                                                                    >  8:98994/1‑149 (MQ=255)
    gCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCATGCGTCGCGAAAGGGCCGCGCAGCGGCCCCCGGATCTCGGCCTAATGCATAGATTTCTGGGGCCGCATTTCGGcc                                                                                                                                >  1:212749/1‑149 (MQ=37)
      aTGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGccc                                                                                                                               >  6:306192/1‑148 (MQ=255)
       tGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCAAAGCGGcccca                                                                                                                             <  7:214635/149‑3 (MQ=255)
       tGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCAAAGCGGcccca                                                                                                                             <  4:171138/149‑3 (MQ=255)
                ccTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGacac                                                                                                                    <  5:306192/149‑1 (MQ=255)
                    gggAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAgg                                                                                                                <  4:295117/149‑1 (MQ=255)
                       aGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCg                                                                                                             <  7:98994/149‑1 (MQ=255)
                           aCTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATc                                                                                                                           <  2:208427/131‑1 (MQ=255)
                           aCTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGAGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATc                                                                                                                           >  1:208427/1‑131 (MQ=255)
                             tGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGcc                                                                                                              <  4:266253/142‑1 (MQ=255)
                             tGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCATGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCgagaacg                                                                                                       >  3:266253/1‑143 (MQ=255)
                                     tCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCAAAGCGGcccca                                                                                                                             <  8:23676/119‑3 (MQ=255)
                                     tCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCAAAGCGGcccca                                                                                                                             >  7:23676/1‑117 (MQ=255)
                                                   cTGGCGCGATCCCTGTGGGAGCGGCCATGCGTCGCGAAAAGGCCCCGCAGCGGCCCCAGGATATCGGCAACAAGCAAAGATTGCTGGGGGCGCTTTGAGGGCCAAAAGCGACACAAGGCCGCGCCCCCCCCGGTGACAAGCCacggtt                                                                                  >  8:150135/1‑143 (MQ=255)
                                                        gcgATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAgcg                                                                            <  3:170823/149‑1 (MQ=255)
                                                              ccTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCttt                                                                                                                                      >  6:40165/1‑85 (MQ=255)
                                                              ccTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCttt                                                                                                                                      <  5:40165/85‑1 (MQ=255)
                                                                  tgGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTgcg                                                                   <  2:209731/148‑1 (MQ=255)
                                                                          gCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGc                                                          <  1:204296/149‑1 (MQ=255)
                                                                            cTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCAt                                                        >  6:269569/1‑149 (MQ=255)
                                                                                     cgAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGt                                               >  6:25144/1‑149 (MQ=255)
                                                                                      gAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGtt                                              >  6:26045/1‑149 (MQ=255)
                                                                                       aaaGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGtt                                              >  6:16624/1‑148 (MQ=255)
                                                                                           ggCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGc                                                                        <  8:265371/118‑1 (MQ=255)
                                                                                             ccGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGAAGGTg                                       <  2:212749/149‑1 (MQ=255)
                                                                                              cgcgCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGAAGGTGa                                      <  7:73899/149‑1 (MQ=255)
                                                                                              cgcgCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGAAGGTGa                                      <  5:269569/149‑1 (MQ=255)
                                                                                                   aGCGGCCGCAGGATCGCGGCATCATGCATAGATTGCTGGGTCCGCTTTGCGGCCCAATCTCGACACATGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGTGGTTGAAGGTGAGGATc                                 <  2:308871/149‑1 (MQ=255)
                                                                                                           cAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGAAGGTGAGGATCGCcagcag                         >  4:12862/1‑149 (MQ=255)
                                                                                                                        tcatGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACGCGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGAAGGTGAGGATCGCCAGCAGGGTGAACAATGCg            <  7:150135/149‑1 (MQ=255)
                                                                                                                          atGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTg                                                    >  7:195118/1‑107 (MQ=255)
                                                                                                                          atGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTg                                                    <  8:195118/107‑1 (MQ=255)
                                                                                                                                 gATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGAAGGTGAGGATCGCCAGCAGGGTGAACAATGCGGTCAAGccc   <  5:16624/149‑1 (MQ=255)
                                                                                                                                  aTTGCTGGGGCCGCTTTGCGGCCAAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCATGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGAAGGTGAGGATCGCCAGCAGGGTGAACAATGCGGTCAAGCCCa  <  7:22396/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                      
GAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGAAGGTGAGGATCGCCAGCAGGGTGAACAATGCGGTCAAGCCCA  >  NC_002947/3054634‑3054912

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: