Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,274,740 G→A 13.4% intergenic (+79/+20) ssb → / ← yjcB single‑stranded DNA‑binding protein/putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,274,7400GA13.4% 183.2 / 18.2 67intergenic (+79/+20)ssb/yjcBsingle‑stranded DNA‑binding protein/putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (31/27);  new base A (5/4);  total (36/31)
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.25e-01

GAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAAC  >  NC_000913/4274623‑4274887
                                                                                                                     |                                                                                                                                                   
gAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCATAACGGAGGc                                                                                                                                           >  4:1737099/1‑128 (MQ=255)
gAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCATAACGGAGGc                                                                                                                                           <  3:1737099/128‑1 (MQ=255)
                   atgatgaCATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATg                                                                                                                                                <  2:163208/104‑1 (MQ=255)
                   atgatgaCATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATg                                                                                                                                                <  4:1848319/104‑1 (MQ=255)
                   atgatgaCATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATg                                                                                                                                                >  1:163208/1‑104 (MQ=255)
                   atgatgaCATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATg                                                                                                                                                >  3:1848319/1‑104 (MQ=255)
                     gatgaCATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATgg                                                                                                                                               >  4:1357652/1‑103 (MQ=255)
                     gatgaCATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATgg                                                                                                                                               <  3:1357652/103‑1 (MQ=255)
                        gaCATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATg                                                                                                                                                >  4:199216/1‑99 (MQ=255)
                        gaCATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATg                                                                                                                                                <  3:199216/99‑1 (MQ=255)
                                gTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAgg                                                                                                                                            <  2:620251/95‑1 (MQ=255)
                                gTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAgg                                                                                                                                            >  1:620251/1‑95 (MQ=255)
                                     gATTTGTCATTAAAAAAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCtcttc                                                                                  <  3:321625/148‑1 (MQ=255)
                                     gATTTGTCATTAAAAAAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCAGCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCtcttc                                                                                  >  4:321625/1‑148 (MQ=255)
                                              ttAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCAtt                                                                                                                  >  1:634979/1‑107 (MQ=255)
                                              ttAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCAtt                                                                                                                  <  2:634979/107‑1 (MQ=255)
                                              ttAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCATAACGGAGGCtt                                                                                                                                         <  2:671129/84‑1 (MQ=39)
                                              ttAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCATAACGGAGGCtt                                                                                                                                         >  1:671129/1‑84 (MQ=38)
                                                         ggTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGAt                                                                                            <  1:1271814/118‑1 (MQ=255)
                                                         ggTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGAt                                                                                            >  2:1271814/1‑118 (MQ=255)
                                                          gTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATc                                                                                                                                >  3:486980/1‑81 (MQ=255)
                                                          gTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATc                                                                                                                                <  4:486980/81‑1 (MQ=255)
                                                          gTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCATAACGGAGGc                                                                                                                                           <  1:1828884/70‑1 (MQ=38)
                                                          gTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCATAACGGAGGc                                                                                                                                           >  2:1828884/1‑70 (MQ=38)
                                                               atTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCa                                                       >  3:1027669/1‑149 (MQ=255)
                                                                   ttttAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGc                                                                                                      <  3:1383112/98‑1 (MQ=255)
                                                                   ttttAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGc                                                                                                      >  4:1383112/1‑98 (MQ=255)
                                                                   ttttAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGgt                                                                                                   >  1:722454/1‑101 (MQ=255)
                                                                   ttttAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGgt                                                                                                   <  2:722454/101‑1 (MQ=255)
                                                                        aGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGt                                              <  4:1027669/149‑1 (MQ=255)
                                                                                            tGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTatcatc                                                                           >  3:1869150/1‑100 (MQ=255)
                                                                                            tGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTatcatc                                                                           <  4:1869150/100‑1 (MQ=255)
                                                                                              ccagccaTAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGa                                            <  2:311736/129‑1 (MQ=255)
                                                                                              ccagccaTAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGa                                            >  1:311736/1‑129 (MQ=255)
                                                                                                 gccaTAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGa                                                                                             >  4:1680794/1‑77 (MQ=255)
                                                                                                 gccaTAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGa                                                                                             <  3:1680794/77‑1 (MQ=255)
                                                                                                     tAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCtcttc                                                                                  <  2:613398/84‑1 (MQ=255)
                                                                                                     tAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCtcttc                                                                                  >  1:613398/1‑84 (MQ=255)
                                                                                                      aaaaaaGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCAtt                                                                                                           <  4:1724955/58‑1 (MQ=255)
                                                                                                      aaaaaaGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCAtt                                                                                                           >  3:1724955/1‑58 (MQ=255)
                                                                                                        aaaaGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGa                             >  3:1880912/1‑134 (MQ=255)
                                                                                                        aaaaGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGa                             <  4:1880912/134‑1 (MQ=255)
                                                                                                            gaagCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGAtgttg                                                            >  4:301421/1‑99 (MQ=255)
                                                                                                            gaagCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGAtgttg                                                            <  3:301421/99‑1 (MQ=255)
                                                                                                             aagCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTcc                                       <  2:754711/119‑1 (MQ=255)
                                                                                                             aagCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTcc                                       >  1:754711/1‑119 (MQ=255)
                                                                                                               gCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCaa                       >  1:1336734/1‑133 (MQ=255)
                                                                                                               gCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCaa                       <  2:1336734/133‑1 (MQ=255)
                                                                                                               gCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGaa       >  3:969353/1‑149 (MQ=255)
                                                                                                               gCCTCCATAACGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGc                     <  1:1068411/135‑1 (MQ=255)
                                                                                                               gCCTCCATAACGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGc                     >  2:1068411/1‑135 (MQ=255)
                                                                                                               gCCTCCATAACGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGaa       >  1:606645/1‑149 (MQ=255)
                                                                                                                ccTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGTTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAAc      >  1:1213425/1‑149 (MQ=255)
                                                                                                                 cTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAAccc                <  1:218681/138‑1 (MQ=255)
                                                                                                                 cTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAAccc                >  2:218681/1‑138 (MQ=255)
                                                                                                                  tCCGTTATGGAGGCTTCTAGGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTc                                        >  3:1568588/1‑113 (MQ=255)
                                                                                                                  tCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGAtgttgt                                                           <  2:1336204/94‑1 (MQ=255)
                                                                                                                  tCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGAtgttgt                                                           >  2:1411753/1‑94 (MQ=255)
                                                                                                                  tCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGAtgttgt                                                           <  1:1411753/94‑1 (MQ=255)
                                                                                                                  tCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGAtgttgt                                                           >  1:1336204/1‑94 (MQ=255)
                                                                                                                  tCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTc                                        <  4:1568588/113‑1 (MQ=255)
                                                                                                                  tCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGcc                         <  1:1628910/128‑1 (MQ=255)
                                                                                                                  tCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGcc                         >  2:1628910/1‑128 (MQ=255)
                                                                                                                  tCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAAccc                >  1:1422189/1‑137 (MQ=255)
                                                                                                                  tCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAAccc                <  2:1422189/137‑1 (MQ=255)
                                                                                                                  tCCGTTATAGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCa    >  3:1742207/1‑149 (MQ=255)
                                                                                                                    cGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAAc  >  1:660819/1‑149 (MQ=255)
                                                                                                                     gTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTat                                                                               <  3:1314271/71‑1 (MQ=255)
                                                                                                                     gTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTat                                                                               >  4:1314271/1‑71 (MQ=255)
                                                                                                                     |                                                                                                                                                   
GAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAAC  >  NC_000913/4274623‑4274887

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