Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,029,042 A→G 12.1% intergenic (+71/‑368) 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,0420AG12.1% 76.4 / 3.4 33intergenic (+71/‑368)PP_4432/PP_4433Xaa‑Pro aminopeptidase/amino acid MFS transporter
Reads supporting (aligned to +/- strand):  ref base A (14/15);  new base G (2/2);  total (16/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 = 5.86e-01

CCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATCAAGGACGACGCTGAATAGAAGTAATACC  >  NC_002947/5028895‑5029188
                                                                                                                                                   |                                                                                                                                                  
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)
                                                                                                                                                tAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATCAAGGACGACGCTGAATAGAAGTAATAcc  <  2:552631/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTGTTCCCCAGGGAGCGAGTCAACCTGTGCCTAACGCCTGAAACTTATCAAGGACGACGCTGAATAGAAGTAATACC  >  NC_002947/5028895‑5029188

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: