Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,675,336 T→C 13.3% intergenic (+35/+69) PP_4137 → / ← chrR outer membrane siderophore receptor/chromate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,675,3360TC13.3% 67.7 / 3.8 30intergenic (+35/+69)PP_4137/chrRouter membrane siderophore receptor/chromate reductase
Reads supporting (aligned to +/- strand):  ref base T (9/17);  new base C (2/2);  total (11/19)
Fisher's exact test for biased strand distribution p-value = 6.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.19e-01

TGCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATT  >  NC_002947/4675193‑4675482
                                                                                                                                               |                                                                                                                                                   
tGCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGaa                                                                                                                                                  <  8:164562/147‑1 (MQ=255)
 gCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAAttt                                                                                                                                               >  8:287036/1‑149 (MQ=255)
            gACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCAtt                                                                                                                                    >  3:17157/1‑149 (MQ=255)
             aCGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCAtt                                                                                                                                    <  7:287036/148‑1 (MQ=255)
               gCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTcgc                                                                                                                                 <  4:17157/149‑1 (MQ=255)
                ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCGGCAAATTCAATGACGCATTcgcg                                                                                                                                <  7:22750/149‑1 (MQ=37)
                 cGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCgg                                                                                                                               <  8:325827/149‑1 (MQ=255)
                 cGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCgg                                                                                                                               <  5:135528/149‑1 (MQ=255)
                           cgacgCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCGGCAAATTCAATGACGCATTCGCGGGt                                                                                                                             <  8:66503/141‑1 (MQ=37)
                           cgacgCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCGGCAAATTCAATGACGCATTCGCGGGt                                                                                                                             >  7:66503/1‑141 (MQ=37)
                             acgCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTc                                                                                                                   <  8:261885/149‑1 (MQ=255)
                                            aTGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTg                                                                                                                                                    >  1:272346/1‑101 (MQ=255)
                                            aTGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTg                                                                                                                                                    <  2:272346/101‑1 (MQ=255)
                                                                gAAGGATTTCGCGTGCATGCCGGGGTTTCCTACAGCGACCGCTGACGGGGGGGCGCGGGTTTACCCGCGCATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCGCCTACACAGGGGCGTTGCACAGCAGGATTAgggggt                                                                                <  8:174875/149‑3 (MQ=38)
                                                                          gCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGTGTTTACCCGTGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGcc                                                                       <  1:109932/148‑1 (MQ=255)
                                                                               tATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGa                                                                            <  5:319707/138‑1 (MQ=255)
                                                                               tATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTc                                                                 <  3:157935/149‑1 (MQ=255)
                                                                                  gCCGGGGTTTCCTACCGCTTCCGCGGACGGGGGGGCGCGGGTTTACCCGCTAATGCGTCATTGAATTTGCCGCCGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCaa                                                               <  3:60119/148‑1 (MQ=38)
                                                                                                      ccGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAAc                                            <  7:61290/147‑1 (MQ=255)
                                                                                                       cGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGAATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAACAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCt                                          <  5:18954/148‑1 (MQ=255)
                                                                                                          tGACGGGGGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCTACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCg                                      <  7:214370/149‑1 (MQ=255)
                                                                                                                          gTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAg                      >  8:301108/1‑148 (MQ=255)
                                                                                                                              aCCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGt                                        >  4:99085/1‑127 (MQ=255)
                                                                                                                              aCCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGt                                        <  3:99085/127‑1 (MQ=255)
                                                                                                                                ccGCGAATGCGTCGGCAAATTCAATGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTc               >  6:31295/1‑149 (MQ=37)
                                                                                                                                ccGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTc               >  6:280945/1‑149 (MQ=255)
                                                                                                                                ccGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTc               >  3:316372/1‑149 (MQ=255)
                                                                                                                                   cgAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcgcg            >  2:209617/1‑149 (MQ=255)
                                                                                                                                         cGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGCACAAACGCCTGCAAGAACGGTCGCGTCtt        <  1:209617/148‑1 (MQ=255)
                                                                                                                                           tCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCa    >  6:330613/1‑149 (MQ=255)
                                                                                                                                              ttGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCtt        <  3:81774/142‑1 (MQ=255)
                                                                                                                                              ttGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCtt        >  4:81774/1‑142 (MQ=255)
                                                                                                                                              ttGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCAtt  >  6:120436/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
TGCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATT  >  NC_002947/4675193‑4675482

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: