Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,559,761 T→C S50S (TCA→TCG)  ddpC ← D,D‑dipeptide ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,559,7610TC98.1% 177.7 / ‑4.3 52S50S (TCA→TCG) ddpCD,D‑dipeptide ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (23/28);  total (23/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATCGCGTTTGTTTTTGTGGGCGTACGGCGGACGTTTCCTCG  >  NC_000913/1559614‑1559899
                                                                                                                                                   |                                                                                                                                          
ggCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGa                                                                                                                                          >  1:752929/1‑150 (MQ=255)
  cGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGaaa                                                                                                                                        >  2:613624/1‑150 (MQ=255)
   gagaATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGaaaaa                                                                                                                                      <  1:504513/151‑1 (MQ=255)
   gagaATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGAGCATGCGTCGCCAGCCACTGCGaaaaa                                                                                                                                      <  1:715637/151‑1 (MQ=255)
    agaATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGaaaaaa                                                                                                                                     <  1:483443/151‑1 (MQ=255)
       attgattgCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAatca                                                                                                                                  <  2:232644/151‑1 (MQ=255)
       attgattgCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAatc                                                                                                                                   >  2:176353/1‑150 (MQ=255)
         tgattgCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAatcat                                                                                                                                 <  2:450455/150‑1 (MQ=255)
            ttgCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAatcatcatc                                                                                                                             >  1:33840/1‑151 (MQ=255)
                 gACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAg                                                                                                                        <  2:773288/151‑1 (MQ=255)
                     gCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCa                                                                                                                      <  2:58063/149‑1 (MQ=255)
                     gCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCAt                                                                                                                     <  2:104861/150‑1 (MQ=255)
                      ccGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTa                                                                                                                   <  1:34621/151‑1 (MQ=255)
                       cGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTaa                                                                                                                  <  1:672583/151‑1 (MQ=255)
                           cAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATAc                                                                                                               >  2:353585/1‑150 (MQ=255)
                            aGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGAGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACaa                                                                                                             <  1:18045/151‑1 (MQ=255)
                               aCGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAatga                                                                                                          <  1:229715/151‑1 (MQ=255)
                                 gcgGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAatgatg                                                                                                        <  2:141982/151‑1 (MQ=255)
                                   ggCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGAc                                                                                                      >  2:124584/1‑151 (MQ=255)
                                    gCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGaaaaaa                                                                                                                                     <  2:741752/119‑1 (MQ=255)
                                    gCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGaaaaaa                                                                                                                                     >  1:741678/1‑119 (MQ=255)
                                           aGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTgg                                                                                              >  2:771876/1‑151 (MQ=255)
                                              tCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTc                                                                                            >  2:706959/1‑150 (MQ=255)
                                               cgcgTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCAc                                                                                          >  1:191040/1‑151 (MQ=255)
                                                gcgTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACg                                                                                         >  1:588263/1‑151 (MQ=255)
                                                    cccACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGa                                                                                     <  2:148424/151‑1 (MQ=255)
                                                      cACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGc                                                                                   <  2:91316/151‑1 (MQ=255)
                                                          tCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGAc                                                                               >  2:34197/1‑151 (MQ=255)
                                                           cATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGAc                                                                               <  2:191073/150‑1 (MQ=255)
                                                            aTCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTg                                                                             <  2:26349/151‑1 (MQ=255)
                                                                  gCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCttt                                                                        <  2:694588/150‑1 (MQ=255)
                                                                      aaCCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGc                                                                    >  1:389954/1‑150 (MQ=255)
                                                                        ccAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATc                                                                 >  1:720669/1‑151 (MQ=255)
                                                                              cgccgcAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATAcc                                                           <  2:345330/151‑1 (MQ=255)
                                                                                   cAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGttt                                                      >  1:414516/1‑151 (MQ=255)
                                                                                    aGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGtttt                                                     <  2:594243/151‑1 (MQ=255)
                                                                                      aaGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTg                                                    <  1:340885/150‑1 (MQ=255)
                                                                                                 aGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATc                                        >  1:220418/1‑151 (MQ=255)
                                                                                                        cggTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATCGCgtttg                                 >  1:588515/1‑151 (MQ=255)
                                                                                                         ggTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATCGCgtttgt                                <  2:708618/151‑1 (MQ=255)
                                                                                                                tCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATCGCGTTTGTTTTtgtg                         >  2:18046/1‑151 (MQ=255)
                                                                                                                 cAATGGCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATCGCGTTTGTTTTTGTgg                        >  1:192734/1‑151 (MQ=255)
                                                                                                                   aTGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCTTTAATACAATGATGACTCCGCCGGTCACTGTTAGCGTACTGCCTTTCAGCATCCATACCATTTTTGCACCGTTAACTCGCGTTTGttttt                            <  1:523874/145‑1 (MQ=255)
                                                                                                                     ggCGTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATCGCGTTTGTTTTTGTGGGCGt                    >  2:137511/1‑151 (MQ=255)
                                                                                                                        gTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATCGCGTTTGTTTTTGTGGGCGTAcg                 <  2:421718/151‑1 (MQ=255)
                                                                                                                        gTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATCGCGTTTGTTTTTGTGGGCGTAcg                 >  1:166386/1‑151 (MQ=255)
                                                                                                                        gTTGGGATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATCGCGTTTGTTTTTGTGGGCGTAcg                 >  1:494762/1‑151 (MQ=255)
                                                                                                                           gggATCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATCGCGTTTGTTTTTGTGGGCGTAcggc               <  1:698497/150‑1 (MQ=255)
                                                                                                                              aTCATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCTCCGTTAAATCGCGTTTGTTTTTGTGGGCGTACGGCGGAc           <  2:631371/151‑1 (MQ=255)
                                                                                                                                cATGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATCGCGTTTGTTTTTGTGGGCGTACGGCGGACGt         >  2:167952/1‑151 (MQ=255)
                                                                                                                                 aTGCGTCGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATCGCGTTTGTTTTTGTGGGCGTACGGCGGACGtt        <  1:166450/151‑1 (MQ=255)
                                                                                                                                       cGCCAGCCACGGCGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATCGCGTTTGTTTTTGTGGGCGTACGGCGGACGTTTCCTCg  <  2:202714/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                          
GGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACGGTGAGCGGACTGCCTTTCAGCATCCATACCAGTTTTGCACCGTTAAATCGCGTTTGTTTTTGTGGGCGTACGGCGGACGTTTCCTCG  >  NC_000913/1559614‑1559899

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑