Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,812,807 (C)7→8 intergenic (‑50/+55) celF ← / ← chbR 6‑phospho‑beta‑glucosidase/transcriptional regulator ChbR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,812,8001.C97.9% 153.6 / ‑3.0 48intergenic (‑43/+62)celF/chbR6‑phospho‑beta‑glucosidase/transcriptional regulator ChbR
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (29/18);  total (30/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 = 4.74e-01

TAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAAT  >  NZ_CP009273/1812711‑1812883
                                                                                          |                                                                                    
taGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTg                                                                                       >  2:1054610/1‑90 (MQ=255)
taGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTg                                                                                       >  1:1318343/1‑90 (MQ=255)
     gctCCCGCCACAAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑‑ccccc                                                                                >  2:312579/1‑90 (MQ=255)
      ctCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑‑cccccc                                                                               >  2:155425/1‑90 (MQ=255)
       tCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑‑ccccccc                                                                              >  1:993759/1‑90 (MQ=255)
        cccGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑‑cccccccc                                                                             <  2:1506036/90‑2 (MQ=255)
             caccaATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACa                                                                         <  2:1208962/90‑1 (MQ=255)
             caccaATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACa                                                                         <  2:187030/90‑1 (MQ=255)
             caccaATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACa                                                                         <  2:1057113/90‑1 (MQ=255)
             caccaATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACa                                                                         <  1:661752/90‑1 (MQ=255)
                      tGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATc                                                                >  1:417785/1‑90 (MQ=255)
                       gacgacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCCGGACATAAATAATc                                                                >  1:99730/1‑90 (MQ=255)
                       gacgacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATcc                                                               >  2:742908/1‑90 (MQ=255)
                        acgacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCa                                                              <  1:1212786/90‑1 (MQ=255)
                                  aTTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGAc                                                    <  2:1314910/90‑1 (MQ=255)
                                  aTTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGAc                                                    <  2:467438/90‑1 (MQ=255)
                                        ggCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGatat                                              >  1:1093431/1‑90 (MQ=255)
                                         gCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATAtg                                             <  2:1487790/90‑1 (MQ=255)
                                              tAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCGACAGGACAGATATGTGAAt                                        <  1:449116/90‑1 (MQ=255)
                                              tAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAAt                                        <  2:478441/90‑1 (MQ=255)
                                                   tctcCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCa                                   >  1:1335902/1‑90 (MQ=255)
                                                   tctcCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCa                                   <  2:976068/90‑1 (MQ=255)
                                                    ctcCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAg                                  <  1:742908/90‑1 (MQ=255)
                                                      cccTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGt                                >  1:1180820/1‑90 (MQ=255)
                                                        cTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGtat                              >  2:314977/1‑90 (MQ=255)
                                                          tCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATaa                            >  2:424164/1‑90 (MQ=255)
                                                          tCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATaa                            >  1:773805/1‑90 (MQ=255)
                                                                 aGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAGTAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTa                     <  2:391296/90‑1 (MQ=255)
                                                                  gAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGAATAACGACTTAc                    >  2:698794/1‑90 (MQ=255)
                                                                   aaTACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACt                   <  1:609947/90‑1 (MQ=255)
                                                                       cTGATATCTGGCATATCTGCCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGAt                >  2:98495/1‑90 (MQ=255)
                                                                       cTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGat                                                >  2:491452/1‑57 (MQ=255)
                                                                       cTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGat                                                <  1:491452/57‑1 (MQ=255)
                                                                       cTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGAtt               >  2:1100831/1‑90 (MQ=255)
                                                                       cTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGAtt               >  2:494122/1‑90 (MQ=255)
                                                                       cTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGAtt               <  2:331715/90‑1 (MQ=255)
                                                                       cTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGAtt               >  1:493677/1‑90 (MQ=255)
                                                                         gATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTaa             >  2:372914/1‑90 (MQ=255)
                                                                            atCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGAc                                                    <  2:865555/48‑1 (MQ=255)
                                                                            atCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGAc                                                    >  1:865555/1‑48 (MQ=255)
                                                                            atCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAAtt          >  2:1495446/1‑90 (MQ=255)
                                                                            atCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAAtt          >  2:837864/1‑90 (MQ=255)
                                                                            atCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAAtt          >  2:324408/1‑90 (MQ=255)
                                                                            atCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAAtt          >  2:1320633/1‑90 (MQ=255)
                                                                            atCTGGCATATCCGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAAtt          >  1:896812/1‑90 (MQ=255)
                                                                              cTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCa        >  1:394083/1‑90 (MQ=255)
                                                                              cTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCa        >  1:1482173/1‑90 (MQ=255)
                                                                              cTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCa        >  2:963434/1‑90 (MQ=255)
                                                                              cTGGCATATCTG‑‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAg       >  1:1219689/1‑90 (MQ=255)
                                                                                ggCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGt      >  2:486661/1‑90 (MQ=255)
                                                                                 gCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTc     >  1:871051/1‑90 (MQ=255)
                                                                                 gCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTc     >  2:1547347/1‑90 (MQ=255)
                                                                                 gCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTc     >  2:899068/1‑90 (MQ=255)
                                                                                    tatCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAAt  <  1:1054610/90‑1 (MQ=255)
                                                                                          |                                                                                    
TAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAAT  >  NZ_CP009273/1812711‑1812883

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