Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 3,486,574 C→A A152E (GCA→GAA)  crp → cAMP‑activated global transcription factor, mediator of catabolite repression

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi3,486,5740CA97.2% 114.7 / NA 36A152E (GCA→GAA) crpcAMP‑activated global transcription factor, mediator of catabolite repression
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (18/17);  total (18/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACAC  >  NC_000913_3_bme_pgi/3486432‑3486719
                                                                                                                                              |                                                                                                                                                 
cAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAaaca                                                                                                                                             >  1:398820/1‑149 (MQ=255)
  aTTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCACAGAAAGGGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAaacaac                                                                                                                                           <  1:348958/149‑1 (MQ=255)
  aTTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGACGCGTCGTCTGCACGTCACTTCAGGGAAAGTGGGCAACCTGGCGTGCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAaaaaacac                                                                                                                                           >  2:227514/1‑145 (MQ=255)
         cAGGTAAACCCGGACATTCGGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGc                                                                                                                                    <  1:530021/149‑1 (MQ=255)
          aGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCt                                                                                                                                   >  2:18938/1‑149 (MQ=255)
            gTAAACCCGGACATTCTGATGCGTTTTTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCACCCCGGCGTTCCTCGACGTGACGGGTCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTAt                                                                                                                                 <  1:382849/149‑1 (MQ=255)
            gTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGCCGCTAt                                                                                                                                 >  2:213301/1‑149 (MQ=255)
              aaaCCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGa                                                                                                                               <  1:18938/149‑1 (MQ=255)
                  ccGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCa                                                                                                                           >  2:483359/1‑149 (MQ=255)
                                   tttGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCaaa                                                                                                          >  1:191947/1‑149 (MQ=255)
                                     tGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATc                                                                                                        <  2:328568/149‑1 (MQ=255)
                                     tGTCGGCGCAGATGGCGCGGCGTCGGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATc                                                                                                        <  2:462872/149‑1 (MQ=255)
                                                          gtcTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAAACAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTAGAGAAAtt                                                                                   >  2:369959/1‑149 (MQ=255)
                                                           tcTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTg                                                                                  <  2:217000/149‑1 (MQ=255)
                                                             tGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGACAGGAAATTgat                                                                                >  1:222703/1‑147 (MQ=255)
                                                               cAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGACCGGTATGCAAGTCAAAATTACCCGCCGGGAAATTGGtaa                                                                              >  2:162069/1‑147 (MQ=255)
                                                                  gTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATg                                                                                                                                <  1:71955/96‑1 (MQ=255)
                                                                  gTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATg                                                                                                                                >  2:71955/1‑96 (MQ=255)
                                                                        tCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCgg                                                                     >  2:148225/1‑149 (MQ=255)
                                                                        tCAGACAGAGTTGGCAACCTGTCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCgg                                                                     <  2:209416/149‑1 (MQ=255)
                                                                         cAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGc                                                                    <  1:384494/149‑1 (MQ=255)
                                                                                  tGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGt                                                           <  1:148225/149‑1 (MQ=255)
                                                                                             gCGGTCCTCGACGTGAGGAGCCGCAATGCACAGACTCTGCTGAATCTGGAAAAGCAACAAGACGCTATAAGTCACCCGGACGGTAATCAAATCAAAATTAGCAATAAGGAACGTGGATCGAGTGTCGGCTGGTCAGCTAAAACCGgggg                                                >  2:221981/1‑149 (MQ=255)
                                                                                                              gggCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGc                               <  1:528627/149‑1 (MQ=255)
                                                                                                                  cGCATTGCACAGACTCTGCTTAATCTGGAAAAACAACCCGACGCGATGACCCACACCGACGGAAGGCAAATCAAAATTACCCGTCAGCAAATTATTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTATGAAGATGCTGGa                           >  2:199396/1‑149 (MQ=255)
                                                                                                                       ggCACAGACGCTGATGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATc                      <  1:369959/148‑1 (MQ=255)
                                                                                                                         cacaGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTATTGTGTGTTCTGGGGATGCATTCTGAAGATGCTGGAAGATcgg                    >  1:423813/1‑147 (MQ=255)
                                                                                                                          acatcctctGCTGTATCGTGCAAAACACCCAGCCGCTACGACTCACCCGGACGGTATGGAAATCAAAATTACCCGTCATGAAATTGGTTAGAGTGACGGCTGTTCTCGTGAAACTGTGGGTCGCATTCTGAAGATGCTGGAAGATCAGa                   <  2:463390/144‑1 (MQ=255)
                                                                                                                             gACTCTGCTGAATCTGGAAAAACACCCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAAcc                <  1:430043/149‑1 (MQ=255)
                                                                                                                                tctGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAAcataa             >  1:36863/1‑145 (MQ=255)
                                                                                                                                tctGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGa             <  1:410593/149‑1 (MQ=255)
                                                                                                                                tctGCTGAATCTGGAAAAACAACCAGACGATATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGa             <  1:337213/149‑1 (MQ=255)
                                                                                                                                  tgctgAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGAtc           >  2:464542/1‑149 (MQ=255)
                                                                                                                                   gctgAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGAtct          >  1:526897/1‑149 (MQ=255)
                                                                                                                                     tgAATCTGGAAAAACAATCAGAGGCTATGGCGAACCTGGACGGTAGGCAAATCAAAATTACGCGTGAGGAAATTGGTGAGATTGTCGGCTGTTCGCGGGAAACCGTGGGTCGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTcc        <  2:516473/149‑1 (MQ=255)
                                                                                                                                           tGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGCTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCacac  >  1:486411/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
CAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACAC  >  NC_000913_3_bme_pgi/3486432‑3486719

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