Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,549,605 T→A 13.4% intergenic (+18/‑154) PP_t47 → / → PP_4881 tRNA‑Leu/iron ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,549,6050TA13.4% 64.6 / 3.9 30intergenic (+18/‑154)PP_t47/PP_4881tRNA‑Leu/iron ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (10/16);  new base A (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 = 5.19e-01

AACGATAGTCGCCTTCTGCAGATGATGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATGTAGCCGTAAGCGAGG  >  NC_002947/5549465‑5549748
                                                                                                                                            |                                                                                                                                               
aaCGATAGTCGCCTTCTGCAGATGATGTTAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCAAGCTAGTg                                                                                                                                         >  1:193770/1‑149 (MQ=255)
aaCGATAGTCGCCTTCTGCAGATGATGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTg                                                                                                                                         <  3:141572/149‑1 (MQ=255)
 aCGATAGTCGCCTTCTGCAGATGATGTTAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCAAGCTAGTgg                                                                                                                                        <  2:193770/149‑1 (MQ=255)
    aTAGTCGCCTTCTGCAGATGATGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGtt                                                                                                                                     >  3:79015/1‑149 (MQ=255)
                           tAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTgg                                                                                                                                        >  1:84935/1‑123 (MQ=255)
                           tAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTgg                                                                                                                                        <  2:84935/123‑1 (MQ=255)
                           tAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCAAGCTAGt                                                                                                                                          >  6:230804/1‑121 (MQ=255)
                           tAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCAAGCTAGt                                                                                                                                          <  5:230804/121‑1 (MQ=255)
                             aatgtgtGCCCAGGTGGCGGAAGTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCAAGTGGGTTTTTTCTTGCCTGCGATTTATGTGc                                                                                                             <  6:40847/142‑1 (MQ=255)
                                                                  gTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGATAGTGGGTTTTTTCTTGCCTGCGATTTAt                                                                                                                 >  8:218492/1‑107 (MQ=255)
                                                                      catgTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTtctc                                                                    <  1:334867/145‑1 (MQ=255)
                                                                      cAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTc                                                                                                       <  8:188782/113‑1 (MQ=255)
                                                                               gTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTTGTCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATTGTCTGACCATGAATTTCTCGTCACACGGc                                                          <  1:102271/149‑1 (MQ=255)
                                                                               gTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGGCCCCGCCGCCCCTGGCACCATGGTCGGCCCATGa                                                                           >  4:332954/1‑132 (MQ=38)
                                                                               gTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGa                                                                           <  3:332954/132‑1 (MQ=255)
                                                                               gTGACAGCAAGGTCGTGGAAGTTCGAGTCATCTCCTGGGCACCATACAAAAACCAACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGGCACGGGCCCCCCGTGCACCATGTGCTGAAAATCAAGATCTCTTCacacag                                                           >  2:102271/1‑146 (MQ=25)
                                                                                          gTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGaaa                                               <  7:2278/149‑1 (MQ=255)
                                                                                              tGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCGCCCTCCCCCCCTTGCACCATGGGCTGACCATGAAGTTAGAGTCGCACGGCCAGTTGTGAAACAAt                                           >  8:260761/1‑149 (MQ=255)
                                                                                                   gTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGt                                      >  6:45289/1‑149 (MQ=255)
                                                                                                       gAGTCTTCTCCTTGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTAc                                  <  3:236404/149‑1 (MQ=255)
                                                                                                       gAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGAACATGAATTTCACGTCAAACTGCCACTTGAGAAACAAATTcatata                                   >  4:236404/1‑143 (MQ=255)
                                                                                                        aGTCATCTCCAGGGCACCATTCCAAAACCCCCTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTAcc                                 <  8:82689/149‑1 (MQ=255)
                                                                                                            ttgtcCTTGTCACCATACAAATACCCACTATCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCGCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTc                                        <  7:73603/135‑1 (MQ=255)
                                                                                                                    ggCACCATACAAAAACCCACTATCTTGCTAGTGGTTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTAcc                                 <  7:148404/137‑1 (MQ=255)
                                                                                                                             cAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATGTAGcc            <  2:160704/149‑1 (MQ=255)
                                                                                                                                     cACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGCTTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATGTAGCCGTAAGCGa    <  4:79015/149‑1 (MQ=255)
                                                                                                                                     cACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGtt                                     >  2:283283/1‑116 (MQ=255)
                                                                                                                                     cACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATGTAGCCGTAAGCGa    <  5:45289/149‑1 (MQ=255)
                                                                                                                                     cACTAGCTTGCTAGTGGGTTATTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGtt                                     <  1:283283/116‑1 (MQ=255)
                                                                                                                                        tAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCGCCCCCCCGTGCCCCATGGTCTGCCCATGAATTTCTCGTCACACGGGCACTTGAGAAACAATTTGGTTTAGCATTGTTTGCATATTTGTAGCCGTAATGGAgg  >  8:153959/1‑148 (MQ=255)
                                                                                                                                            |                                                                                                                                               
AACGATAGTCGCCTTCTGCAGATGATGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATGTAGCCGTAAGCGAGG  >  NC_002947/5549465‑5549748

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: