Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 2,198,177 A→G intergenic (‑40/+35) U712_11010 ← / ← U712_11015 Chorismate synthase/Chemotaxis protein methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,198,1770AG96.9% 89.4 / NA 32intergenic (‑40/+35)U712_11010/U712_11015Chorismate synthase/Chemotaxis protein methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (22/9);  minor base T (1/0);  total (23/9)
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

GATATCTTCCTCCGTTATGTAAAGCCCGGCAGGCACACCTTCAATAATGGTTGTCAGTTGGGGGCCGTGTGATTCTCCGGCTGTTAAATATCTCATGACTTCTTCTCCCTTCATCGCATTAACGCTTTTAAGTAAAACTACTATAACATATTTTTCAGCAAAAAGTCAGTCTATCTTTTTTGATAAAAGAATGTATCGGCTGGAGCCAGGCCAAATTTCTCCGGATTGAAGATTTGCTCGGTGCTGCCGACAAATAGGACTCCATT  >  CP006881/2198042‑2198307
                                                                                                                                       |                                                                                                                                  
gATATCTTCCTCCGTTATGTAAAGCCCGGCAGGTACACCCTCAATAATGGTTGTCAGTTGGGGGCCGTGTGATTCTCCGTCTGTTAAATATCTCATGACTTCTTCTCCCTTCATCGCATTAACGCTTTTAAGTAAGac                                                                                                                                  <  5:174228/138‑1 (MQ=255)
                 tGTAAAGCCCGGCAGGTACACCCTCAATAATGGTTGTCAGTTGGGGGCCGTGTGATTCTCCGGCTGTTAAATATCTCATGACTTCTTCTCCCGTCATCGCATTAACGCTTTTAAGTAAGACTACTATAACATAttttt                                                                                                                 >  3:379362/1‑138 (MQ=255)
                        cccGGCAGGTACACCCTCAATAATGGTTGTCAGTTGGGGGCCGTGTGATTCTCCGGCTGTTAAATATATCATGACTTGTTCTCCCTTCATACCAATAACGCTTTTAATTAAGACTACTATAACATATTTTTCAGCaaa                                                                                                          >  7:357518/1‑138 (MQ=255)
                           ggcaggtaCACCCTCAATAATGGTTGTCAGTTGGGGGCCGTGTGATTCTCCGGCTGTTTAATTTCTCATGACTTATTGTACCTTGATAGCATTAACGCTTTTAAGTAAGACTACTATAACATATTTTTCAGCAAAAagt                                                                                                      >  4:168790/1‑139 (MQ=255)
                           ggcaggtaCACCCTCAATAATGGTTGTCAGTTGGGGGCCGTGTGATTCTCCGGCTGTTAAATATCTCATGACTTCTTCTCCCTTCATCGCATTAACGCTTTTAAGTAAGACTACTATAACATATTTTTCAGCAAAAagt                                                                                                      >  8:243000/1‑139 (MQ=255)
                                       ctCAATAATGGTTGTCAGTTGGGGGCCGTGTGATTCTCCGGCTGTCAAATATCTCATGAAGTCTTGTCCCTTCATCGCATTAACGCTTTTTAGTAAGCATACTATAACATATTTTTCAGCAAAAAGTCATTGTATAttt                                                                                          >  5:368623/2‑139 (MQ=255)
                                          aataatGGTTGTCAGTTGGGGGCCGTGTGATTCTCCGGCTGTTAAATATCTCATGACTTCTTATCCCTTCTTCGCATTAAAGCTTTTAATTAAGACTACTATAACATATATTTACGCAAAAAGTCAGTCTATCtttttt                                                                                       >  5:177379/1‑139 (MQ=255)
                                               tGGTTGTCAGTTGGGGGCCGTGTGATTCTCCGGCTGTTAAATATATCATTACTTATTCGCCCTTCATAGCATTAACGCTTTTTAGTAAGAATACTATAACATATTTTTAAGCAAAAAGTCAGTCTATATTTTTTGATa                                                                                   >  2:231756/1‑138 (MQ=255)
                                                  ttGTCAGTTGGGGGCCGTGTGATTCTCCGGCTGTTAAATATCTCATGACTTCTTCTCCCTTCATCGCATTAACGCTTTTAAGTAAGACTACTATAACATATTTTTCAGCAAAAAGTCAGTCTATCTTTTTTGATaa                                                                                  >  7:326692/1‑136 (MQ=255)
                                                   tGTCAGTTGGGGGCCGTGTGATTCTCCGGCTGTTAAATATCTCATGACTTCTTCTCCCTTCATCGCATTAACGCTTTTAAGTAAGACTACTATAACATATTTTTCAGCAAAAAGTCAGTCTATCTTTTTTGATAAAAg                                                                               >  4:184406/1‑138 (MQ=255)
                                                    gTCAGTTGGGGGCCGTGTGGTTGTCCGGCGGTTTAATATCTCATGGCTTCTTCTCCCGTCATAGCATTTACGCGTTTAAGTTAGACTACTATAACATATTTTTAAGCAAAAAGTCAGTCTATATTTTTTGTTAAAATaa                                                                             >  6:364503/1‑139 (MQ=255)
                                                    gTCAGTTGGGGGCCGTGTGATTCTCCGGCTGTTAAATATATCATGACTTCTTCTCCCTTCATCGCATTAACGCTTTTAAGTAATACTAAGATAACATATTTTTGAGCAAAAATTCAGTCTATCTTTTTTGTTAAAAGaa                                                                             >  8:343875/1‑139 (MQ=255)
                                                    gTCAGTTGGGGGCCGTGTGATTATACGGCTGTTAAATATCTCATGACTTCTTCTCCCTTCATCGCATTAACGCTTTTAAGTAAGACTACTATAACATATTTTTCAGCAAAAAGTCAGTCTATCTTTTTTGATAAAAGaa                                                                             >  1:39613/1‑139 (MQ=255)
                                                        gTTGGGGGCCGTGTGATTCTCCGGCTGTTAAATATCTCATGACTTCTTCTCCCTTCATCGCATTAACGCTTTTAAGTAAGACTACTATAACATATTTTTCAGCAAAAAGTCAGTCTATCTTTTTTGATAAAAGAATGTa                                                                         <  3:184406/139‑1 (MQ=255)
                                                                   tgtgATTCTCCGGCTGTTAAATATCTCATGACTTCTTCTCACTTCATCGCATTAACGCTTTTAAGTAAGACTACTATAACATATTTTTCAGCAAAAAGTCAGTCTATCTTTTTTGATAAAAGAATGTATCGGCTgggac                                                              <  1:69437/139‑4 (MQ=255)
                                                                             cGGCTGTTAAATATCTCATGACTTCTTCTCCCTTCATCGCATTAACGCTTTTAAGTAAGACTACTATAACATATTTTTCAGCAAAAAGTCAGTCTATCTTTTTTGATAAAAGAATGTATCGGCTgggac                                                              <  4:172289/129‑4 (MQ=255)
                                                                             cGGCTGTTAAATATCTCATGACTTCTTCTCCCTTCATCGCATTAACGCTTTTAAGTAAGACTACTATAACATATTTTTCAGCAAAAAGTCAGTCTATCTTTTTTGATAAAAGAATGTATCGGCTgggac                                                              >  3:172289/1‑126 (MQ=255)
                                                                             cGGCTGTTAAATATCTCATGACTTCTTCTCCCTTCATCGCATTAACGCTTTTAAGTAAGACTACTATAACATATTTTTCAGCAAAAAGTCAGTCTATCTTTTTTGATAAAAGAATGTATCGGCTGGGACTAGGCCAAAt                                                    >  8:282002/1‑139 (MQ=255)
                                                                               gCTGTTAAATATCTCATGACTTCTTCTCCCTTCATCGCATTAACGCTTTTAAGTAAGACTACTATAACATATTTTTCAGCAAAAAGTCAGTCTATCTTTTTTGATAAAAGAATGTATCGGCTGGGACTAGGCCAAAtt                                                   >  1:178772/1‑138 (MQ=255)
                                                                                 tGTTAAATATCTCATGACTTCTACACCCTTCATCGCATTAACGCTTTTAAGTAAGACTACTATAACATATTTTTAAGCAAAAAGTCAGACTAGCTTTTTTGATAAAAGAAAGTATCTGATGGTACAACGCAAAATTtct                                                >  8:405160/1‑139 (MQ=255)
                                                                                          tctcATGACTTCTTCTCCCTTCATCGCATTAACGCTTTTAAGTAAGACTACTATAACATATTTTTCAGCAAAAAGTCAGTCTATCTTTTTTGATAAAAGAATGTATCGGCTGGGACTAGGCCAAATTTCTCGGGATTg                                        <  4:379362/138‑1 (MQ=255)
                                                                                             cATGACTTCTTCTCCCTTCATCGCATTAACGCTTTTAAGTAAGACTACTATAACATATTTTTCAGCAAAAAGTCAGTCTATCTTTTTTGATAAAAGAATGTATCGGCTGGGACTAGGCCAAATTTCTCGGGATTg                                        >  2:134173/1‑135 (MQ=255)
                                                                                                     cttctCCCTTCATCGCATTAACGCTTTTAAGTAAGACTACTATAACATATTTTTCAGCAAAAAGTCAGTCTATCTTTTTTGATAAAAGAATGTATCGGCTGGGACTAGGCCAAATTTCTCGGGATTGAAGATTTGCt                              <  2:178772/137‑1 (MQ=255)
                                                                                                     cttctCCCTTCATCGCATTAACGCTTTTAAGTAAGACTACTATAACATATTTTTCAGCAAAAAGTCAGTCTATCTTTTTTGATAAAAGAATGTATCGGCTGGGACTAGGCCAAATTTCTCGGGATTGAAGATTTGCTCg                            >  8:239346/1‑139 (MQ=255)
                                                                                                     cttctCCATTCATCGCATTAACGCTTTTAAGTAAGACTACTATAACATATTTTTCAGCAAAAAGTCAGTCTATCTTTTTTGATAAAAGAATGTATCGGCTGGGACTAGGCCAAATTTCTCGGGATTGAAGATTTGCt                              <  5:28367/137‑1 (MQ=255)
                                                                                                      ttctCCCTTCATCGCATTAACGCTTTTAAGTAAGACTACTATAACATATTTTTCAGCAAAAAGTCAGTCTATCTTTTTTGATAAAAGAATGTATCGGCTGGGACTAGGCCAAATTTCTCGGGATTGAAGATTTGCTCgg                           >  6:313704/1‑139 (MQ=255)
                                                                                                             ttCATCGCATTAACGCTTTTAAGTAAGACTACTATAACATATTTTTCAGCAAAAAGTCAGTCTATCTTTTTTGATAAAAGAATGTATCGGCTGGGACTAGGCCAAATTTCTCGGGATTGAAGATTTGCTCGGTGCTGcc                    >  7:154119/1‑139 (MQ=255)
                                                                                                                    cATTAACGCTTTTAAGTAAGACTACTATAACATATTTTTCAGCAAAAAGTCAGTCTATCTTTTTTGATAAAAGAATGTATCGGCTGGGACTAGGCCAAATTTCTCGGGATTGAAGATTTGCTCGGTGCTGCCGACaaaa             <  7:282002/139‑2 (MQ=255)
                                                                                                                               ttAAGTAAGACTACTATAACATATTTTTCAGCAAAAAGTCAGTCTATCTTTTTTGATAAAAGAATGTATCGGCTGGGACTAGGCCAAATTTCTCGGGATTGAAGATTTGCTCGGTGCTGCCGACAAAAAGAACTCCAtt  >  5:153629/1‑139 (MQ=255)
                                                                                                                               ttAAGTAAGACTACTATAACATATTTTTCAGCAAAAAGTCAGTCTATCTTTTTTGATAAAAGAATGTATCGGCTGGGACTAGGCCAAATTTCTCGGGATTGAAGATTTGCTCGGTGCTGCCGACAAAAAGAACTCCAtt  >  2:29975/1‑139 (MQ=255)
                                                                                                                               ttAAGTAAGACTACTATAACATATTTTTCAGCAAAAAGTCAGTCTATCTTTTTTGATAAAAGAATGTATCGGCTGGGACTAGGCCAAATTTCTCGGGATTGAAGATTTGCTCGGTGCTGCCGACAAAAAGAACTCCAtt  >  1:433449/1‑139 (MQ=255)
                                                                                                                                taagtaagTCTATTTTAACATATTTTTCAGAATAAAGTCAGTATTTCTTTTTTGATAAAAGAATGTAGCGGCTGGGTCTAGGCCATATTTCTGGGGATTGAAGATTTGCTCGGTGCTGCCGACAAATAGAACTGCAtt  <  5:200524/138‑1 (MQ=255)
                                                                                                                                taagtaagACTACTATAACATATTTTTCAGCAAAAAGTCAGTCTATCTTTTTTGATAAAAGAATGTATCGGCTGGGACTAGGCCAAATTTCTCGGGATTGAAGATTTGCTCGGTGCTGCCGACAAAAAGAACTCCAtt  <  8:154119/138‑1 (MQ=255)
                                                                                                                                taagtaagACTACTATAACATATTTTTCAGCAAAAAGTCAGTCTATCTTTTTTGATAAAAGAATGTATCGGCTGGGACTAGGCCAAATTTCTCGGGATTGAAGATTTGCTCGGTGCTGCCGACAAAAAGAACTCCAtt  <  1:231756/138‑1 (MQ=255)
                                                                                                                                tAATTAAGACTACTATAACATATTTTTCAGCAAAAAGTCAGTCTATCTTTTTTGATAAAAGAATGTATCGGCTGGGACTAGGCCAAATTTCTCGGGATTGAAGATTTGCTCGGTGCTGCCGACAAAAAGAACTCCAtt  <  8:357518/138‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                  
GATATCTTCCTCCGTTATGTAAAGCCCGGCAGGCACACCTTCAATAATGGTTGTCAGTTGGGGGCCGTGTGATTCTCCGGCTGTTAAATATCTCATGACTTCTTCTCCCTTCATCGCATTAACGCTTTTAAGTAAAACTACTATAACATATTTTTCAGCAAAAAGTCAGTCTATCTTTTTTGATAAAAGAATGTATCGGCTGGAGCCAGGCCAAATTTCTCCGGATTGAAGATTTGCTCGGTGCTGCCGACAAATAGGACTCCATT  >  CP006881/2198042‑2198307

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: