Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,634,895 +A coding (416/870 nt) 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,634,8931.A98.2% 210.9 / ‑2.8 54coding (418/870 nt)robright oriC‑binding transcriptional activator, AraC family
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base A (23/30);  total (23/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 = 8.36e-01

CGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACC‑AATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTT  >  NC_000913/4634753‑4635034
                                                                                                                                             |                                                                                                                                             
cGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAgcgg                                                                                                                                      <  2:88446/151‑1 (MQ=255)
cGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAgc                                                                                                                                        <  2:244669/149‑1 (MQ=255)
 ggACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAgcggc                                                                                                                                     <  1:392073/151‑1 (MQ=255)
 ggACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAgcg                                                                                                                                       >  1:866014/1‑149 (MQ=255)
     gcgTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTAt                                                                                                                                 >  2:257025/1‑151 (MQ=255)
       gTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCt                                                                                                                               >  2:520732/1‑151 (MQ=255)
        tttCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCtt                                                                                                                              >  2:472872/1‑151 (MQ=255)
           cattcaGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCa                                                                                                                           <  1:146401/151‑1 (MQ=255)
           cattcaGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCa                                                                                                                           <  2:465114/151‑1 (MQ=255)
               caGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGt                                                                                                                       >  2:577764/1‑151 (MQ=255)
                aGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTg                                                                                                                      <  1:530074/151‑1 (MQ=255)
                   ccGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACa                                                                                                                   <  1:872757/151‑1 (MQ=255)
                     gTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACC‑AATCAGCGGCGTATCTTCCAGGGTGACAAAt                                                                                                                <  1:237750/151‑1 (MQ=255)
                            aCCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTtgt                                                                                                            <  2:436097/149‑1 (MQ=255)
                               ggcggAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGctct                                                                                                       >  1:240474/1‑151 (MQ=255)
                               ggcggAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGctct                                                                                                       <  2:35361/151‑1 (MQ=255)
                                    aaTGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCAt                                                                                                  <  1:820371/151‑1 (MQ=255)
                                     aTGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATa                                                                                                 <  2:145464/151‑1 (MQ=255)
                                           ggCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAAt                                                                                           >  2:757402/1‑151 (MQ=255)
                                             cgcgTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTc                                                                                         <  2:278027/151‑1 (MQ=255)
                                               cgTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCAc                                                                                       >  2:699650/1‑151 (MQ=255)
                                                    ccGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAgg                                                                                  <  2:840311/151‑1 (MQ=255)
                                                      gagaAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCg                                                                                <  2:213188/151‑1 (MQ=255)
                                                         aaaaTCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGt                                                                               >  2:76548/1‑149 (MQ=255)
                                                            aTCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAgcgg                                                                          >  2:590054/1‑151 (MQ=255)
                                                                     gAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAAt                                                                 <  1:962351/151‑1 (MQ=255)
                                                                      aaCTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCggcgggcggaa                                                                 >  2:209784/1‑147 (MQ=255)
                                                                      aaCTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATa                                                                <  1:609871/151‑1 (MQ=255)
                                                                           aTAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCggcggggcgaaa                                                                >  2:261468/1‑140 (MQ=255)
                                                                           aTAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCaa                                                            >  2:41070/1‑150 (MQ=255)
                                                                                    ttCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTcc                                                   <  1:237635/150‑1 (MQ=255)
                                                                                      cATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAGTACCACAAgcg                                                       >  2:491391/1‑144 (MQ=255)
                                                                                       aTGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAAAAACCAAAGCGCCCCAtt                                                >  1:120098/1‑150 (MQ=255)
                                                                                         ggCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAg                                             <  2:331557/151‑1 (MQ=255)
                                                                                           cGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGGATACCCAAAGCGCTCCATTCAgg                                            >  1:901529/1‑150 (MQ=255)
                                                                                              aaaTCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGaa                                        <  1:855033/151‑1 (MQ=255)
                                                                                              aaaTCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAAGCGCCCCATTCAGgag                                          >  1:491679/1‑149 (MQ=255)
                                                                                               aaTCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAAc                                       <  1:831446/151‑1 (MQ=255)
                                                                                                aTCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAAcg                                      <  2:216023/151‑1 (MQ=255)
                                                                                                  cAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAAGCGCTCCATTCAGGAGaa                                        >  1:23959/1‑147 (MQ=255)
                                                                                                      gaTTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTa                                <  1:794150/151‑1 (MQ=255)
                                                                                                      gaTTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTa                                <  2:10180/151‑1 (MQ=255)
                                                                                                           gCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATAACAAAAGCGCTCCATTCAGGAGAACGGCAGGAAAggg                           >  1:386040/1‑149 (MQ=255)
                                                                                                            cTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGGATACCAAAGGAGCGCCATTCAGGAGAAcagcgg                                  >  2:761206/1‑143 (MQ=255)
                                                                                                                       aGGAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTggg               <  1:779264/151‑1 (MQ=255)
                                                                                                                        ggAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGc              <  2:381323/151‑1 (MQ=255)
                                                                                                                         gAGTAGCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAAGCGCTCCATTCAGGAGAACGGCGGGGAAAAGGAGGAGGCTggg               >  1:942891/1‑149 (MQ=255)
                                                                                                                              gCTCTGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTg        <  1:248043/151‑1 (MQ=255)
                                                                                                                                 ctGGGTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGc       <  2:886250/149‑1 (MQ=255)
                                                                                                                                   gggTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCATGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGcttct   <  1:480547/151‑1 (MQ=255)
                                                                                                                                   gggTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGctt     <  1:703239/149‑1 (MQ=255)
                                                                                                                                   gggTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAAGCGCTCCCTTCAGGAGAACGGGGGGAAAAGGCAGGAGGCTGGGGAAAACGctt     >  1:55913/1‑149 (MQ=255)
                                                                                                                                    ggTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCCTTCAGGAGAACGGCGGGAAAAAGCAGGGGGCTGGGCCAAATGcttctt  >  1:890426/1‑151 (MQ=255)
                                                                                                                                    ggTAACACCAAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGcttctt  <  2:223045/151‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                             
CGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACC‑AATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTT  >  NC_000913/4634753‑4635034

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