Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,265,386 G→T 18.3% intergenic (‑417/‑13) yjcB ← / → yjcC putative inner membrane protein/putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,265,3860GT18.3% 48.6 / 4.1 22intergenic (‑417/‑13)yjcB/yjcCputative inner membrane protein/putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (9/9);  new base T (2/2);  total (11/11)
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.64e-01

TGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGAAATTAATTTTGAATATCAATGAATTATTTTCATCCAGGTGACGATTAGAAAGGTATCAATTTCAAATCAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGAACCTTTGCTCAACTGGCACTGGATAAATCCGAGCTGGTCATTCGCCAGGCAGATTTAGTGAGCGATGCAGCTGAACGCTATCAGGGGCAAGTTTGCACTC  >  CP009273/4265145‑4265627
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
tGACTGTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATCTCTTCCCGCCGCCTCTGCATTCCTGTAGGAAATTAATTTTGAATATCAATGAATTATTTTCATCCAGGTGACGATTAGAAAGGTATCAATTTCAAATCAGGCAAAAGTGCTATTTATACTGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACTTTGCAATGTCTGGTTATAAACAAGCAGGAAATTGATATCAATGAGTCATCGTGCCCGaaac                                                                                                                                                                                                                  >  2:108549/1‑275 (MQ=255)
                        tgtgGTTTTTCAGTTGCATGTAGGGCATATTTTCCCTCCGCCTCTTCATTCCTGTATGAAATTAATTTTGAATATCAATGAATTATTTTCATCCAGGTGACGATTAGAAAGGTATCAATTTCAAATCAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCAtt                                                                                                                                                                                          <  1:172922/275‑1 (MQ=255)
                                                                                ggAAATTAATTTTGAATATCAATGAATTATTTTCATCCAGGTGACGATTAGAAAGGTATCAATTTCAAATCAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGTAAATTGATAGCaa                                                                                                                                                                                                                                      >  1:348440/1‑175 (MQ=255)
                                                                                    attaattTTGAATATCAATGAATTATTTTCATCCAGGTGACGATTAGAAAGGTATCAATTTCAAATCAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGTAAATTGATAGCAATGAg                                                                                                                                                                                                                                  <  2:348434/175‑1 (MQ=255)
                                                                                                             ttttCATCCAGGTGACGATTAGAAAGGTATCAATTTCAAATCAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTAGCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGTATTACCTTCCTTTGGCCAACATCAGAAGTGCATAATCAGCTcctcc                                                                                                     >  2:347003/1‑272 (MQ=255)
                                                                                                                 cATCCAGGTGACGATTAGAAAGGTATTAATTTCAAATCAGGCAAAAGTGCTATTTATACCGTAAGAATTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAAGNGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGaataat                                                                                                               <  1:347009/261‑1 (MQ=255)
                                                                                                                                                         gCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGAACCTTTGCTCAACTGGCACTGGATAAATCCGAGCTGGTCATTCGc                                                         >  1:24412/1‑275 (MQ=255)
                                                                                                                                                                               aaGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTATCTTTCCCATCATTCTTTCGCTATGGATTTCCTTCCTTTGTGCAAAATCAgca                                                                                                                        >  2:32849/1‑188 (MQ=255)
                                                                                                                                                                                aGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTATCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTTGGCAAAATCAGAAGTGAATAATCAGc                                                                                                           >  2:23396/1‑202 (MQ=255)
                                                                                                                                                                                   tttATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGa                                                                                                                    <  1:32850/190‑1 (MQ=255)
                                                                                                                                                                                    ttATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCg                                                                                                       <  1:23397/202‑1 (MQ=255)
                                                                                                                                                                                      aTCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGAACCTTTGCTCAACTGGCACTGGATAAATCCGAGCTGGTCATTCGCCAGGCAGATTTAGTGAGCGATGCAGCTGa                            <  1:108553/275‑1 (MQ=255)
                                                                                                                                                                                                 tCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCttttt                                                                                                                                                                                  >  2:70107/1‑114 (MQ=255)
                                                                                                                                                                                                     tACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGtt                                                                                                                                                                              <  1:70110/114‑1 (MQ=255)
                                                                                                                                                                                                        ccAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGTAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTtc                                                                                                                                                                             >  1:19467/1‑112 (MQ=255)
                                                                                                                                                                                                            ggTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGTAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCtt                                                                                                                                                                         <  2:19466/112‑1 (MQ=255)
                                                                                                                                                                                                               tttCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGAACCTTTGCTCAACTGGCACTGGATAAATCCGAGCTGGTCATTCGCCAGGCAGATTTAg                                            >  1:159703/1‑234 (MQ=255)
                                                                                                                                                                                                                   aCCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGAACCTTTGCTCAACTGGCACTGGATAAATCCGAGCTGGTCATTCGCCAGGCAGATTTAGTGAg                                        <  2:159698/234‑1 (MQ=255)
                                                                                                                                                                                                                               ccGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGAACCTTTGCTCAACTGGCACTGGATAAATCCGAGCTGGTCATTCGCCAGGCAGATTTAg                                            >  1:61604/1‑218 (MQ=255)
                                                                                                                                                                                                                                   gTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTtc                                                                                                                                                                             >  1:39627/1‑85 (MQ=255)
                                                                                                                                                                                                                                   gTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGAACCTTTGCTCAACTGGCACTGGATAAATCCGAGCTGGTCATTCGCCAGGCAGATTTAGTGAg                                        <  2:61601/218‑1 (MQ=255)
                                                                                                                                                                                                                                       aaaCAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCtt                                                                                                                                                                         <  2:39626/85‑1 (MQ=255)
                                                                                                                                                                                                                                             gcaggAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGAACCTTTGCTCAACTGGCACTGGATAAATCCGAGCTGGTCATTCGCCAGGCAGATTTAGTGAGCGATGCAGCTGAACGCTATCAGGGGCAAGTTTGc      >  1:329202/1‑242 (MQ=255)
                                                                                                                                                                                                                                                 gAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGAACCTTTGCTCAACTGGCACTGGATAAATCCGAGCTGGTCATTCGCCAGGCAGATTTAGTGAGCGATGCAGCTGAACGCTATCAGGGGCAAGTTTGCACTc  <  2:329196/242‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
TGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGAAATTAATTTTGAATATCAATGAATTATTTTCATCCAGGTGACGATTAGAAAGGTATCAATTTCAAATCAGGCAAAAGTGCTATTTATACCGTAAGATTTATCTAAAGACGTCGGTACCCAGGGTTTTCACCTTGCAATGGCCGGGTATAAACAGGCAGGAAATTGATAGCAATGAGTCATCGTGCACGACACCAATTACTGGCGTTGCCGGGCATTATCTTTTTAGTTCTCTTTCCCATCATTCTTTCGCTATGGATTGCCTTCCTTTGGGCAAAATCAGAAGTGAATAATCAGCTCCGAACCTTTGCTCAACTGGCACTGGATAAATCCGAGCTGGTCATTCGCCAGGCAGATTTAGTGAGCGATGCAGCTGAACGCTATCAGGGGCAAGTTTGCACTC  >  CP009273/4265145‑4265627

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