Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,557,343 C→T 17.6% T232T (ACG→ACA ddpF ← putative D,D‑dipeptide ABC transporter ATP‑binding subunit DdpF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,557,3430CT17.6% 48.1 / 3.2 28T232T (ACG→ACAddpFputative D,D‑dipeptide ABC transporter ATP‑binding subunit DdpF
Reads supporting (aligned to +/- strand):  ref base C (7/16);  new base T (3/2);  total (10/18)
Fisher's exact test for biased strand distribution p-value = 3.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

AACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCA  >  NC_000913/1557201‑1557485
                                                                                                                                              |                                                                                                                                              
aaCGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCATATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGcc                                                                                                                                                <  1:631498/143‑1 (MQ=255)
     tcgtcgTAAAACGAGACTTGCGGCAAGGGGCGGTTGCCCGGCAGATCCGTCTTTCGTGATGCCCATTCTTCCGCCAGCGGTTTATGAATGGCGGGGAGGGAATCCAGCAATAAGCGGGTGTATGGATGTGCAGGTGCCGTCAg                                                                                                                                           <  2:592703/143‑1 (MQ=255)
     tCGTAGAAAAACCGGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGTCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAg                                                                                                                                           <  2:504305/143‑1 (MQ=255)
          gAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACct                                                                                                                                      <  1:545896/143‑1 (MQ=255)
           aaagaaCAGCCTTGCGGCAACGCGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTTTCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCTGACGGCACctg                                                                                                                                     <  2:357460/143‑1 (MQ=255)
           aaaatgcAGCCTTGCGGGAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCTGACGGCACctg                                                                                                                                     <  2:658959/137‑1 (MQ=255)
               accAGCCTTGCGGCAACGTGCGGTTCCCCGGCCGCTCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGCGGGAATCCAGCAATAATCGCGTGTATGGATGTGCAGGTGCCGTCAGCACCtgctgc                                                                                                                                 <  1:511243/141‑1 (MQ=255)
                 cAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTgcgc                                                                                                                               <  1:477031/143‑1 (MQ=255)
                      ttGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGCAATCCAGCAATAATCGGGGGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGc                                                                                                                          >  2:110855/1‑143 (MQ=255)
                            cAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCACATGATTATATCGTTGAAGAATCACTGGCAGCCATTGTATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGtt                                                                                                                    <  2:660190/143‑1 (MQ=255)
                                gcgCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATGCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTAc                                                                                                                <  1:27311/141‑1 (MQ=255)
                                     gTTTCCCGGCAGATCCGTGTTACGTAATGCCCATTCTTCCTCCAGCGGATTATCAATGGCGCGGAGGGACTCCAGCAATCAGCGGGTGGATGGGTGCGCAGGGGCCGCCAGCACCAGACGCGCGGCGCCCGGTTCTCCAAGCt                                                                                                           >  2:261924/1‑143 (MQ=255)
                                           cagCAGATCCGTTTTACGGAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGAAGGGAATCCAGCAATAATCGGGGGTATGGATGTGCGGGTGCCGTCGGCACCTGCGGCACGTCGCCCAGTTGTACAGTCTGCCCga                                                                                                     >  2:555170/3‑143 (MQ=255)
                                                                 gtccATTCTCCTTCCAGCGTTCTGTCAATGTCGGGGACGGTATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCTAGATACATCACCGCCACCCGATc                                                                               <  2:382860/141‑1 (MQ=255)
                                                                 gCCCATTCTTCCTCCAGCGGTTTATCGATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGTATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATc                                                                               <  2:546726/143‑1 (MQ=255)
                                                                   ccATTCTTCCTCCATCGGTTTATCAATGGCGGGGAGGGAATCCCGCAATAATCGGGTGTATGGCTGTGCAGGTGCTGTCAGCAGCTGCTGCGCGTCGGCCAGTTCTACAATCTGCACGAGGCACATCGCGGCCGCCCgccagc                                                                             >  2:434199/1‑138 (MQ=255)
                                                                   ccATTCTGCCTCCAACAGTTGGTCGGTGGCGGGGAGGGAATCCAGCAATAAGCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGc                                                                             <  1:120647/143‑1 (MQ=255)
                                                                     aTTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGAATGGATGTGCAGCTGCCGTCAGCACCTGCCGGGAGTCGCCCAGTTCCACAATCTCCCCGGGCTGCGGCACCGCCACCCGAGCGCtg                                                                           >  2:77072/1‑142 (MQ=255)
                                                                      ttcttcGTCCAGCGGGTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCa                                                                          <  2:306082/143‑1 (MQ=255)
                                                                                    gTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGCAAAGTCCTTGCCTGTGGTTTTTGCGGGTCGCCCAGTTCTACAATCTGCCCGAGATACGTCACCGCCACCCGATCGCTCATATGACGTATcccc                                                            >  1:126586/1‑143 (MQ=255)
                                                                                     tttATCAATGGCGGGGCAGCAATCGAGCAATAATAGGGTGTATGGATGCGCAGGTGCTGTCAGCACGAGCGGCGCGTCGCCCAGTACTACAATCTGCCCGAGATACAACACCGGCGCGCGGTGGCGCATAGGCCATATCaacc                                                           >  2:468010/1‑140 (MQ=255)
                                                                                           aaTGGCCGAGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCGGCTGCGAGTCGCCCAGTTCTACAATCGGCCCGAGAAACAACACCGCCACCCGATCGCAGCTATGACGTAGCACCGAGAAAt                                                     >  2:116789/1‑143 (MQ=255)
                                                                                                    gagggaATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACAGCACCGCCACCCGATCGTTCATAGGACGGATCACGGAGAAACTGGGCGaaa                                            >  2:616960/1‑143 (MQ=255)
                                                                                                         aaTCCAGCAATAATCGTGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGc                                       <  1:626972/143‑1 (MQ=255)
                                                                                                                   taatCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCa                             <  1:35776/143‑1 (MQ=255)
                                                                                                                        gggTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACACTCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCg                        <  2:126586/143‑1 (MQ=255)
                                                                                                                             gaTGGAGGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGAtt                   <  2:650817/142‑1 (MQ=255)
                                                                                                                                         ggTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCAATCGCTCATATGAGGTATCAGCGAGACATTGCGAGAAAACAGCCCACAGGGCCGCACGTGAGTGTCCTggcacgg       >  2:491799/1‑137 (MQ=255)
                                                                                                                                            gCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCGCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTAc    <  2:156116/143‑1 (MQ=255)
                                                                                                                                              cGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCa  <  1:110855/143‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                              
AACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCA  >  NC_000913/1557201‑1557485

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