Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,675,342 T→C 13.3% intergenic (+41/+63) 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,3420TC13.3% 66.4 / 3.5 30intergenic (+41/+63)PP_4137/chrRouter membrane siderophore receptor/chromate reductase
Reads supporting (aligned to +/- strand):  ref base T (10/16);  new base C (2/2);  total (12/18)
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

GCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCT  >  NC_002947/4675194‑4675487
                                                                                                                                                    |                                                                                                                                                  
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)
                                                               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                                             <  3:99085/127‑1 (MQ=255)
                                                                                                                             aCCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGt                                             >  4:99085/1‑127 (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)
                                                                                                                                                 aTTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCt  >  4:175797/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCT  >  NC_002947/4675194‑4675487

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: