Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,029,037 C→T 11.1% intergenic (+66/‑373) PP_4432 → / → PP_4433 Xaa‑Pro aminopeptidase/amino acid MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,029,0370CT11.1% 91.3 / 3.4 36intergenic (+66/‑373)PP_4432/PP_4433Xaa‑Pro aminopeptidase/amino acid MFS transporter
Reads supporting (aligned to +/- strand):  ref base C (17/15);  new base T (2/2);  total (19/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

GATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATCAAGGACGACGCTGAATAGAAGTAA  >  NC_002947/5028890‑5029184
                                                                                                                                                   |                                                                                                                                                   
gaTTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCtt                                                                                                                                                   >  1:233283/1‑150 (MQ=255)
gaTTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCtt                                                                                                                                                   >  1:414426/1‑150 (MQ=255)
     ccGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGc                                                                                                                                              >  1:662593/1‑150 (MQ=255)
        agagCATCCTGATAACTGATACTGGGATCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCAcc                                                                                                                                           <  1:1076184/150‑1 (MQ=255)
                 tGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCTTTAAGGCACCAt                                                                                                                                         <  2:787913/143‑1 (MQ=255)
                 tGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCTTTAAGGCACCAt                                                                                                                                         >  1:787913/1‑143 (MQ=255)
                 tGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCa                                                                                                                                  >  1:311737/1‑150 (MQ=255)
                        tGATACTGGGGCCGACACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCCCCATCACTCCAGCCACTc                                                                                                                           >  2:457801/1‑150 (MQ=255)
                           tACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGCCAATGGCTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCtt                                                                                                                        <  1:965777/150‑1 (MQ=255)
                             cTGCGGTCGAGACGTTATGAGCGCCTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTg                                                                                                                      <  2:233283/150‑1 (MQ=255)
                              tGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGt                                                                                                                     <  2:311737/150‑1 (MQ=255)
                                               gCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACg                                                                                                    >  2:464400/1‑150 (MQ=255)
                                                        ctcGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCTTTAAGGCACCATCa                                                                                                                                       <  2:1015223/106‑1 (MQ=255)
                                                        ctcGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCTTTAAGGCACCATCa                                                                                                                                       >  1:1015223/1‑106 (MQ=255)
                                                        ctcGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGa                                                                                           >  1:660826/1‑150 (MQ=255)
                                                         tcGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGaa                                                                                          >  1:84598/1‑150 (MQ=255)
                                                            cAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAagag                                                                                       <  1:842156/150‑1 (MQ=255)
                                                             aaGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAgg                                                                                      <  2:662593/150‑1 (MQ=255)
                                                                                 aCGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTcccc                                                                  >  1:118370/1‑150 (MQ=255)
                                                                                  cGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCa                                                                 <  2:118370/150‑1 (MQ=255)
                                                                                   gAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTaaa                                                                                                                                                <  1:327655/70‑1 (MQ=255)
                                                                                   gAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTaaa                                                                                                                                                >  2:327655/1‑70 (MQ=255)
                                                                                   gAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCttt                                                                                                                       <  2:697311/95‑1 (MQ=255)
                                                                                   gAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCttt                                                                                                                       >  1:697311/1‑95 (MQ=255)
                                                                                          ttGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGa                                                         <  2:660826/150‑1 (MQ=255)
                                                                                                                   cACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGaaa                                <  1:380171/150‑1 (MQ=255)
                                                                                                                    aCGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGa                                                                           >  2:140148/1‑106 (MQ=255)
                                                                                                                    aCGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGa                                                                           <  1:140148/106‑1 (MQ=255)
                                                                                                                      ggCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACtt                             >  1:464432/1‑150 (MQ=255)
                                                                                                                       gCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACCGCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCTTAACGCCTGAGACTTa                            <  2:911009/150‑1 (MQ=255)
                                                                                                                            cTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATCAAg                       >  2:201541/1‑150 (MQ=255)
                                                                                                                            cTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATCAAg                       >  1:992795/1‑150 (MQ=255)
                                                                                                                              tCAGCGACAATGGGATGGAGCCTTAAAGCCCCAGCACTGCAGCCCCTGCTATGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGGCTGAAACTTATCAAGga                     <  2:586121/150‑1 (MQ=255)
                                                                                                                                     cAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATCAAGGACGACGCt              >  1:552631/1‑150 (MQ=255)
                                                                                                                                             tGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATCAAGGACGACGCTGAATAGaa      <  1:719278/150‑1 (MQ=255)
                                                                                                                                                 gCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATCAAGGACGACGCTGAATAGAAGTaa  >  1:314070/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATCAAGGACGACGCTGAATAGAAGTAA  >  NC_002947/5028890‑5029184

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: