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.C94.5% 193.0 / 4.5 55intergenic (‑43/+62)celF/chbR6‑phospho‑beta‑glucosidase/transcriptional regulator ChbR
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base C (28/24);  total (30/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.76e-01

GCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATT  >  NZ_CP009273/1812713‑1812884
                                                                                        |                                                                                    
gctgctCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGcc                                                                                    <  1:390917/90‑1 (MQ=255)
   gctCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGccccc                                                                                 >  2:2922934/1‑90 (MQ=255)
   gctCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGccccc                                                                                 >  2:2923564/1‑90 (MQ=255)
   gctCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGccccc                                                                                 >  2:487187/1‑90 (MQ=255)
    ctCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGcccccc                                                                                >  2:1010417/1‑90 (MQ=255)
    ctCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGcccccc                                                                                >  1:161619/1‑90 (MQ=255)
     tCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGccccccc                                                                               >  2:2160544/1‑90 (MQ=255)
     tCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGccccccc                                                                               >  2:1071581/1‑90 (MQ=255)
     tCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGccccccc                                                                               >  1:1177082/1‑90 (MQ=255)
      cccGCCACCAATAGTGACGATTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGcccccccc                                                                              <  1:2456164/90‑2 (MQ=255)
      cccGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGcccccccc                                                                              <  1:924116/90‑2 (MQ=255)
      cccGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGcccccccc                                                                              <  2:339157/90‑2 (MQ=255)
      cccGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGcccccccc                                                                              <  1:1494380/90‑2 (MQ=255)
                 tAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACataaata                                                                    <  1:1430596/90‑1 (MQ=255)
                    tGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATc                                                                 >  2:1358185/1‑90 (MQ=255)
                     gacgacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATcc                                                                >  1:2478365/1‑90 (MQ=255)
                      acgacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCa                                                               <  1:1358185/90‑1 (MQ=255)
                       cgacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAg                                                              >  2:2667346/1‑90 (MQ=255)
                       cgacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAg                                                              >  2:2693513/1‑90 (MQ=255)
                        gacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGc                                                             >  2:2950454/1‑90 (MQ=255)
                        gacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGc                                                             >  2:2510078/1‑90 (MQ=255)
                        gacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGc                                                             >  2:2780914/1‑90 (MQ=255)
                         acTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCa                                                            <  1:100854/90‑1 (MQ=255)
                           ttttAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAAc                                                          <  1:1010417/90‑1 (MQ=255)
                           ttttAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAAc                                                          <  2:2182773/90‑1 (MQ=255)
                           ttttAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAAc                                                          <  1:498124/90‑1 (MQ=255)
                           ttttAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAAc                                                          <  2:2255958/90‑1 (MQ=255)
                                aTTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGAc                                                     <  1:2922934/90‑1 (MQ=255)
                                aTTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGAc                                                     <  1:1071581/90‑1 (MQ=255)
                                aTTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGAc                                                     <  2:887396/90‑1 (MQ=255)
                                aTTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGAc                                                     <  2:1352458/90‑1 (MQ=255)
                                  ttCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAg                                                   <  2:513433/90‑1 (MQ=255)
                                  ttCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAg                                                   <  2:324875/90‑1 (MQ=255)
                                            tAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAAt                                         <  1:2160544/90‑1 (MQ=255)
                                                ttCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTc                                     <  2:1133889/90‑1 (MQ=255)
                                                 tctcCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCa                                    <  1:1103060/90‑1 (MQ=255)
                                                  ctcCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAg                                   <  2:1529093/90‑1 (MQ=255)
                                                  ctcCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAg                                   <  1:1578484/90‑1 (MQ=255)
                                                    cccTTCAGTTCAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGAc                                                     >  2:1472424/1‑70 (MQ=255)
                                                    cccTTCAGTTCAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGAc                                                     <  1:1472424/70‑1 (MQ=255)
                                                      cTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGtat                               >  1:103686/1‑90 (MQ=255)
                                                       ttCAGTACAGAATACTGATATCTGGCATATCTCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAAc                            >  2:506462/1‑90 (MQ=255)
                                                               aGAATACTGATATCTGGCATATCTCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACt                    <  1:506462/90‑1 (MQ=255)
                                                                gAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTAc                     >  2:94533/1‑90 (MQ=255)
                                                                gAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTAc                     >  1:625974/1‑90 (MQ=255)
                                                                gAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTAc                     <  1:487187/90‑1 (MQ=255)
                                                                gAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTAc                     >  1:1010828/1‑90 (MQ=255)
                                                                     cTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGa                                                      <  2:1697860/52‑1 (MQ=255)
                                                                     cTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGa                                                      >  1:1697860/1‑52 (MQ=255)
                                                                     cTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGAtt                >  1:1048844/1‑90 (MQ=255)
                                                                     cTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGAtt                >  1:2527543/1‑90 (MQ=255)
                                                                     cTGATATCTGGCATATCTGCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTa               >  1:1363958/1‑90 (MQ=255)
                                                                       gATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTaa              >  2:705738/1‑90 (MQ=255)
                                                                       gATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTaa              >  1:2374362/1‑90 (MQ=255)
                                                                          atCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGat                                                 <  2:945454/52‑1 (MQ=255)
                                                                          atCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGat                                                 >  1:945454/1‑52 (MQ=255)
                                                                          atCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAAtt           >  1:2428644/1‑90 (MQ=255)
                                                                          atCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAAtt           >  1:2898748/1‑90 (MQ=255)
                                                                          atCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAAtt           >  1:168942/1‑90 (MQ=255)
                                                                          atCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAAtt           >  2:971676/1‑90 (MQ=255)
                                                                          atCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAAtt           >  2:1316956/1‑90 (MQ=255)
                                                                            cTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAg                                   >  2:1227639/1‑64 (MQ=255)
                                                                            cTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAg                                   <  1:1227639/64‑1 (MQ=255)
                                                                            cTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTc          >  1:1104816/1‑89 (MQ=255)
                                                                            cTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTc          <  2:1104816/89‑1 (MQ=255)
                                                                            cTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCa         >  1:253439/1‑90 (MQ=255)
                                                                                   atCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGCCAGATATTTGAATTGTCAGGGATAACGACTTACTGATTAAATTCAGTCAAtt  >  1:1736406/1‑90 (MQ=255)
                                                                                   atCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAAtt  >  1:1323175/1‑90 (MQ=255)
                                                                                        |                                                                                    
GCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATT  >  NZ_CP009273/1812713‑1812884

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