Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,635,048 A→C F88C (TTT→TGT)  rob ← right oriC‑binding transcriptional activator, AraC family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,635,0480AC96.0% 166.3 / ‑6.1 50F88C (TTT→TGT) robright oriC‑binding transcriptional activator, AraC family
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (19/29);  total (20/30)
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.85e-01

GCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACATTCTCT  >  NC_000913/4634902‑4635195
                                                                                                                                                  |                                                                                                                                                   
gcgTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTACATGt                                                                                                                                                 <  2:964565/151‑1 (MQ=255)
    aTCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTACATGTCTg                                                                                                                                              <  1:611594/150‑1 (MQ=255)
        tCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTACATGTCTGt                                                                                                                                             <  2:715539/147‑1 (MQ=255)
        tCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTACATGTCTGTTg                                                                                                                                           <  1:297120/149‑1 (MQ=255)
         ccAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTACATGTCTGTTgaga                                                                                                                                        <  1:320478/151‑1 (MQ=255)
           aGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTACATGTCTGTTGAGAGt                                                                                                                                      <  2:909874/151‑1 (MQ=255)
                       ttGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCCCCCTTCAGGAGAACGGCGGGAAAATGGAGGGGGCTGGGGAAAACGCTTCTTTTAAGGGCGGGGAAAAGTCTTTTTAGAGGGGAAGCgggtat                                                                                                                          >  1:225567/1‑147 (MQ=255)
                        tgtgCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTAtt                                                                                                                          <  2:849749/150‑1 (MQ=255)
                           gCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCa                                                                                                                       <  2:422403/150‑1 (MQ=255)
                           gCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAg                                                                                                                      <  1:119227/151‑1 (MQ=255)
                              ctGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAgc                                                                                                                     <  2:714650/149‑1 (MQ=255)
                               tGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGCATAACAAAAGCGGTCCATTCAGGAGAACGGCGGGAAAGAGGCGGGGGCTGGGGAAACTTCTTCTTGAATGCGCGGGGACCTAGCTTTTGAGAGGCGAAGCGGGATTGTAgcgc                                                                                                                   >  2:606009/1‑150 (MQ=255)
                               tGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAgc                                                                                                                     <  2:892036/148‑1 (MQ=255)
                                ggCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGa                                                                                                                 <  1:414065/151‑1 (MQ=255)
                                 gCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGAt                                                                                                                <  1:879851/151‑1 (MQ=255)
                                  cATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAAGCGCTCCATTTAGGAGAACGGGGGGAAAGAGCAGGGGGGTGGGGAAAAAGCTTCTTTAAAACGCGGGGGCAAATCTGTTTAGAGACGAAACGGGAtt                                                                                                                          >  2:146320/1‑140 (MQ=255)
                                     aGTGAATTCACCCAGGCGTAGCGGCGGGCGGATAACAAAGGCGCTCCATTCAGGAGAACGGCGGGAAAAAGGAGGGGGCTGGGCAAAAATCTTTTTTAAAGCGCGGGGGCAAGTCTTTTTAGAGGCGAAGGGGGATTTTAGCGCGGTGTcc                                                                                                            >  1:930240/1‑151 (MQ=255)
                                     aGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTcc                                                                                                            <  2:292259/151‑1 (MQ=255)
                                      gTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCa                                                                                                           <  1:43743/151‑1 (MQ=255)
                                                    gCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCa                                                                                             <  1:575384/151‑1 (MQ=255)
                                                        aGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTc                                                                                          <  1:533646/150‑1 (MQ=255)
                                                          cggcggGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGGCTGGGCAAACTGCTTCTTGAATGCGCGGGGACATGTCTGTTGAGAGTCGAAGCGGGATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAgg                                                                                       >  1:824662/1‑151 (MQ=255)
                                                               ggcgAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAg                                                                                  <  1:565261/151‑1 (MQ=255)
                                                                   aaTACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGGCTGGGCAAACTGCTTTTTGAATGCGCGGGGACATGACTGTTGAGAGTCGAAGCGGGATTTCAACGCGATGTCCAAAATAGGACGCGCAgcacagcgcagt                                                                                 >  2:682811/1‑137 (MQ=255)
                                                                       ccAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGAcc                                                                           >  1:694779/1‑150 (MQ=255)
                                                                              cgcTCCATTCAGGAGAACGGCGGTAAAGTGCGGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTACATGTCGGTTGAGGGTCGAAGCGGTATTGCAGCGCGATGCCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGAtt                                                                    <  2:485571/150‑1 (MQ=255)
                                                                               gcTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTc                                                                  >  1:988503/1‑151 (MQ=255)
                                                                                cTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGTTACATGTCTGTTGAGAGTCGAAGCGGTATGGCAGCGCGATGTCCAGAATCGGACGGGCAGTCAGGCGTAGTGCGACCGCCGATTTc                                                                  <  2:48163/150‑1 (MQ=255)
                                                                                    aTTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACa                                                              >  1:440231/1‑150 (MQ=255)
                                                                                      tCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGGCCAGCATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACaaa                                                            >  1:824292/1‑150 (MQ=255)
                                                                                        aGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAacg                                                          >  2:264950/1‑150 (MQ=255)
                                                                                         ggAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAac                                                           <  1:782652/148‑1 (MQ=255)
                                                                                          gagaACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAGGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAacgacg                                                       <  1:806232/151‑1 (MQ=255)
                                                                                               cggcggTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCAc                                                  >  2:889444/1‑151 (MQ=255)
                                                                                                ggcggTAAAGTGCAGGATTCGGGGCAAACTGCTTTTTGAAGGCGCGGTTAAATGTCTGTTGAGAGTCGAGGCGGTATTGCAGCGGGACGTCCAGAGTCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACg                                                 <  2:469392/151‑1 (MQ=255)
                                                                                                ggcggTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGAAGTGCGACCGCCGATTTCGACAAACGACGAGCACg                                                 >  2:907203/1‑151 (MQ=255)
                                                                                                       aaGTGCAGGAGTCGGGGCAACCTGCTTTTTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATAc                                          <  1:823816/151‑1 (MQ=255)
                                                                                                              ggAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCa                                    >  1:295093/1‑150 (MQ=255)
                                                                                                                      ggCAAACTGTTTCTTGAATGCGCGGGTACTTGTCTTTGGAGAGTCGAAGCGTTATTGCAGCGCGATTTCCAGAATCGGCCGCGCAGTCAGGCGTAGTCCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATgg                           <  2:242837/151‑1 (MQ=255)
                                                                                                                         aaaCTGCTTCTTGAATGCGCGGGTACATGTCTGTTGAGAGTCAAAGCGGTATTGCAGCGCGATGTCCAGAATCGGCCGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCa                        <  1:768431/151‑1 (MQ=255)
                                                                                                                             tGCTTCTTGAATGCGCGGGTACATGTCTGTGGAGAGTCGAGGCGTTATTGCAGCGCGATGCCCAGAACCGGACGCGCAGTAAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGt                      <  2:80912/149‑1 (MQ=255)
                                                                                                                                ttcttGAATGCGGGGGTACATGCTTGTTGAGGGTCGAGGCGGTATTGCAGCGCGATGTCCAGAAGCGGACGCGCAGTCGGGCGTGGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGa                    <  2:694792/148‑1 (MQ=255)
                                                                                                                                ttcttGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACAt                 >  1:100109/1‑151 (MQ=255)
                                                                                                                                 tcttGAATGCGCGGGTACATGTCTGTTGAGAGTCGAGGCGTTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACAt                 <  1:396245/150‑1 (MQ=255)
                                                                                                                                   ttGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCtt              >  2:26548/1‑151 (MQ=255)
                                                                                                                                    tGAATGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCt               <  1:849501/149‑1 (MQ=255)
                                                                                                                                       aTGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTaa           >  2:292372/1‑150 (MQ=255)
                                                                                                                                       aTGCGCGGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTaa           >  1:96567/1‑150 (MQ=255)
                                                                                                                                           gcgGGTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACAtt      >  2:347329/1‑151 (MQ=255)
                                                                                                                                               gTACATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACATTctct  >  1:4103/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACATTCTCT  >  NC_000913/4634902‑4635195

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