Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,922,908 G→T 12.5% intergenic (‑100/+59) menE ← / ← USA300HOU_1785 o‑succinylbenzoic acid‑‑CoA ligase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,922,9080GT12.5% 117.5 / 6.0 40intergenic (‑100/+59)menE/USA300HOU_1785o‑succinylbenzoic acid‑‑CoA ligase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (15/20);  new base T (2/3);  total (17/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

CCATTTTGTTGTGCTTGTTTATATAACCAAAAGTCCATGCGTTATTCCTCCAAAATCATTTACATTATAATTATAACGATTTTATGACATTCTAGCAGTGGTTATGTTTAAAAATATAAAAAAGTAGACGAATTGATGCATTGATATGATTGTTATAATGCTCAATACATATCGTTATATCATTCGTCTACTATTATCAGTTATTTTTATTTAATTTTAGTGTCATTCTGTCATTTTGATGTGGTGATTTACCCATTGTTGCCACATCATCTGCAAT  >  CP000730/1922771‑1923047
                                                                                                                                         |                                                                                                                                           
ccATTTTGTTGTGCTTGTTTATATAACCAAAAGTCCATGCGTTATTCCTCCAAAATCATTTACATTATAATTATAACGATTTTATGACATTCTAGCAGTGGTTATGTTTAAAAATATAAAAAAGTAGACGAATTGATGCAt                                                                                                                                          >  1:280938/1‑141 (MQ=255)
  aTTTTGTTGTGCTTGTTTATATAACCAAAAGTCCATGCGTTATTCCTCCAAAATCATTTACATTATAATTATAACGATTTTATGACATTCTAGCAGTGGTTATGTTTAAAAATATAAAAAAGTAGACGAATTGATGCATTg                                                                                                                                        >  2:207968/1‑141 (MQ=255)
    tttGTTGTGCTTGTTTATATAACCAAAAGTCCATGCGTTATTCCTCCAAAATCATTTACATTATAATTATAACGATTTTATGACATTCTAGCAGTGGTTATGTTTAAAAATATAAAAAAGTAGACGAATTGATGCATTGat                                                                                                                                      <  2:399715/141‑1 (MQ=255)
         tgtgCTTGTTTATATAACCAAAAGTCCATGCGTTATTCCTGCAAACTCATATAGATTATAATTATAACGATTTTATGACATTCTAGCAGTGGTTATGTTTAAAAATATAAAAAAGTAGACGAATTGATGCATTGATATGAt                                                                                                                                 <  2:284775/141‑1 (MQ=255)
         tgtgCTTGTTTATATAACCAAAAGTCCATGCGTTATTCCTCCAAAATCATTTACATTATAATTATAACGATTTTATGACATTCTAGCAGTGGTTATGTTTAAAAATATAAAAAAGTAGACGAATTGATGCATTGATATGAt                                                                                                                                 <  2:361402/141‑1 (MQ=255)
         tgtgCTTGTTTATAAAACGAAAGGGCCATGCGTTATTCCTCGAAAATCATTTACATTATAATTATAAGGATTTTATGACATTCTAGCAGTAGTTATGTTTAAAAATATAAAAAAGTAGACGAATTGATGCATTGATATGAt                                                                                                                                 <  2:34022/141‑1 (MQ=255)
            gCTTGTTTATATAACCAAAAGTCCATGCGTTATTCCTCCAAAATCATTTACATTATAATTATAACGATTTTATGACATTCTAGCAGTGGTTATGTTTAAAAATATAAAAAAGAAGACGAATTGATGCATTGATATGATTGt                                                                                                                              >  1:396876/1‑141 (MQ=255)
                ttttatttaACCAAAAGTCCATGCGTTATTCCTCCAAAATCATTTACATTATAATTATAACAATTTTAGGACATTCTAGCAGTGGTTAAGTTTAAAAATATAAAAAAGTAGACAAATTGATGCATTGATATGATTGTtata                                                                                                                          <  2:61277/140‑1 (MQ=255)
                    atataACCAAAAGTCCATGCGTAATTCCTCCAAAATCATTTACATTATAATTATAACGATTTTATGACATTCTAGCAGTGGTTATGTTTAAAAAGATAAAAAAGTAGACGAATTGATTCATTGATATGATTGTTATAATGc                                                                                                                      <  1:379021/141‑1 (MQ=255)
                      ataACCAAAAGTCCATGCGTTATTCCTCCAAAATCATTTACATTATAATTATAACGATTTTATGACATTCTAGCAGTGGTTATGTTTAAAAATATAAAAAAGTAGACGAATTGATGCATTGATATGATTGTTATAATGCTc                                                                                                                    >  1:76100/1‑141 (MQ=255)
                                   cATGCGTTATTCCTCCAAAATCATTTACATTATAATTATAACGATTTTATGACATTCTAGCAGTGGTTATGTTTAAAAATATAAAAAAGTAGACGAATTGATGCATTGATATGATTGTTATAATGCTCAATACATATCGtt                                                                                                       >  2:249325/1‑141 (MQ=255)
                                   cATGCGTTATTCCTCCAAAATCATTTACATTATAATTATAACGATTTTATGACATTCTAGCAGTGGTTATGTTTAAAAATATAAAAAAGTAGACGAATTGATGCATTGATATGATTGTTATAATGCTCAATACATATCGtt                                                                                                       >  1:424613/1‑141 (MQ=255)
                                   cATGCGTTATTCCTCCAAAATCATTTACATTATAATTATAACGATTTTATGACATTCTAGCAGTGGTTATGTTTAAAAATATAAAAAAGTAGACGAATTGATGCATTGATATGATTGTTATAATGCTCAATACATATCGtt                                                                                                       <  1:189604/141‑1 (MQ=255)
                                    aTGCGTTATTCCTCCAAAATCATTTACATTATAATTATAACGATTTTATGACATTCTAGCAGTGGTTATGTTTAAAAATATAAAAAAGTAGACGAATTGATTCATTGATATGATTGTTATAATGCTCAATACATATCGTta                                                                                                      <  1:449551/141‑1 (MQ=255)
                                    aTGCGTTATTCCTCCAAAATCATTTACATTATAATTATAACGATTTTATGACATTCTAGCAGTGGTTATGTTTAAAAATATAAAAAAGTAGACGAATTGATTCATTGATATGATTGTTATAATGCTCAATACATATCGTta                                                                                                      <  1:449548/141‑1 (MQ=255)
                                              cctccAAAATCATTTACATTATAATTATAACGATTTTATGACATTCTAGCAGTGGTTATGTTTAAAAATATAAAAAAGTAGACGAATTGATGCATTGATATGATTGTTATAATGCTCAATACATATCGTTATATCATTCgg                                                                                            >  1:210741/1‑140 (MQ=255)
                                                 ccAAAATCATTTACATTATAATTATAACGATTTTATGACATTCTAGCAGTGGTTATGTTTAAAAATATAAAAAAGTAGACGAATTGATGCATTGATATGATTGTTATAATGCTCAATACATATCGTTATATCATTCGTcta                                                                                         >  2:116818/1‑141 (MQ=255)
                                                       tCATTTACATTATAATTATAACGATTTTATGACATTCTAGCAGTGGTTATGTTTAAAAATATAAAAAAGTAGACGAATTGATGCATTGATATGATTGTTATAATGCTCAATACATATCGTTATATCATTCGTCTACtatta                                                                                   >  2:291261/1‑141 (MQ=255)
                                                        cATTTACATTATAATTATAACGATTTTATGACATTCTAGCAGTGGTTATGTTTAAAAATATAAAAAAGTAGACGAATTGATGCATTGATATGATTGTTATAATGCTCAATACATATCGTTATATCATTCGTCTACtattat                                                                                  <  2:70219/141‑1 (MQ=255)
                                                        cATTTACATTATAATTATAACGATTTTATGACATTCTAGCAGTGGTTATGTTTAAAAATATAAAAAAGTAGACGAATTGATGCATTGATATGATTGTTATAATGCTCAATACATATCGTTATATCATTCGTCTACtattat                                                                                  >  2:16173/1‑141 (MQ=255)
                                                             aCATTATAATTATAACGATTTTATGACATTCTAGCAGTGGTTATGTTTAAAAATATAAAAAAGTAGACGAATTGATGCATTGATATGATTGTTATAATGCTCAATACATATCGTTATATCATTCGTCTACTATTATCAGtt                                                                             <  2:104361/141‑1 (MQ=255)
                                                                  ataattataaCGATTTTATGACATTCTAGCAGTGGTTATGTTTAAAAATATAAAAAAGTAGACGAATTGATGCATTGATATGATTGTTATAATGCTCAATACATATCGTTATATCATTCGTCTACTATTATCAGttatttt                                                                        <  1:134832/141‑1 (MQ=255)
                                                                         taaCGATTTTATGACATTCTAGCAGTGGTTATGTTTAAAAATATAAAAAAGTAGACGAATTGATGCATTGATATGATTGTTATAATGCTCAATACATATCGTTATATCATTCGTCTACTATTATCAGTTATTTTTATTTaa                                                                 >  1:151279/1‑141 (MQ=255)
                                                                           aCGATTTTATGACATTCTAGCAGTGGTTATGTTTAAAAATATAAAAAAGTAGACGAATTGATGAATTGATATGATTGTTATAATGCTCAATACATATCGTTATATCATTCGTCTACTATTATCAGTTATTTTTATTTAAtt                                                               <  2:298151/141‑1 (MQ=255)
                                                                             gATTTTATGACATTCTAGCAGTGGTTATGTTTAAAAATATAAAAAAGTAGACGAATTGATGCATTGATATGATTGTTATAATGCTCAATACATATCGTTATATCATTCGTCTACTATTATCAGTTATTTTTATTTAAtttt                                                             <  1:310665/141‑1 (MQ=255)
                                                                              aTTTTATGACATTCTAGCAGTGGTTATGTTTAAAAATATAAAAAAGTAGACGAATTGATGCATTGATATGATTGTTATAATGCTCAATACATATCGTTATATCATTCGTCTACTATTATCAGTTATTTTTATTTAATTTTa                                                            <  1:81962/141‑1 (MQ=255)
                                                                                      aCATTCTAGCAGTGGTTATGTTTAAAAATATAAAAAAGTAGACGAATTGATGCATTGATATGATTGTTATAATGCTCAATACATATCGTTATATCATTCGTCTACTATTATCAGTTATTTTTATTTAATTTTAGTGTCAtt                                                    >  2:402921/1‑141 (MQ=255)
                                                                                                 gTGGTTATGTTTAAAAATATAAAAAAGTAGACGAATTGATGCATTGATATGATTGTTATAATGCTCAATACATATCGTTATATCATTCGTCTACTATTATCAGTTATTTTTATTTAATTTTAGTGTCATTCTGTCATTTTg                                         >  1:173979/1‑141 (MQ=255)
                                                                                                       aagTTCAAAAATATAAAATAGTAGTCGAATTGATGCATTGATATGATTATTATAATGCTCAATACATCTCGTTATATCATTCGTCTCCTATTATTAGTTATTCTTATTTAATTTTAGTGTCATCCTGTCATTTTGATgtgg                                   <  2:105310/139‑1 (MQ=255)
                                                                                                            tAAAAATATAAAAAAGTAGACGAATTGATGCATTGATATGATTGTTATAATGCTCAATACATATCGTTATATCATTCGTCTACTATTATCAGTTATTTTTATTTAATTTTAGTGTCATTCTGTCATTTTGATGTGGTGAtt                              <  2:151279/141‑1 (MQ=255)
                                                                                                              aaaaTATAAAAAAGTAGACGAATTGATGCATTGATATGATTGTTATAATGCTCAATACATATCGTTATATCATTCGTCTACTATTATCAGTTATTTTTATTTAATTTTAGTGTCATTCTGTCATTTTGATGTGGTGATTTa                            <  2:280938/141‑1 (MQ=255)
                                                                                                                aaTATAAAAAAGTAGACGAATTGATGCATTGATATGATTGTTATAATGCTCAATACATATCGTTATATCATTCGTCTACTATTATCAGTTATTTTTATTTAATTTTAGTGTCATTCTGTCATTTTGATGTGGTGATTTAcc                          <  2:7335/141‑1 (MQ=255)
                                                                                                                aaTATAAAAAAGTAGACGAATTGATGCATTGATATGATTGTTATAATGCTCAATACATATCGTTATATCATTCGTCTACTATTATCAGTTATTTTTATTTAATTTTAGTGTCATTCTGTCATTTTGATGTGGTGATTTAcc                          <  1:427126/141‑1 (MQ=255)
                                                                                                                   ataAAAAAGTAGACGAATTGATTCATTGATATGATTGTTATAATGCTCAATACATGTCGTTATATCATTCGTCTACTATTATCAGTTATTTTTATTTAATTTTAGTGTCATTCTGTCATTTTGATGTGGTGATTTAcccgt                       >  1:43265/1‑139 (MQ=255)
                                                                                                                   ataAAAAAGTAGACGAATTGATTCATTGATATGATTGTTATAATGCTCAATACATATCGTTATATCATTCGTCTACTATTATCAGTTATTTTTATTTAATTTTAGTGTCATTCTGTCATTTTGATGTGGTGATTTAcccca                       >  1:43023/1‑139 (MQ=255)
                                                                                                                   ataAAAAAGTAGACGAATTGATGCATTGATATGATTGTTATAATGCTCAATACATATCGTTATATCATTCGTCTACTATTATCAGTTATTTTTATTTAATTTTAGTGTCATTCTGTCATTTTGATGTGGTGATTTACCCAt                       >  1:336991/1‑141 (MQ=255)
                                                                                                                           gTAGACGAATTGATGCATTGATATGATTGTTATAATGCTCAATACATATCGTTATATCATTCGTCTACTATTATCAGTTATTTTTATTTAATTTTAGTGTCATTCTGTCATTTTGATGTGGTGATTTACCCATTGTTGCca               >  2:217641/1‑141 (MQ=255)
                                                                                                                             aGACGAATTGATGCATTGATATGATTGTTATAATTCTCAATACATATCGTTATATCATTCGTCTACTATTATCAGTTATTTTTATTTAATTTTAGTGTCATTCTGTCATTTTGATGTGGTGATTTACCCATTGTTGCcaca             <  1:128374/141‑1 (MQ=255)
                                                                                                                                     tGATGCATTGATATGATTGTTATAATGCTCAATACATATCGTTATATCATTCGTCTACTATTATCAGTTATTTTTATTTAATTTTAGTGTCATTCTGTCATTTTGATGTGGTGATTTACCCATTGTTGCCACATCATCTGc     <  1:193502/141‑1 (MQ=255)
                                                                                                                                        tGCATTGATATGATTGTTATAATGCTCAATACATATCGTTATATCATTCGTCTACTATTATCAGTTATTTTTATTTAATTTTAGTGTCATTCTGTCATTTTGATGTGGTGATTTACCCATTGTTGCCACATCATCTGCAAt  <  1:207968/141‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                           
CCATTTTGTTGTGCTTGTTTATATAACCAAAAGTCCATGCGTTATTCCTCCAAAATCATTTACATTATAATTATAACGATTTTATGACATTCTAGCAGTGGTTATGTTTAAAAATATAAAAAAGTAGACGAATTGATGCATTGATATGATTGTTATAATGCTCAATACATATCGTTATATCATTCGTCTACTATTATCAGTTATTTTTATTTAATTTTAGTGTCATTCTGTCATTTTGATGTGGTGATTTACCCATTGTTGCCACATCATCTGCAAT  >  CP000730/1922771‑1923047

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑