Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,951,386 2 bp→TC coding (555‑556/567 nt) yecD → isochorismatase family protein YecD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,951,3860CT89.5% 84.2 / ‑0.4 40L185L (CTC→CTT) yecDisochorismatase family protein YecD
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base T (12/22);  total (13/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009131,951,3870AC89.5% 100.1 / 3.3 38N186H (AAC→CAC) yecDisochorismatase family protein YecD
Reads supporting (aligned to +/- strand):  ref base A (1/3);  new base C (12/22);  total (13/25)
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.46e-01

GGAGTTGCAATTACGCCGCCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑CAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAACC  >  NC_000913/1951167‑1951525
                                                                                                                                                                                                                            ||                                                                                                                                           
cgatgtgcggttACGCCGCCGAGGTATCGATACAAGTGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑TCACGCGTTATTATGTACATCGGTCTGCCGc                                                                                                                <  1:133693/240‑1 (MQ=255)
       cggttACGCCGCCGAGTTATCGATACAATTGTGTTATGTGGGATCTCGACCAATAGCGGTGGTGAATCCTCCGCCCGCAATGCCTGGGAACTGGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAAGAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCTTAGCGGGGAAGAGATCCT‑TCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCg                                                                                                        <  1:65340/248‑1 (MQ=255)
         gttACGCCGCCGAGGTATCGATACAATTGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑TCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGc                                                                                                       >  1:82473/2‑250 (MQ=255)
             cgccgactaGGTATCGATACAATTGTGGTATGTGGGATCTCGACCAATATCGGTGTTGAATCCCCCGCCCGCAATGCCAGGGATCGTGGCTTTAAGCTGGTGATGGCCGAAGACGCCATTAGTGCCTCCCGTGCCGAGCAACACAATAACAGCATGAATCATAACGAGCCGCGCATCGACCGTGTGCGTAGCGTGGAAGAGATCCT‑TCACGCGTTATGATTTACATc                                                                                                                          <  1:176469/219‑1 (MQ=255)
                    gAGGTATCGATACAAGTGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGGGGAAGAGATCCT‑TCACGCGTTATGATTTACAGCGTTCTGCCGCAATGGTCGCAT‑CCTAAATggn                                                                                          <  1:18911/251‑2 (MQ=255)
                      ggTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑CAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCAT‑CCTAAATGGGTg                                                                                        <  1:138799/251‑1 (MQ=255)
                              cagactgtgctcttccgatctctcGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑TCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTggg                                                                                <  1:79384/231‑1 (MQ=255)
                               tcAATTGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCGACCCGCTCATCGCCCGTGTGCGTAGCGGGGAAGAGATCCT‑TCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTggg                                                                                <  1:132636/249‑1 (MQ=255)
                               aCAATTGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTTGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑TCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTggg                                                                                <  1:73708/250‑1 (MQ=255)
                               aCAATTGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑TCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTggg                                                                                <  1:7045/250‑1 (MQ=255)
                                      tgtTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTTGAAGAGATCCT‑TCACGCGGTATGATTTACATCGGTCTGCCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTGGGTATCAc                                                                          <  1:101796/249‑1 (MQ=255)
                                                 aTCTCGACCAATATCGGTGTTGAATCCACCGCCCGCACTGCCTTGTAACTTGTCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCATCACAATAACAGCATTAATAATAGCTACCCGCGCAGCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑TCACGCGTTAt                                                                                                                                    <  1:168708/181‑1 (MQ=255)
                                                  tctcGACCAATATCGGTGTTGAATCCNCCGCCCNNAANGCCTGGGAACTTGGCTTTANTCNNNTNNNTNCCGAANACNNCTGTAGTGCNGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCNCGTGTGCGTNNCGTGGAAGAGATCCT‑TCACGCNTTANNANNTNCATNgg                                                                                                                        >  2:40008/1‑192 (MQ=255)
                                                     cGACCAATCTCGGTGGTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTTGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑GCACGAGTTATGATTTACATCGGTCTGCCTCAATGGTCGCAT‑CCTAAATTGGTTCGGGTGGGTATCACCATCCTTGAAGAGTAt                                                          <  1:49294/250‑1 (MQ=255)
                                                     cGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑TCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATg                                                         >  1:28264/1‑251 (MQ=255)
                                                      gACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑TCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATGc                                                        >  1:26449/1‑251 (MQ=255)
                                                       aCCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTTCCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑GCACGCGTTATGATTTACATCTGTCTGCCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATGcc                                                       <  1:69877/251‑1 (MQ=255)
                                                        ccAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCACCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑TCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTggg                                                                                >  1:53542/1‑225 (MQ=255)
                                                        ccAAGATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑TCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTGGGTATCAc                                                                          <  1:72487/231‑1 (MQ=255)
                                                        ccAAGATCGGTGTTGAAGCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGA‑GACGACTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGGGCGTAGCGTGGAAGAGATCCT‑TCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATGccct                                                     <  1:15784/251‑2 (MQ=255)
                                                                     tGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCTAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑TCACGCGTTATGATTTACATCGGTCTTCCGCAATGGTCg                                                                                                        <  1:25610/189‑1 (MQ=255)
                                                                               gcccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGGGCGTAGCGTGGAATAGATCCT‑TCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTGGGTATCAcc                                                                         <  1:125776/209‑1 (MQ=255)
                                                                                   gcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑TCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCAT‑CCTAAATGGGTg                                                                                        <  1:55780/190‑1 (MQ=255)
                                                                                    cAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑TCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGAGTgacgcgggcatt                          >  1:48018/1‑240 (MQ=255)
                                                                                         ccTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑TCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGCAGAGTATGCCCGCCACTTTAACTGCGTgacgcgggcattttaa                      >  1:151208/1‑235 (MQ=255)
                                                                                                 cTTGGATTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑TCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTgacgcgggcattttaaaaatcact              >  1:4047/1‑227 (MQ=255)
                                                                                                     gCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑TCACGCGTTATGATTTACATCGGTc                                                                                                                      <  1:134465/143‑1 (MQ=255)
                                                                                                        ttAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAAGCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑TCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTGGGTATCACCAGCCTTg                                                                  <  1:142672/191‑1 (MQ=255)
                                                                                                             tggtggATTGCCAAAGACGCCGGTAGGCCCGCCAGCGCCGAGCAGCACAAAAACAGTAGCAACCATACCCCCCCGCGCATCGCCGGTGTGCGTAGCGGGGAAGATACCCT‑CACCGCGTTATGTTTTACACCGTTCTCCCGCAAGGGTC‑CATTCCTAAAGGGGTGGGTTTGGGGACCACCAGCCTTGAAGAGTATCCCCGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAAcc  <  1:4390/245‑1 (MQ=255)
                                                                                                                                ccTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑CAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAag                                                               <  1:139612/170‑1 (MQ=255)
                                                                                                                                 cTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑TCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCAT‑CCTAAATGGGTGCGGGTGTGTATCACCAGCCTTGAAGAGTATGct                                                       <  1:51332/177‑2 (MQ=255)
                                                                                                                                           gccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑TCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTg                                                                                  <  1:9408/140‑1 (MQ=255)
                                                                                                                                           gccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑TCACGCGTNNTGAttt                                                                                                                               >  2:9408/1‑96 (MQ=255)
                                                                                                                                             cAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑TCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCAT‑CCTAAATggg                                                                                          >  1:167136/1‑130 (MQ=255)
                                                                                                                                                       naacacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑TCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAagcg                                                             >  1:19354/4‑147 (MQ=255)
                                                                                                                                                                          atcATATCTCCCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCAA‑GCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTGGGTATCACCa                                                                        >  1:40797/3‑119 (MQ=255)
                                                                                                                                                                          atcATATCTACCCGCGCATCGCCCGTGTGCNTAGCGTGGANNAGNTCCT‑TNANNNGNNNTGNTTTACATCGGTCTGCCGCANTGNNCGCAT‑CNTAAATGGGTGCGGTTGGGTATCACCa                                                                        <  2:40797/117‑1 (MQ=255)
                                                                                                                                                                                        cgcATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑TCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGAAT‑CATAAATGGGTGCGGTTGGGTATCACCAg                                                                       >  1:128772/1‑106 (MQ=255)
                                                                                                                                                                                        cgcATCGCCCGNGTGCGTAGCGNGGANNAGATNCN‑TNACGCGTTATGATTTACATCGGTCTGCCGCNATGGTCGCAT‑CCTAAATGGGTGCNNTTNNNTANNACCAg                                                                       <  2:128772/106‑1 (MQ=255)
                                                                                                                                                                                                    gtgCGTAGCGGGGAAGAGGTCCT‑TCACGCGTTATGATTTGCATCGGTCTGCCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTGGGTATc                                                                            <  1:160080/89‑1 (MQ=255)
                                                                                                                                                                                                                            ||                                                                                                                                           
GGAGTTGCAATTACGCCGCCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCT‑CAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCAT‑CCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAACC  >  NC_000913/1951167‑1951525

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