Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,326,373 T→C 7.0% intergenic (‑632/+26) yjdN ← / ← yjdM metalloprotein superfamily protein/zinc‑ribbon family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,326,3730TC 7.0% 192.2 / 3.6 57intergenic (‑632/+26)yjdN/yjdMmetalloprotein superfamily protein/zinc‑ribbon family protein
Reads supporting (aligned to +/- strand):  ref base T (28/25);  new base C (2/2);  total (30/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.76e-01

TCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATCGATGTTATGGTCGCCTTCAACCAGGCGGATGTTTTTCACTTTGGTGCCAATTTTCAGCA  >  NC_000913/4326240‑4326520
                                                                                                                                     |                                                                                                                                                   
tCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGAt                                                                                                                                      >  2:1194298/1‑149 (MQ=255)
                 tGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTc                                                                                                                     <  1:1194298/149‑1 (MQ=255)
                 tGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTc                                                                                                                     >  4:1274195/1‑149 (MQ=255)
                    ggcgTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCacaatacaa                                                                                                                            <  3:182421/139‑1 (MQ=255)
                    ggcgTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCacaatacaa                                                                                                                            >  4:182421/1‑139 (MQ=255)
                           aCGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCacaatacaat                                                                                                                           <  2:1861482/133‑1 (MQ=255)
                           aCGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCacaatacaat                                                                                                                           >  1:1861482/1‑133 (MQ=255)
                           aCGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCacaat                                                                                                                                >  1:1162290/1‑128 (MQ=255)
                           aCGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCacaat                                                                                                                                <  2:1162290/128‑1 (MQ=255)
                           aCGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAg                                                                                                                        >  2:1183592/1‑136 (MQ=255)
                           aCGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAg                                                                                                                        <  1:1183592/136‑1 (MQ=255)
                                ttATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCa                                                                                                      <  3:1274195/149‑1 (MQ=255)
                                     cGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGAttt                                                                                                 >  3:983005/1‑149 (MQ=255)
                                       gCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCa                                                                                               <  4:983005/149‑1 (MQ=255)
                                                 ggTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCgg                                                                                     >  1:1387304/1‑149 (MQ=255)
                                                        gAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCaaaa                                                                              >  3:1362695/1‑149 (MQ=255)
                                                                  gCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCacaa                                                                                                                                 <  3:858558/88‑1 (MQ=255)
                                                                  gCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCacaa                                                                                                                                 >  4:858558/1‑88 (MQ=255)
                                                                         aaGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGAt                                                                                                                                      >  1:675754/1‑76 (MQ=255)
                                                                         aaGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGAt                                                                                                                                      <  2:675754/76‑1 (MQ=255)
                                                                           gACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCtt                                                                                                                   >  4:1091945/1‑93 (MQ=255)
                                                                           gACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCtt                                                                                                                   <  3:1091945/93‑1 (MQ=255)
                                                                           gACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTc                                                                         <  1:1401944/135‑1 (MQ=255)
                                                                           gACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTc                                                                         >  2:1401944/1‑135 (MQ=255)
                                                                           gACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCaa                                                               <  4:148469/145‑1 (MQ=255)
                                                                           gACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCaa                                                               >  3:148469/1‑145 (MQ=255)
                                                                            aCGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATAAGGCGCTCGCGCCGCATcc                                                                                                                                                 >  2:394902/1‑62 (MQ=17)
                                                                            aCGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATAAGGCGCTCGCGCCGCATcc                                                                                                                                                 <  1:394902/62‑1 (MQ=17)
                                                                              gcgcAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATCGATGt                                                        >  1:1947/1‑149 (MQ=255)
                                                                                gcAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAAc                                                                                                         <  3:256378/98‑1 (MQ=255)
                                                                                gcAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAAc                                                                                                         >  4:256378/1‑98 (MQ=255)
                                                                                gcAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATCGATGTTa                                                      >  4:999700/1‑149 (MQ=255)
                                                                                gcAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATCGATGTTa                                                      <  3:999700/149‑1 (MQ=255)
                                                                                 cAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATCGATGTTAt                                                     <  2:1387304/149‑1 (MQ=255)
                                                                                               ggCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCaa                                                               >  3:1335875/1‑125 (MQ=255)
                                                                                               ggCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCaa                                                               <  4:1335875/125‑1 (MQ=255)
                                                                                                gCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTc                                                                                                       >  4:1593045/1‑84 (MQ=255)
                                                                                                gCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTc                                                                                                       <  3:1593045/84‑1 (MQ=255)
                                                                                                  aTCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATCGATGTTATGGTCGCCTTCaa                                         >  4:1507162/1‑144 (MQ=255)
                                                                                                  aTCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATCGATGTTATGGTCGCCTTCaa                                         <  3:1507162/144‑1 (MQ=255)
                                                                                                   tCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATCGATGTTATGGTCGCCTTCAACCAGGc                                   <  4:1362695/149‑1 (MQ=255)
                                                                                                           cGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATCGATGTTATGGTCGCCTTCAACCAGGCGGATGttt                           <  2:1947/149‑1 (MQ=255)
                                                                                                              agCCGGATAAGGCGCTCGCGCCGCATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATCGATGTTATGGTCGcc                                              <  2:470163/126‑1 (MQ=39)
                                                                                                              agCCGGATAAGGCGCTCGCGCCGCATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATCGATGTTATGGTCGcc                                              >  1:470163/2‑127 (MQ=38)
                                                                                                                  ggATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATCGATGTTATGGTCGCCTTCAACCAGGCGGa                                <  2:660012/137‑1 (MQ=255)
                                                                                                                  ggATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATCGATGTTATGGTCGCCTTCAACCAGGCGGa                                >  1:660012/1‑137 (MQ=255)
                                                                                                                  ggATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATCGATGTTATGGTCGCCTTCAACCAGGCGGATGTTTTTCACtt                    >  3:1248311/1‑149 (MQ=255)
                                                                                                                        ggcgCGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATCGATGTTATGGTCGCCTTCAACCAGGCGGATGTTTTTCACTTTGGTGc              >  3:1411645/1‑149 (MQ=255)
                                                                                                                               gcgcCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATc                                                             >  2:938933/1‑95 (MQ=255)
                                                                                                                               gcgcCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATc                                                             <  1:938933/95‑1 (MQ=255)
                                                                                                                                cgcCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATCGATGTTATGGTCGCCTTCAACCAGGCGGATGTTTTTCACTTTGGTGCCaa           <  2:218628/144‑1 (MQ=255)
                                                                                                                                cgcCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATCGATGTTATGGTCGCCTTCAACCAGGCGGATGTTTTTCACTTTGGTGCCaa           >  1:218628/1‑144 (MQ=255)
                                                                                                                                 gcCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATCGATGTTATGGTCGCCTTCAACCAGGCGGATGTTTTTCACTTTGGTGCCAATTTTCa     >  3:1455445/1‑149 (MQ=255)
                                                                                                                                   cTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAATCCGTCGATTTTGCAATCGATGTTATGGTCGCCTTCaa                                         <  4:697960/111‑1 (MQ=255)
                                                                                                                                   cTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATCGATGTTATGGTCGCCTTCaa                                         >  3:697960/1‑111 (MQ=255)
                                                                                                                                   cTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATCGATGTTATGGTCGCCTTCAACCAGGCGGATGTTTTTCa                       <  2:1405129/129‑1 (MQ=255)
                                                                                                                                   cTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATCGATGTTATGGTCGCCTTCAACCAGGCGGATGTTTTTCa                       >  1:1405129/1‑129 (MQ=255)
                                                                                                                                    ttATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATCGATGTTATgg                                                   <  3:1373355/100‑1 (MQ=255)
                                                                                                                                    ttATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATCGATGTTATgg                                                   >  4:1373355/1‑100 (MQ=255)
                                                                                                                                    ttATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATCGATGTTATGGTCGCCTTCAACCAGGCGGATGTTTTTCACTTTGGTGCCAATTTTCAGCa  >  1:1074204/1‑149 (MQ=255)
                                                                                                                                     |                                                                                                                                                   
TCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATCTGGCACCGTTGCCGGATGCGGCGCGAGCGCCTTATCCGGCTTACCGATCACAATACAATCAGTTCTTTTTCACAAACTCAGATTTCAGTTTCATCGGACCAAAACCGTCGATTTTGCAATCGATGTTATGGTCGCCTTCAACCAGGCGGATGTTTTTCACTTTGGTGCCAATTTTCAGCA  >  NC_000913/4326240‑4326520

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