Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 778,021 G→C 7.7% L94L (CTG→CTC)  tolB → periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913778,0210GC 7.7% 132.9 / 3.2 52L94L (CTG→CTC) tolBperiplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (23/25);  new base C (2/2);  total (25/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 = 9.23e-01

CGGGGCCTGGTGCGGCACCTGAAGATATTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGTTGTCGGTCAGGTCACTCCGAATCCGGATGGTTCTTACAATGTTGCTTATCAACTTGTTGACACTGGCGGCGCACCGGGTACTGTACTTGCTCAGAACTCGTACAAAGTGAACAAGCAGTGGCTGCG  >  NC_000913/777873‑778164
                                                                                                                                                    |                                                                                                                                               
cGGGGCCTGGTGCGGCACCTGAAGATATTGGCGGCATCGTTGCTTCTGACTTGCGTAACAGCGGTTAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTg                                                                                                                                                 <  7:393636/149‑1 (MQ=255)
    gCCTGGTGCGGCACCTGAAGATATTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCa                                                                                                                                             >  8:113867/1‑149 (MQ=255)
      cTGGTGCGGCACCTGAAGATATTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCAtt                                                                                                                                           <  8:134992/149‑1 (MQ=255)
                aCCTGAAGATATTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAg                                                                                                                                 >  3:282710/1‑149 (MQ=255)
                  cTGAAGATATTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGATAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGtt                                                                                                                               >  1:260071/1‑149 (MQ=255)
                   tGAAGATATTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCCGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCCCTCGGCATTGACGCTGTAGTTg                                                                                                                              >  1:242131/1‑149 (MQ=255)
                   tGAAGATATTGGCGGCATAGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTACGCACTCGGCATTGACGCTGTAGTTg                                                                                                                              >  1:247781/1‑149 (MQ=255)
                    gAAGATATTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGTTGt                                                                                                                             >  6:188274/1‑149 (MQ=255)
                     aaGATATTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGATAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGTTGt                                                                                                                             <  2:260071/148‑1 (MQ=255)
                      atatatTGGCGGCATCTTTGCTGCTTACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGACGGGTAGTGCGCAGGAAGTACAACCAGCTGAATGGTCCGCACTCGGCATTGACGCTGTAGTTGTCg                                                                                                                           <  2:247781/147‑1 (MQ=255)
                      aGATATTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTCGGCATTGACGCTGTAGTTGTCg                                                                                                                           <  2:242131/149‑1 (MQ=255)
                          atTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGTTGTCggt                                                                                                                         >  2:239145/1‑147 (MQ=255)
                           ttGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGt                                                                                                                                <  7:69223/139‑1 (MQ=255)
                           ttGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGt                                                                                                                                >  8:69223/1‑139 (MQ=255)
                                 gcATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGTTGTCGGTCAGGTCACt                                                                                                                >  4:100681/1‑149 (MQ=255)
                                  cATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGTTGTCGGTCAGGTCACTc                                                                                                               <  1:10718/149‑1 (MQ=255)
                                    tCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAg                                                                                                                                 <  4:242665/129‑1 (MQ=255)
                                    tCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAg                                                                                                                                 >  3:242665/1‑129 (MQ=255)
                                       ttGCTGCTGCCTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGTTGTCGGTCAGGTCACTCCGAAt                                                                                                          <  4:251733/149‑1 (MQ=255)
                                        tgctgcTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCTCTCGTCTGCCACATCAGACGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGTTGTCGGTCAGGTCACTCCGAATc                                                                                                         >  7:84547/1‑149 (MQ=255)
                                          ctgctgACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGTTGTCGGTCAGGTCACTCCGAATCCg                                                                                                       <  1:215207/149‑1 (MQ=255)
                                           tgctgACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGTTGTCGGTCAGGTCACTCCGAATCCg                                                                                                       >  5:132981/1‑148 (MQ=255)
                                             ctgACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGTTGTCGGTCAGGTCACTCCGAATCCGGa                                                                                                     <  2:216288/148‑1 (MQ=255)
                                                aCTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACACCCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGTTGTCGGTCAGGTCACTCCGAATCCGGATGGt                                                                                                 >  6:350516/1‑149 (MQ=255)
                                                         aCAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGTTGTCGGTCAGGTCACTCCGAATCCGGATGGTTCTTACAAt                                                                                        <  5:243865/149‑1 (MQ=255)
                                                          cAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGTTGTCGGTCAGGTCACTCCGAATCCGGATGGTTCTTACAATg                                                                                       >  6:330117/1‑149 (MQ=255)
                                                           aGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGTTGTCGGTCAGGTCACTCCGAATCCGGATGGTTCTTACAATga                                                                                      >  5:196906/1‑148 (MQ=255)
                                                           aGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGTTGTCGGTCAGGTCACTCCGAATCCGGATGGTTCTTACAATGt                                                                                      >  3:354626/1‑149 (MQ=255)
                                                             cGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGTTGTCGGTCAGGTCACTCCGAATCCGGATGGTTCTTACAATGTTg                                                                                    <  1:239145/149‑1 (MQ=255)
                                                               gTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACCACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGTTGTCGGTCAGGTCACTCAGAATCCGGATGGTTCTTACAATGTTGCt                                                                                  >  5:341658/1‑149 (MQ=255)
                                                                   aTTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGTTGTCGGTCAGGTCACTCCGAATCCGGATGGTTCTTACAATGTTGCTTATc                                                                              <  3:100681/149‑1 (MQ=255)
                                                                         tCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGTTGTCGGTCAGGTCACTCCGAATCCGGATGGTTCTTACAATGTTGCTTATCAACttg                                                                        >  1:296397/1‑149 (MQ=255)
                                                                           cGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGTTGTCGGTCAGGTCACTCCGAATCCGGATGGTTCTTACAATGTTGCTTATCAACttgtt                                                                      <  3:265518/149‑1 (MQ=255)
                                                                           cGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGTTGTCGGTCAGGTCACTCCGAATCCGGATGGTTCTTACAATGTTGCTTATCAACttgt                                                                       <  6:132981/148‑1 (MQ=255)
                                                                            gTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGTTGTCGGTCAGGTCACTCCGAATCCGGATGGTTCTTACAATGTTGCTTATCAACttgttg                                                                     <  7:247263/149‑1 (MQ=255)
                                                                                gATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGTTGTCGGTCAGGTCACTCCGAATCCGGATGGTTCTTACAATGTTGCTTATCAACTTGTTGacac                                                                 >  5:23537/1‑149 (MQ=255)
                                                                                 aTCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGTTGTCGGTCAGGTCACTCCGAATCCGGATGGTTCTTACAATGTTGCTTATCAACTTGTTGACACt                                                                >  8:270651/1‑149 (MQ=255)
                                                                                   cgcgCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGTTGTCGGTCAGGTCACTCCGAATCCGGATGGTTCTTACAATGTTGCTTATCAACTTGTTGACACTgg                                                              >  2:341870/1‑149 (MQ=255)
                                                                                               ccACAGCAGCCGGGTAGTGCGCAGGAAGTACACCCAGCTGCATGGTCCGCACTGGGCATTGACGCGGTAGTTGTCGGTCAGGTCACTCCGAATCCGGATGGTTCTTACAATGTTGCTTATCAACTTGTTGACACTGGCGGCGCACCg                                                    <  6:379391/147‑1 (MQ=255)
                                                                                               ccACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACAGGGCATTGCCGCTGTAGTTGTCGGTCAGCTCACTCCGAATCCGGATGGTTCTTACAATGTTGCTTATCAACTTTTTGACACTGGCGGCGCACCggg                                                  >  1:268937/1‑149 (MQ=255)
                                                                                                cacaGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGTTGTCGGTCAGGTCACTCCGAATCCGGATGGTTCTTACAATGTTGCTTATCAACTTGTTGACACTGGCGGCGCACCggg                                                  >  6:375040/1‑148 (MQ=255)
                                                                                                  cagcagCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGTTGTCGGTCAGGTCACTCCGAATCCGGATGGTTCTTACAATGTTGCTTATCAACTTGTTGACACTGGCGGCGCACCGGgtac                                               <  6:23537/149‑1 (MQ=255)
                                                                                                        ccGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGTTGTCGGTCAGGTCACTCCGAATCCGGATGGTTCTTACAATGTTGCTTATCAACTTGTTGACACTGGCGGCGCACCGGgtactgtact                                         <  6:196906/149‑1 (MQ=255)
                                                                                                                    cAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGTTGTCGGTCAGGTCACTCCGAATCCGGATGGTTCTTACAATGTTGCTTATCAACTTGTTGACACTGGCGGCGCACCGGGTACTGTACTTGCTCAGAACTc                             <  3:76527/149‑1 (MQ=255)
                                                                                                                      ggAAGTACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGTTGTCGGTCAGGTCACTCCGAATCCGGATGGTTCTTACAATGTTGCTTATCAACTTGTTGACACTGGCGGCGCACCGGGTACTGTACTTGCTCAGAACTc                             >  7:281635/1‑147 (MQ=255)
                                                                                                                      ggAAGTACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGTTGTCGGTCAGGTCACTCCGAATCCCGATGGTTCTTACAATGTTGCTTATCAACTTGTTGACACTGGCGGCGCACCGGGTACTGTACTTGCTCAGAACTc                             <  8:281635/147‑1 (MQ=255)
                                                                                                                           tACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGTTGTCGGTCAGGTCACTCCGAATCCGGATGGTTCTTACAATGTTGCTTATCAACTTGTTGACACTGGCGGCGCACCGGGTACTGTACTTGCTCAGAACTCGTACaaa                      <  5:330117/149‑1 (MQ=255)
                                                                                                                              aaCCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGTTGTCGGTCAGGTCACTCCGAATCCGGATGGTTCTTACAATGTTGCTTATCAACTTGTTGACACTGGCGGCGCACCGGGTACTGTACTTGCTCAGAACTCGTACAAAGTg                   <  4:325796/149‑1 (MQ=255)
                                                                                                                              aaCCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGTTGTCGGTCAGGTCACTCCGAATCCGGATGGTTCTTACAATGTTGCTTATCAACTTGTTGACACTGGCGGCGCACCGGGTACTGTACTTGCTCAGAACTCGTACAAAGTg                   <  4:282710/149‑1 (MQ=255)
                                                                                                                               tccAGCTGCTTGGTCCGCACTGGGCATTGACTCTGTAGTTGTCGGTCAGGTCACTCCGAATCCGGATGTTTCTTACAATGTTGCTTATCAACTTGTTGACACTGGCGGCGCACCGTGTACTGTACTTGCTCAGAACTCTTACAAAGtta                  <  7:198363/148‑3 (MQ=255)
                                                                                                                                    cTGCATGGTCCGCACTGGGCATTGACGCTGTAGTTGTCGGTCAGGTCACTCCGAATCCGGATGGTTCTTACAATGTTGCTTATCAACTTGTTGACACTGGCGGCGCACCGGGTACTGTACTTGCTCAGAACTCGTACAAAGTGAACAAg             <  5:188274/149‑1 (MQ=255)
                                                                                                                                       cATGGTCCGCACTGGGCATTGACGCTGTAGTTGTCGGTCAGGTCACTCCGAATCCGGATGGTTCTTACAATGTTGCTTATCAACTTGTTGACACTGGCGGCGCACCGGGTCCTGTACTTGCTCAGAACTCGTACAAAGTGAACAAGCAg          >  1:163819/1‑149 (MQ=255)
                                                                                                                                               gCACTGGGCATTGACGCTGTAGTTGTCGGTCAGGTCACTCCGAATCCGGATGGTTCTTACAATGTTGCTTATCAACTTGTTGACACTGGCGGCGCACCGGGTACTGTACTTGCTCAGAACTCGTACAAAGTGAACAAGCAGTGGCTGCg  <  1:341870/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CGGGGCCTGGTGCGGCACCTGAAGATATTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCATTGACGCTGTAGTTGTCGGTCAGGTCACTCCGAATCCGGATGGTTCTTACAATGTTGCTTATCAACTTGTTGACACTGGCGGCGCACCGGGTACTGTACTTGCTCAGAACTCGTACAAAGTGAACAAGCAGTGGCTGCG  >  NC_000913/777873‑778164

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