Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,549,932 T→A 8.9% intergenic (+222/+21) fimH → / ← gntP minor component of type 1 fimbriae/fructuronate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,549,9320TA 8.9% 220.3 / 10.4 67intergenic (+222/+21)fimH/gntPminor component of type 1 fimbriae/fructuronate transporter
Reads supporting (aligned to +/- strand):  ref base T (32/29);  new base A (3/3);  total (35/32)
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.29e-01

CCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAA  >  NC_000913/4549788‑4550078
                                                                                                                                                |                                                                                                                                                  
ccAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGc                                                                                                                                                <  1:721709/149‑1 (MQ=255)
         aCGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGTACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCAATGGCAgg                                                                                                                                             >  2:1200836/1‑143 (MQ=255)
         aCGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCAATGGCAgg                                                                                                                                             <  1:1200836/143‑1 (MQ=255)
                 gCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGgt                                                                                                                                      <  2:513080/142‑1 (MQ=255)
                 gCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGgt                                                                                                                                      >  1:513080/1‑142 (MQ=255)
                   gTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGc                                                                                                                                                <  3:554651/130‑1 (MQ=255)
                   gTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGc                                                                                                                                                >  4:554651/1‑130 (MQ=255)
                                 aaaCAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCAATGGCAGGc                                                                                                                                            >  2:446845/1‑120 (MQ=255)
                                 aaaCAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCAATGGCAGGc                                                                                                                                            <  1:446845/120‑1 (MQ=255)
                                   aCAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTgg                                                                                                                                                 <  1:837102/113‑1 (MQ=255)
                                   aCAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTgg                                                                                                                                                 >  2:837102/1‑113 (MQ=255)
                                                   tGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGctc                                                                                                                                          <  4:1408919/104‑1 (MQ=255)
                                                   tGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGctc                                                                                                                                          >  3:1408919/1‑104 (MQ=255)
                                                     gTCTGTTCCCCTCACCCTTACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACGTTGTCATCAGTCTGAGCCTGCCAATGGCAGGctc                                                                                                                                          <  1:747896/102‑1 (MQ=255)
                                                     gTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCAATGGCAGGctc                                                                                                                                          >  2:747896/1‑102 (MQ=255)
                                                               cTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTgg                                                                                                                                       <  2:1160363/95‑1 (MQ=255)
                                                               cTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTgg                                                                                                                                       >  1:1160363/1‑95 (MQ=255)
                                                                       aaCCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCtaa                                                                                                               >  4:1454414/1‑111 (MQ=255)
                                                                       aaCCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCtaa                                                                                                               <  3:1454414/111‑1 (MQ=255)
                                                                             ctcCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGca                                                                                                      <  4:153620/114‑1 (MQ=255)
                                                                             ctcCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGca                                                                                                      >  3:153620/1‑114 (MQ=255)
                                                                                cccGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGgtgt                                                                                                                                    <  3:1394772/81‑1 (MQ=255)
                                                                                cccGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGgtgt                                                                                                                                    <  1:1180459/81‑1 (MQ=255)
                                                                                cccGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGgtgt                                                                                                                                    >  2:1180459/1‑81 (MQ=255)
                                                                                cccGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGgtgt                                                                                                                                    >  4:1394772/1‑81 (MQ=255)
                                                                                cccGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACaataat                                                                                               >  3:690832/1‑118 (MQ=255)
                                                                                cccGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACaataat                                                                                               <  4:690832/118‑1 (MQ=255)
                                                                                cccGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTc                                                                    <  3:793642/145‑1 (MQ=255)
                                                                                cccGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTc                                                                    >  4:793642/1‑145 (MQ=255)
                                                                                cccGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCcagc                                                                >  4:1221343/1‑149 (MQ=255)
                                                                                 ccGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGa                                                                           >  2:125525/1‑137 (MQ=255)
                                                                                 ccGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGa                                                                           <  1:125525/137‑1 (MQ=255)
                                                                                 ccGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCcagca                                                               >  3:306934/1‑149 (MQ=255)
                                                                                         gCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCtaat                                                                                                              <  4:326813/94‑1 (MQ=255)
                                                                                         gCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCtaat                                                                                                              >  3:326813/1‑94 (MQ=255)
                                                                                                                   tAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAg                             <  3:1221343/149‑1 (MQ=255)
                                                                                                                        tttGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGca                                                                                                      <  3:1027166/71‑1 (MQ=255)
                                                                                                                        tttGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGca                                                                                                      >  4:1027166/1‑71 (MQ=255)
                                                                                                                        tttGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGa                                                             <  1:82090/112‑1 (MQ=255)
                                                                                                                        tttGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGa                                                             >  2:82090/1‑112 (MQ=255)
                                                                                                                        tttGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCaa                        >  3:1041460/1‑149 (MQ=255)
                                                                                                                        tttGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCaa                        >  1:52548/1‑149 (MQ=255)
                                                                                                                          tGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAg                      <  4:1041460/149‑1 (MQ=255)
                                                                                                                             catcaGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTaa                   <  2:52548/149‑1 (MQ=255)
                                                                                                                              atcaGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATc                                                                                                          >  3:1220062/1‑61 (MQ=255)
                                                                                                                              atcaGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATc                                                                                                          <  4:1220062/61‑1 (MQ=255)
                                                                                                                              atcaGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACaataat                                                                                               >  1:381287/1‑72 (MQ=255)
                                                                                                                              atcaGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACaataat                                                                                               <  2:381287/72‑1 (MQ=255)
                                                                                                                              atcaGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAAc                  >  1:1412853/1‑149 (MQ=255)
                                                                                                                                caGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAg                                                                                            >  2:297373/1‑73 (MQ=255)
                                                                                                                                caGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAg                                                                                            <  1:297373/73‑1 (MQ=255)
                                                                                                                                   tCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGAc                                                                          >  2:1463140/1‑88 (MQ=255)
                                                                                                                                   tCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGAc                                                                          <  1:1463140/88‑1 (MQ=255)
                                                                                                                                         ccTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGaacaac       <  4:306934/149‑1 (MQ=255)
                                                                                                                                           tGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTa                                    >  2:1128121/1‑118 (MQ=255)
                                                                                                                                           tGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTa                                    <  1:1128121/118‑1 (MQ=255)
                                                                                                                                            gCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAg                                                                                                        <  4:839104/49‑1 (MQ=255)
                                                                                                                                            gCCATTGGCAGCCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAg                                                                                                        >  3:839104/1‑49 (MQ=255)
                                                                                                                                            gCCATAGGCAGGCTCTGGTGTCCTTTTCCGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGGGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGGCAGGTCAAAGTAACCTTTGAACAACCAg    >  4:698545/1‑149 (MQ=255)
                                                                                                                                             ccATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACgt                                            >  2:1176006/1‑108 (MQ=255)
                                                                                                                                             ccATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACgt                                            <  1:1176006/108‑1 (MQ=255)
                                                                                                                                             ccATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAg                      <  2:1122069/130‑1 (MQ=255)
                                                                                                                                             ccATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAg                      >  1:1122069/1‑130 (MQ=255)
                                                                                                                                             ccATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGa   >  2:22659/1‑149 (MQ=255)
                                                                                                                                              cATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTGCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGaa  >  3:1209115/1‑149 (MQ=255)
                                                                                                                                               aTTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAAc                                              >  2:581377/1‑104 (MQ=255)
                                                                                                                                               aTTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAAc                                              <  1:581377/104‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
CCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAA  >  NC_000913/4549788‑4550078

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