Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,063,221 T→G 11.6% intergenic (+78/‑129) PP_1838 → / → PP_1839 membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,063,2210TG11.6% 75.4 / 3.2 34intergenic (+78/‑129)PP_1838/PP_1839membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (13/16);  new base G (2/2);  total (15/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.05e-01

TATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGT  >  NC_002947/2063075‑2063359
                                                                                                                                                   |                                                                                                                                          
tataTCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTcc                                                                                                                                          <  3:31846/149‑1 (MQ=255)
   atCCAGACCGCGTCCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCtgt                                                                                                                                       >  5:298444/1‑149 (MQ=255)
           cgcgACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCgg                                                                                                                               <  3:251098/149‑1 (MQ=18)
           cgcgACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCgg                                                                                                                               <  1:36059/149‑1 (MQ=255)
              gACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTTCCCGCTACCACAAGCAATTCACACGGTTCGGACCTTGTGCAGGTCCTCTGCGAGCCGGc                                                                                                                             >  3:250450/1‑148 (MQ=255)
                 cAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCGTGCCCGCGCCCACAGGCATTGCACACGGGTCGGACCGTGTGGAGGTGCGGTGCCAGCGGGCAAGc                                                                                                                         >  1:135659/1‑149 (MQ=255)
                     agCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCcgc                                                                                                                     <  6:298444/149‑1 (MQ=255)
                          gctgGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATg                                                                                                                >  8:80293/1‑149 (MQ=255)
                           ctgGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTGGCCCGCTCCCACCGCCATTGCACACGGTTCGGCCCTTGTCCTGGTCCGGTGCGAGCGGGCAAGCCCGCGAATg                                                                                                                <  5:192616/148‑1 (MQ=255)
                                   cacCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGca                                                                                                       >  3:88228/1‑149 (MQ=255)
                                    acCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCGGCCCCAGGCAGTGGCCAGGGGTGGGGCCTTGTGGCGGTCCCGTGCGGGCGGGCAAGCCCGCGCAAGGGGCCCggc                                                                                                       >  3:90760/1‑146 (MQ=255)
                                             gtggttCAGGGAACACCGTATTAACCCGGTGCGGCGTTGCGGCCCATTCGCGGGTTTGCCCGCTCCCACAGGCTTTGCACAGTGTTCGGACTTTGTGCTGTTCGTGTGGGAGCGGGAAAGCCCGCGAATGGGCCCCGCACAATTTTccc                                                                                             <  4:250450/149‑1 (MQ=255)
                                                 tgCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGAc                                                                                         >  2:42590/1‑149 (MQ=255)
                                                     cGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGaa                                                                                     <  1:42590/149‑1 (MQ=255)
                                                            ccGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATc                                                                              <  7:80293/149‑1 (MQ=255)
                                                             cGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGc                                                                                                                             <  2:319423/101‑1 (MQ=255)
                                                             cGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGc                                                                                                                             >  1:319423/1‑101 (MQ=255)
                                                             cGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATcc                                                                             <  3:291578/149‑1 (MQ=255)
                                                                    aaCCGGGGCCGCGTTGCGGCCCATTAGCGGGCATTCCCCGAGCCCCCCGGCATGCACAGGGTTAGGGCCATTTTCAGGGCCGGGGCGGCCCGGGAAGGCCGCGGATGGGGCCCGCACAATTTTACCGGCCGGGACATTTTAcccccccc                                                                      >  2:69044/1‑149 (MQ=255)
                                                                           gCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGCCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCa                                                               <  3:121152/149‑1 (MQ=255)
                                                                                   gCGGCCCATCCGCGGGCTTGCCCGCGCCCACCGCCCTTGCACCGGGTCCGCACCTCGTCCAGGACCTGGGCTCGCGGGCAACCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCa                                                       <  6:192593/149‑1 (MQ=255)
                                                                                                  gCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTg                                        <  4:88228/149‑1 (MQ=255)
                                                                                                   cTTGCCCGCTCCCTCAGGCATTGCACAGGGTGCGGACCTTGTGCAGGTCCTGTGCGAGCGTGCAAGCCCGCGAATGTGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTgg                                       <  2:43128/149‑1 (MQ=255)
                                                                                                           gCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTACCGGACTGAACATTATCCCCGCCCGAAGGTCAGTTATTCAGGCATAAAACACCTGGGGTTCGt                                >  6:90969/1‑149 (MQ=39)
                                                                                                                cacaGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACt                           >  4:49478/1‑149 (MQ=255)
                                                                                                                       cctGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTga                         >  1:322497/3‑144 (MQ=25)
                                                                                                                       cctGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTga                         <  2:322497/142‑1 (MQ=25)
                                                                                                                          gCACAGGGTTCGGACTTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATca                                                <  1:130919/118‑1 (MQ=255)
                                                                                                                          gCACAGGGTTCGGACTTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATca                                                >  2:130919/1‑118 (MQ=255)
                                                                                                                            acaGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTc                                                        <  5:45314/108‑1 (MQ=255)
                                                                                                                            acaGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTc                                                        >  6:45314/1‑108 (MQ=255)
                                                                                                                                  ttCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGgt                                      <  7:252649/120‑1 (MQ=255)
                                                                                                                                  ttCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGgt                                      >  8:252649/1‑120 (MQ=255)
                                                                                                                                         cTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGt  <  3:49478/149‑1 (MQ=255)
                                                                                                                                                   tCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTcc                 <  2:101989/124‑1 (MQ=255)
                                                                                                                                                   tCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTcc                 >  1:101989/1‑124 (MQ=255)
                                                                                                                                                   |                                                                                                                                          
TATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAATGAAAGAGT  >  NC_002947/2063075‑2063359

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: