Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,371,967 G→A 9.5% G20D (GGC→GAC)  ycjN → putative ABC transporter periplasmic binding protein YcjN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,371,9670GA 9.5% 128.2 / 4.6 53G20D (GGC→GAC) ycjNputative ABC transporter periplasmic binding protein YcjN
Reads supporting (aligned to +/- strand):  ref base G (30/18);  new base A (3/2);  total (33/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.27e-01

TGATTAGTGAAGTTGATATATTGGGTAATGAAATAACGCTGCGCCCCTGGCAGGTTATGTGGATTAAATAAAAAGGAACATCTCATGATTAAATCAAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCG  >  NC_000913/1371825‑1372107
                                                                                                                                              |                                                                                                                                            
tGATTGGGGAAGTGGAGATATTGGGTAACGAAATAACGCTGCGCCCCTGGCAGGTTATGTGGATTAAATAAAAAGGAACATCTCATGATTAAATCAAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAgg                                                                                                                                              <  2:261296/143‑1 (MQ=255)
tGATTAGTGAAGTTGATATATTGGGTAATGAAATAACGCTGCGCCCCTGGCAGGTTATGTGGATTAAATAAAAAGGAACATCTCATGATTAAATCAAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGCTTTCAgg                                                                                                                                              >  2:85536/1‑143 (MQ=255)
       tcacggtGATATATTGGGTAATGAAATAACGCTGCGCCCCTGGCAGGTTATGTGGATTAAATAAAAAGGAACATCTCATGATTAAATCAAAAAGCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTaaa                                                                                                                                       <  1:74724/137‑1 (MQ=255)
       tGAAGTTGATATATTGGGTAATGAAATAACGCTGCGCCCCTGGCAGGTTATGTGGATTAAATAAAAAGGAACATCTCATGATTAAATCAAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTaaa                                                                                                                                       >  1:802795/1‑143 (MQ=255)
       tGAAGTTGATATATTGGGTAATGAAATAACGCTGCGCCCCTGGCAGGTTATGTGGATTAAATAAAAAGGAACATCTCATGATTAAATCAAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTaaa                                                                                                                                       >  1:74155/1‑143 (MQ=255)
          aGTTGATATATTGGGTAATGAACTAACGCTGCGCCCCTGGCAGGTTATGTGGATTAAATAAAAAGGAACATCTCATGATTAAATCAAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAaagaa                                                                                                                                    <  1:264316/143‑1 (MQ=255)
           gTTGATATATTGGGTAATGAAATAACGCTGCGCCCCTGGCAGGTTATGTGGATTAAATAAAAAGGAACATCTCATGATTAAATCAAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCCGGCTGTAaagaag                                                                                                                                   >  2:759716/1‑143 (MQ=255)
            ttGATATATTGGGTAATGAAATAACGCTGCGCCCCTGGCAGGTTATGTGGATTAAATAAAAAGGAACATCTCATGATTAAATCAAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAagaaga                                                                                                                                  >  1:535365/1‑143 (MQ=255)
               atatatTGGGTAATGAAATAACGCTGCGCCCCTGGCAGGTTATGTGGATTAAATAAAAAGGAACATCTCATGATTAAATCGAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAACGAGGgaaa                                                                                                                               >  1:496145/1‑143 (MQ=255)
                  gatTGAGTAATGAAATAACGCTGCGCCCCGGGCAGGATATGTGGATTAAATAAAAAGGAACATCTCATGATTAAATCAAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATaa                                                                                                                            <  2:683757/142‑1 (MQ=255)
                           aTGAAATAACGCTGCGCCCCTGGCAGGTTATGTGGATTAAATAAAAAGGAACATCTCATGATTAAATCAAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGt                                                                                                                   >  1:671575/1‑143 (MQ=255)
                                 tAACGCGGCGCGCCTCGCAGGTTATGTGGATTAAATAAAAAGGAACATCTCATGATTAAATCAAAAATCCTGCTGTTATCAGCACTGCTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCAt                                                                                                             <  1:104730/143‑1 (MQ=255)
                                   aCGCTGCGCCCCTGGCAGGTTATGTGGATTAAATAAAAAGGAACATCTCATGATTAAATCAAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCg                                                                                                           >  1:590939/1‑143 (MQ=255)
                                       tGCGCCCCTGGCAGGTTATGTGGATTAAATAGAAAGAAACATCTCATGATTAAATCAAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAAtt                                                                                                       <  2:256311/143‑1 (MQ=255)
                                           ccccTGGCAGGTTATGTGGATTAAATAAAAAGGAACATCTCATGATTAAATCAAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATg                                                                                                   <  1:382494/143‑1 (MQ=255)
                                           ccccTGGCAGGTTATGTGAATTAAATAAAAAGGAACATCTCATGATTAAATCAAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATg                                                                                                   <  1:489927/143‑1 (MQ=255)
                                               tGGCAGGTTATGTGGATTAAATAAAAAGGAACATCTCATGATTAAATCAAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCAtt                                                                                               >  1:106337/1‑143 (MQ=255)
                                                   aGGTTATGTGGATTAAATAAAAAGGAACATCTCATGATTAAATCAAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGGAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCAttcttc                                                                                           >  2:547641/1‑143 (MQ=255)
                                                          gtgGATTAAATAAAAAGGAACATCTCATGATTAAATCAAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCGTCGAATTTATGCATTCTTCGGTGGAg                                                                                    >  1:946160/1‑143 (MQ=255)
                                                           tgGATTAAATAAAAAGGAACATCTCATGATTAAATCAAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGc                                                                                   >  2:656352/1‑143 (MQ=255)
                                                               ttAAATAAAAAGGAACATCTCATGATTAAATCAAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGa                                                                               <  1:425994/143‑1 (MQ=255)
                                                               ttAAATAAAAAGGAACATCTCATGATTAAATCAAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGa                                                                               >  1:798410/1‑143 (MQ=255)
                                                               ttAAATAAAAAGGAACATCTCATGATTAAATCAAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGAGGACCAGGa                                                                               >  2:933886/1‑143 (MQ=255)
                                                               ttAAATAAAAAGGAACATCTCATGATTAAATCAAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATCTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGa                                                                               <  1:955930/143‑1 (MQ=255)
                                                                taaataaaAAGGAACATCTCATGATTAAATCAAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAg                                                                              <  1:773570/143‑1 (MQ=255)
                                                                     aaaaaGGAACATCTCATGATTAAATCAAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCa                                                                         >  1:40087/1‑143 (MQ=255)
                                                                           gAACATCTCATGATTAAATCAAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGt                                                                   <  1:94579/143‑1 (MQ=255)
                                                                              aagttcATGATTAAAGAAAAAATCGAGTAATTATAAGCAATGATTTCATGCTTAAATATTTAAGGCTGTAAATAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTAt                                                                <  2:537680/139‑1 (MQ=255)
                                                                               aTCTCATGATTAAATCAAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATc                                                               >  2:540687/1‑143 (MQ=255)
                                                                                  tcATGATTAAATCAAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGACTGTAACGAAGAAAATAACACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGt                                                            >  1:860937/1‑143 (MQ=255)
                                                                                  tcATGATTAAATCAAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGACTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGt                                                            >  2:986542/1‑143 (MQ=255)
                                                                                  tcATGATTAAATCAAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGACTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGt                                                            >  2:779186/1‑143 (MQ=255)
                                                                                         tAAATCAAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAattga                                                     <  1:540687/143‑1 (MQ=255)
                                                                                            aTCAAACATCGTGCTGGTATCGGCAATGGAATCATGCGCACTAATTTCAGGCTATAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAattgattg                                                  <  2:239241/143‑1 (MQ=255)
                                                                                               aaaaaTCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGACTGTATCCATCGCATTTATGCATTCTTCGGTGGACGAGTAGCGCTAGGCCGTTATCAGTCAACTGTTTGccc                                               >  2:233637/1‑143 (MQ=255)
                                                                                                   aTCGTGCCGTAAACAGAACAAGTTAAATACGCCATGAATACAGGCTCTAAATAAAAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGttt                                           <  2:719725/143‑1 (MQ=255)
                                                                                                         cTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAACTTGATTGCCCGTTTTGcaaa                                     >  2:405625/1‑143 (MQ=255)
                                                                                                            tgaTCAGCACTGGGTTCATGAGCCCCGAGCTGCAGCAGTAAGGAAGAGAATAAAACGAAGGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGa                                  <  1:179857/141‑1 (MQ=255)
                                                                                                              aTCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGCGCAGGAGCGCCACGCCGTTATCAGCAAATTGATTGCCCGTTTTGCAAAAGaac                                >  2:510714/1‑142 (MQ=255)
                                                                                                                  gCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGACAGGCCGTTATCAGTAAATTGATTGGCCGTTTTGAAAAAGACAAccc                            >  2:218201/1‑143 (MQ=255)
                                                                                                                   cACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCt                           >  1:644825/1‑143 (MQ=255)
                                                                                                                   cACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTACCAGCAAAATGATCGCCCGTTTTGAAAAAGAAAAcccg                           >  1:765816/1‑142 (MQ=255)
                                                                                                                        gTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGACTGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAAGACGATATCAGTAAATTGATTGCCCGCTTCGAAAAAGAAAGCACTGGCac                      >  2:431842/1‑142 (MQ=255)
                                                                                                                         tttCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGTCCGTCTTGAAAAAGAAAAACCTGGCATc                     >  1:176881/1‑143 (MQ=255)
                                                                                                                          ttCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTCTGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTGATCAGCAAATTGATTACCAGTTTTGAAGAAGGAAAGCCTGGCATcc                    >  1:236541/1‑142 (MQ=255)
                                                                                                                           tcaggagcCTTGATTTCAGGAGGTAAGGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCAc                   <  1:836130/137‑1 (MQ=255)
                                                                                                                              tGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGCAAAAGAAAACCCTGGCATCACCGt                >  2:709451/1‑143 (MQ=255)
                                                                                                                               gcgcCCTGATTTCAGGCTGTACAGCAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCATGCGCGCCAGGCCGATATCAGTAAATTGATTGCCGGTTTTGCAACAGAAAACCCTGGCATCACCggt               >  2:43203/1‑141 (MQ=255)
                                                                                                                                  cccTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTGTCGGTGAATTGATTGTCCGTTTTGGAAATGTAAGCCCTGGCATCAGCGTTAAg            >  1:295435/1‑143 (MQ=255)
                                                                                                                                  cccTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAg            >  1:368638/1‑143 (MQ=255)
                                                                                                                                    cTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCa          >  1:465218/1‑143 (MQ=255)
                                                                                                                                           gaGGCGGTAAGGGAGAAAATAAAACGAAGGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGccc   <  2:259820/142‑1 (MQ=255)
                                                                                                                                           cAGACTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCTTTAAGCAAGTGccc   <  1:779186/143‑1 (MQ=255)
                                                                                                                                           cAGACTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGccc   <  1:986542/143‑1 (MQ=255)
                                                                                                                                            aGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCATACGTGGAGCGGTAGGGCCGGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCg  >  1:270151/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                            
TGATTAGTGAAGTTGATATATTGGGTAATGAAATAACGCTGCGCCCCTGGCAGGTTATGTGGATTAAATAAAAAGGAACATCTCATGATTAAATCAAAAATCGTGCTGTTATCAGCACTGGTTTCATGCGCCCTGATTTCAGGCTGTAAAGAAGAAAATAAAACGAATGTATCCATCGAATTTATGCATTCTTCGGTGGAGCAGGAGCGCCAGGCCGTTATCAGTAAATTGATTGCCCGTTTTGAAAAAGAAAACCCTGGCATCACCGTTAAGCAAGTGCCCG  >  NC_000913/1371825‑1372107

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