Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,376,207 C→A 9.1% A58A (GCC→GCA ycjR → putative ketohexose‑3‑epimerase YcjR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,376,2070CA 9.1% 124.9 / 4.4 52A58A (GCC→GCAycjRputative ketohexose‑3‑epimerase YcjR
Reads supporting (aligned to +/- strand):  ref base C (25/22);  new base A (2/3);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.39e-01

TTCCGGAAAACATTCTGGAGAAATTTCGTTATATCAAAGAGATGGGCTTCGATGGTTTTGAGATTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGG  >  NC_000913/1376065‑1376340
                                                                                                                                              |                                                                                                                                     
ttCCGGAAAACATTCTGGAGAAATTTCGTTATATCAAAGAGATGGGCTTCGATGGTTTGGAGATTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGGGACCACCGcc                                                                                                                                       >  2:525179/1‑143 (MQ=255)
  ccGGAAAACATTCTGGAGAAATTTCGTTATATCAAAGAGATGGGCTTCGATGGTTTTGAGATTGACGGCAAACTGCTGGTTAACAACATCGACGCAGACAAAGCGGCGCTCGAAGAACCCGGCTTACCGGTGACCACCGCCTg                                                                                                                                     >  1:656629/1‑143 (MQ=255)
   cGGAAAACATTCTGGAGAAATTTCGTTATATCAAAGAGATGGGCTTCGATGGTTTTGAGATTAACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAAAAACCGGTTTACCGGTGACCACCGCATGc                                                                                                                                    >  2:8665/1‑143 (MQ=255)
     gAAAACATTATGGAGAAATTTCGTTATATCAAAGAGATGGGCTTCGATGGTTTTGAGATTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCgg                                                                                                                                  >  1:832000/1‑143 (MQ=255)
       aaaCATTCTGGAGAAATTTCGTTATATCAAAGAGATGGGCTTCGATGGTTTTGAGATTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTg                                                                                                                                >  2:199734/1‑143 (MQ=255)
       aaaCATTCTGGAGAAATTTCGTTATATCAAAGAGATGGGCTTCGATGGTTTTGAGATTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATACAAGACACCGGTTTACCGGTGACCACCGCCTgggggg                                                                                                                                >  1:382492/1‑141 (MQ=255)
         aggtgcTGGCGAAATTTCGTTATATCAAAGAGATGGGCTTCGATGGTTTTGAGATTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCAGGCGGTggt                                                                                                                              <  2:832415/138‑2 (MQ=255)
         aCATTCTGGAGAAATTTCGTTATATCAAAGAGATGGGCTTCGATGGTTTTGAGATTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCAGGCGGTggt                                                                                                                              <  2:903534/143‑2 (MQ=255)
            ttCTGGAGAAATTTCGTTATATCAAAGAGATGGGCTTCGATGGTTTTGAGATTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTAt                                                                                                                           >  2:818105/1‑143 (MQ=255)
                  agaAATTTCGTTATATCAAAGAGATGGGCTTCGATGGTTTTGAGATTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACgga                                                                                                                     >  1:158375/1‑143 (MQ=255)
                  agaAAGTTCTTTATATCAAAGAGATGGGCTTCGATGGTTTTGAGATTGACGGCAAACTGCTGTTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACgga                                                                                                                     <  1:531186/143‑1 (MQ=255)
                           gTTATATCGCAGAGAGGGGCTGCGGTGGTTTGCAGCTCGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAGCCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGc                                                                                                            <  2:34442/143‑1 (MQ=255)
                              atatCAAAGAGATGGGCTTCGATGGTTTTGAGATTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGAc                                                                                                         >  1:887031/1‑143 (MQ=255)
                              atatCAAAGAGATGGGCTTCGATGGTTTTGAGATTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGAc                                                                                                         >  1:643668/1‑143 (MQ=255)
                                    atgagaTGGGCTTTAATGTTTTTGAGATTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATc                                                                                                   <  2:808727/141‑1 (MQ=255)
                                       agaTGGGCTTCGATGGTGTTGAGATTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGaa                                                                                                <  2:556222/143‑1 (MQ=255)
                                       agaTGGGCTTCGAGGGTTTGGGGATTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCACAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATCACGGATGGATTGGCGACTTTATCGaa                                                                                                <  2:1002098/143‑1 (MQ=255)
                                            ggCTTCGATGGTTTTGAGATTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTCTCGAAGAGcg                                                                                           >  2:980694/1‑143 (MQ=255)
                                              cTTCGATGGTTTTGAGATTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGcgtc                                                                                         >  1:237372/1‑143 (MQ=255)
                                               ttCGATGGTTTTGAGATTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGAGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGcgtcg                                                                                        >  1:79627/1‑143 (MQ=255)
                                                tCGATGGTTTTGAGATTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGagtcct                                                                                       <  2:430859/143‑3 (MQ=255)
                                                    tGGTTCTGAGATGGACGGCAAACCGCTGGTTACAAACATCGTAGATGGCGAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTa                                                                                   <  1:593344/143‑1 (MQ=255)
                                                     ggTTTTGAGATTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTaa                                                                                  <  1:199734/143‑1 (MQ=255)
                                                        tttGAGATTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGCCTTAATgg                                                                               >  2:446263/1‑143 (MQ=255)
                                                        tgagagaTGGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAACCCGGTTTACCGGTGACCGCCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATgg                                                                               <  1:265784/140‑1 (MQ=255)
                                                         tggaggTTGACGGCAAGCTGGTGGTTAACATCATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGCGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGc                                                                              <  2:527616/141‑1 (MQ=255)
                                                          tGAGATTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGCTTGGCGACTTTATCGAAGAGCGTCGTCTTCATGGCt                                                                             >  2:916373/1‑143 (MQ=255)
                                                                        aaCTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAg                                                               <  2:158375/143‑1 (MQ=255)
                                                                          ctgctgGTTAACAACATCGAAGAAGTCAACGCGGCGATCAAAGAAAACGGTTTACCGGTGACCCCCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAgcg                                                             <  1:538762/143‑1 (MQ=255)
                                                                                 ttAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTAGTCTTAATGGCTTAAAGCAGATCGAGCGCATtctc                                                      >  1:349218/1‑143 (MQ=255)
                                                                                   aacaacATCGAAGAAGTCAAGGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGa                                                    <  2:656629/143‑1 (MQ=255)
                                                                                   aacaacATCGAAGAAGTCAAAGCGGCGAGCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGa                                                    <  1:388199/143‑1 (MQ=255)
                                                                                     caacaTCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTCAAAGCAGAGCGAGGGCATTCTCGAAg                                                  >  2:6344/1‑143 (MQ=255)
                                                                                          tCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGACCGAGCGCATTCTCGAAGCGCTg                                             >  2:199140/1‑143 (MQ=255)
                                                                                           cGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTgg                                            <  1:1003403/143‑1 (MQ=255)
                                                                                           cGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTgg                                            <  2:832000/143‑1 (MQ=255)
                                                                                            gaagaaGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGc                                           >  1:589678/1‑143 (MQ=255)
                                                                                             aagaagTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGGCGTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCa                                          >  1:1041961/1‑143 (MQ=255)
                                                                                                aagTCAAAGCGGCTATCAAAGATATCTGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGaa                                       <  2:334673/143‑1 (MQ=255)
                                                                                                   tCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCCGATCGAGCGCATTCTCGAGGCGCTGGCAGAAGTg                                    >  1:14015/1‑143 (MQ=255)
                                                                                                     aaaGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTggg                                  >  1:610571/1‑143 (MQ=255)
                                                                                                       aGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCg                                >  1:1016779/1‑143 (MQ=255)
                                                                                                        gcggcgATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCgg                               <  2:643668/143‑1 (MQ=255)
                                                                                                         cggcgATCAAAGAAACCGGTTTACCGGTGACCACCGGCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGt                              >  1:682052/1‑143 (MQ=255)
                                                                                                            cgATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTaaa                           <  1:585796/143‑1 (MQ=255)
                                                                                                                cAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTACTGGCTTAAGGCAGATCGAGCGCATTCTCGGAGCGCTGGCAGAAGTGGGCGGTAAAGGCa                       >  2:518256/1‑143 (MQ=255)
                                                                                                                  aagaagCCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCCGATCCAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCAtc                     <  2:845024/143‑1 (MQ=255)
                                                                                                                  aaGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGACGAGCGCCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGCAGCGGGCGGTAAAGGAAtc                     >  2:913570/1‑143 (MQ=255)
                                                                                                                   aGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCAGTCTCGAAGCGCTGGCAGATGTGGGCGGTAAAGGCAtcg                    >  2:250944/1‑143 (MQ=255)
                                                                                                                        ccGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGtt               <  1:318910/143‑1 (MQ=255)
                                                                                                                                    tGACCTCCGCCCTCGGTGGCTGGGTCGGATGGCTTGGCGACTTTATCGAAGAGCGTCGGCTTAGTTGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTgg   <  2:933066/143‑1 (MQ=255)
                                                                                                                                     gcccgccgccgGCGGGGGCGAGGACGGAGGGAGTGGCGCCGGTACCGAAGAGCGTGGTCTTAATGGCTTCAAGCAGAGCGAGCGCATTGTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTggg  <  2:382492/138‑1 (MQ=255)
                                                                                                                                       ccaccGCAGGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGaa                                       >  2:384082/1‑104 (MQ=255)
                                                                                                                                       ccaccGCAGGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGaa                                       <  1:384082/104‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                     
TTCCGGAAAACATTCTGGAGAAATTTCGTTATATCAAAGAGATGGGCTTCGATGGTTTTGAGATTGACGGCAAACTGCTGGTTAACAACATCGAAGAAGTCAAAGCGGCGATCAAAGAAACCGGTTTACCGGTGACCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGG  >  NC_000913/1376065‑1376340

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