Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,995,231 (A)7→8 intergenic (‑139/+83) ygeK ← / ← ygeN pseudogene, response regulator family, part of T3SS PAI ETT2 remnant, putative 2‑component transcriptional regulator/pseudogene, orgB family, part of T3SS PAI ETT2 remnant

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,995,2241.A88.2% 44.4 / 0.2 17intergenic (‑132/+90)ygeK/ygeNpseudogene, response regulator family, part of T3SS PAI ETT2 remnant, putative 2‑component transcriptional regulator/pseudogene, orgB family, part of T3SS PAI ETT2 remnant
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base A (8/7);  total (9/8)
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

ATCCTGAGTGTAAATGCTTTAACGCATTCATTTATAATCATCCCTTCCATTATTATTATAAGCAAAATCCAAAGAATACATTGATGAAATAATAATGAAATATAATTAAAAATAAAATTTTTGCGTAAAAAAATACCACAGGCATTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAA  >  NC_000913/2995099‑2995344
                                                                                                                              |                                                                                                                        
aTCCTGAGTGTAAATGCTTTAACGCATTCATTTATAATCATCCCTTCCATTATTATTATAAGCAAAATCCAAAGAATACATTGATGAAATAATAATGAAATATAATTAAAAATAAAATTTTTGCGTaaaaaaaatacaa                                                                                                             >  1:1451705/1‑133 (MQ=255)
          tAAATGCTTTAACGCATTCATTTATAATCATCCCTTCCATTATTATTATAAGCAAAATCCAAAGAATACATTGATGAAATAATAATGAAATATAATTAAAAATAAAATTTTTGCGTAAAAAAATACCACAGGCATTaaa                                                                                                   >  1:431052/1‑139 (MQ=255)
            aaTGCTTTAACGCATTCATTTATAATCATCCCTTCCATTATTATTATAAGCAAAATCCAAAGAATACATTGATGAAATAATAATGAAATATAATTAAAAATAAAATTTTTGCGTAAAAAAAATACCACAGGCATTaaaa                                                                                                  >  2:43168/1‑139 (MQ=255)
                  ttAACGCATTCATTTATAATCATCCCTTCCATTATTATTATAAGCAAAATCCAAAGAATACATTGATGAAATAATAATGAAATATAATTAAAAATAAAATTTTTGCGTAAAAAAAATACCACAGGCATTAAAAAATCAt                                                                                            >  2:291688/1‑139 (MQ=255)
                     aCGCATTCATTTATAATCATCCCTTCCATTATTATTATAAGCAAAATCCAAAGAATACATTGATGAAATAATAATGAAATATAATTAAAAATAAAATTTTTGCGTAAAAAAAATACCACAGGCATTAAAAAATCATgag                                                                                         <  1:291688/139‑1 (MQ=255)
                                  taATCATCCCTTCCATTATTATTATAAGCAAAATCCAAAGAATACATTGATGAAATAATAATGAAATATAATTAAAAATAAAATTTTTGCGTAAAAAAAATACCACAGGCATTAAAAAATCATGAGATGATTAAAatat                                                                            <  1:951140/139‑1 (MQ=255)
                                       atcCCTTCCATTATTATTATAAGCAAAATCCAAAGAATACATTGATGAAATAATAATGAAATATAATTAAAAATAAAATTTTTGCGTAAAAAAAATACCACAGGCATTAAAAAATCATGAGATGATTa                                                                                  <  1:102205/128‑1 (MQ=255)
                                       atcCCTTCCATTATTATTATAAGCAAAATCCAAAGAATACATTGATGAAATAATAATGAAATATAATTAAAAATAAAATTTTTGCGTAAAAAAAATACCACAGGCATTAAAAAATCATGAGATGATTa                                                                                  >  2:102205/1‑128 (MQ=255)
                                         cccTTCCATTATTATTATAAGCAAAATCCAAAGAATACATTGATGAAATAATAATGAAATATAATTAAAAATAAAATTTTTGCGTAAAAAAAATACCACAGGCATTAAAAAATCATGAGATGATTAAAATATTACAAtt                                                                     <  2:1451705/139‑1 (MQ=255)
                                               cATTATTATTATAAGCAAAATCCAAAGAATACATTGATGAAATAATAATGAAATATAATTAAAAATAAAATTTTTGCGTAAAAAAAATACCACAGGCATTAAAAAATCATGAGATGATTAAAATATTACAATTAGATta                                                               >  2:446108/1‑139 (MQ=255)
                                                      attatAAGCAAAATCCAAAGAATACATTGATGAAATAATAATGAAATATAATTAAAAATAAAATTTTTGCGTAAAAAAAATACCACAGGCATTAAAAAATCATgag                                                                                         <  1:40542/106‑1 (MQ=255)
                                                      attatAAGCAAAATCCAAAGAATACATTGATGAAATAATAATGAAATATAATTAAAAATAAAATTTTTGCGTAAAAAAAATACCACAGGCATTAAAAAATCATgag                                                                                         >  2:40542/1‑106 (MQ=255)
                                                            aGCAAAATCCAAAGAATACATTGATGAAATAATAATGAAATATAATTAAAAATAAAATTTTTGCGTAAAAAAAATACCACAGGCATTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCAtt                                                  >  2:1493832/1‑139 (MQ=255)
                                                                                    tgaAATAATAATGAAATATAATTAAAAATAAAATTTTTGCGTAAAAAAAATACCACAGGCATTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAg                          >  1:1428167/1‑139 (MQ=255)
                                                                                               tGAAATATAATTAAAAATAAAATTTTTGCGTAAAAAAAATACCACAGGCATTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGGGTTATAAGTCCTCAGGGCt               <  2:694385/139‑1 (MQ=255)
                                                                                                       aaTTAAAAATAAAATTTTTGCGTAAAAAAAATACCACAGGCATTAAAAAATCAAGAAATTATTAAAAAAATAAAATTAAATTATATTCAAATCATTAAACTTGAACCAAGGAATTAAAAAGCCTCAGGGCTTGCCaaaa       >  1:1252572/1‑139 (MQ=255)
                                                                                                         ttAAAAATAAAATTTTTGCGTAAAAAAAATACCACAGGCATTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACa     <  1:43168/139‑1 (MQ=255)
                                                                                                           aaaaaTAAAATTTTTGCGTAAAAAAATACCACAGGCATTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGaa  <  2:431052/139‑1 (MQ=255)
                                                                                                                              |                                                                                                                        
ATCCTGAGTGTAAATGCTTTAACGCATTCATTTATAATCATCCCTTCCATTATTATTATAAGCAAAATCCAAAGAATACATTGATGAAATAATAATGAAATATAATTAAAAATAAAATTTTTGCGTAAAAAAATACCACAGGCATTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAA  >  NC_000913/2995099‑2995344

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