Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 502,653 C→T intergenic (‑191/+47) ybaL ← / ← fsr putative transport protein/fosmidomycin resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913502,6530CT98.8% 238.9 / NA 81intergenic (‑191/+47)ybaL/fsrputative transport protein/fosmidomycin resistance protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (50/30);  minor base G (0/1);  total (50/31)
Fisher's exact test for biased strand distribution p-value = 3.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGC  >  NC_000913/502582‑502724
                                                                       |                                                                       
aaaTGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACt                                                                     >  2:1580706/1‑76 (MQ=255)
aaaTGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACt                                                                     <  2:3467774/76‑1 (MQ=255)
 aaTGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGCGGAAAGGCTTATGGTGACTg                                                                    >  2:4837519/1‑76 (MQ=255)
 aaTGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTg                                                                    <  1:2357077/76‑1 (MQ=255)
  aTGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGt                                                                   <  2:2297294/76‑1 (MQ=255)
  aTGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGt                                                                   >  1:3615435/1‑76 (MQ=255)
  aTGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGt                                                                   <  1:158738/76‑1 (MQ=255)
   tGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTa                                                                  <  2:1145209/76‑1 (MQ=255)
    ggAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTaa                                                                 >  2:1782836/1‑76 (MQ=255)
    ggAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTaa                                                                 <  2:2286933/76‑1 (MQ=255)
     gAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGTCTTATGGTGATGGTaaa                                                                >  1:296073/1‑76 (MQ=255)
     gAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTaaa                                                                >  1:1031253/1‑76 (MQ=255)
     gAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTaaa                                                                >  2:939696/1‑76 (MQ=255)
      aaTTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTaaaa                                                               >  1:1616437/1‑76 (MQ=255)
      aaTTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTaaaa                                                               >  1:1519822/1‑76 (MQ=255)
      aaTTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTaaaa                                                               >  1:4098813/1‑76 (MQ=255)
       aTTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAAc                                                              >  1:4687683/1‑76 (MQ=255)
         tAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGa                                                            >  1:2930115/1‑76 (MQ=255)
          aGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGCGGAAAGGCTTATGGTGACTGTAAAACGAt                                                           >  2:2854543/1‑76 (MQ=255)
          aGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGAt                                                           >  2:2828258/1‑76 (MQ=255)
          aGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGAt                                                           <  2:2479456/76‑1 (MQ=255)
          aGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGAt                                                           <  1:4349154/76‑1 (MQ=255)
           gCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATg                                                          >  1:2626702/1‑76 (MQ=255)
            cTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGc                                                         >  2:930911/1‑76 (MQ=255)
               cAAAAACAGCACGGATCGGAACACAGAAGGCGGTTTAACAGAGGCAAGGCTTATGGTGACTGTAAAACGATGCAGc                                                      >  2:2855745/1‑76 (MQ=255)
               cAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGc                                                      <  2:4470138/76‑1 (MQ=255)
               cAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGc                                                      >  2:2837349/1‑76 (MQ=255)
                aaaaaaaGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAgac                                                     >  1:2334016/1‑74 (MQ=255)
                  aaaaaGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCaa                                                   <  1:4604194/76‑1 (MQ=255)
                   aaaaGCACGGATCGGAAAACAGAAGGCGGTTTAACGGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCaaa                                                  >  2:4539571/1‑76 (MQ=255)
                   aaaaGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCaaa                                                  <  2:2937129/76‑1 (MQ=255)
                     aaGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGt                                                >  1:249731/1‑76 (MQ=255)
                     aaGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGt                                                >  2:691160/1‑76 (MQ=255)
                     aaGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGt                                                <  2:3473051/76‑1 (MQ=255)
                     aaGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGt                                                >  1:1338988/1‑76 (MQ=255)
                       gCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGttt                                              >  1:3644598/1‑76 (MQ=255)
                        cACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTg                                             >  1:1444746/1‑76 (MQ=255)
                        cACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTg                                             >  1:768503/1‑76 (MQ=255)
                         aCGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGc                                            >  1:1134187/1‑76 (MQ=255)
                         aCGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGc                                            <  1:635366/76‑1 (MQ=255)
                          cGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGCTTGcc                                           >  2:2623811/1‑76 (MQ=255)
                           ggATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCt                                          >  1:1122866/1‑76 (MQ=255)
                            gATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCtt                                         <  2:832332/76‑1 (MQ=255)
                             aTCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCttt                                        >  2:1132011/1‑76 (MQ=255)
                              tCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTg                                       >  1:2933716/1‑76 (MQ=255)
                               cGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTgg                                      <  2:4490800/76‑1 (MQ=255)
                               cGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTgg                                      >  2:3961319/1‑76 (MQ=255)
                                ggAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGc                                     <  1:4185467/76‑1 (MQ=255)
                                 gAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCt                                    >  1:1604741/1‑76 (MQ=255)
                                  aaaaCAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTg                                   <  2:3189558/76‑1 (MQ=255)
                                  aaaaCAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTg                                   <  2:2005383/76‑1 (MQ=255)
                                  aaaaCAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTCTGCCTTTGGCtc                                   >  2:1564942/1‑75 (MQ=255)
                                    aaCAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCa                                 >  1:2284340/1‑76 (MQ=255)
                                     aCAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCCGCCCAAGTGTGCCTGTCGCTGCaa                                >  1:3199551/1‑76 (MQ=255)
                                     aCAGAAGGCGGTTTAACAGAGGAAAGGATTATGGTGACTGTAAAACGGTGCAGCCAAAGCTTGCCTGTCGCTGCaa                                >  1:3163346/1‑76 (MQ=255)
                                         aaGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGaa                            >  2:3527266/1‑76 (MQ=255)
                                         aaGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGaa                            <  2:249731/76‑1 (MQ=255)
                                         aaGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGaa                            >  2:2638119/1‑76 (MQ=255)
                                           ggCGGTTTAACATAGGAAAGGCTTATGGTGACTGTAAAACTATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAAt                          <  1:253419/76‑1 (MQ=255)
                                             cGGTTTAGCAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCa                        <  2:2937116/76‑1 (MQ=255)
                                             cGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCa                        <  1:3917677/76‑1 (MQ=255)
                                              ggTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAg                       >  1:393041/1‑76 (MQ=255)
                                               gTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGt                      <  1:3648587/76‑1 (MQ=255)
                                               gTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGt                      <  2:3521201/76‑1 (MQ=255)
                                                tttAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTc                     <  1:790870/76‑1 (MQ=255)
                                                 ttAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGt                      >  1:2008938/1‑74 (MQ=255)
                                                 ttAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGt                      <  2:2008937/74‑1 (MQ=255)
                                                   aaCAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCttt                  >  1:996163/1‑76 (MQ=255)
                                                     cAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCCTTGGCTGCAATGAAATCAGTCTTTAt                >  2:223800/1‑76 (MQ=255)
                                                      agagGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGCCTTTATg               >  2:1824134/1‑76 (MQ=255)
                                                        agGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGcc             >  1:1423782/1‑76 (MQ=255)
                                                        agGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGcc             >  2:76507/1‑76 (MQ=255)
                                                         ggAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCg            <  1:768288/76‑1 (MQ=255)
                                                            aaGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGtt         >  1:3219023/1‑76 (MQ=255)
                                                             aGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTa        <  2:1733341/76‑1 (MQ=255)
                                                              ggCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTAt       <  2:199516/76‑1 (MQ=255)
                                                               gcctATGGTGACTGGAAAACGAGCCAGCCAAAGTGTGCCTTGGACTGCAATGAAACCCGTCTTTATGCCGGTTATc      >  2:759295/4‑76 (MQ=255)
                                                               gCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATc      <  2:4805337/76‑1 (MQ=255)
                                                               gCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATc      >  1:3640232/1‑76 (MQ=255)
                                                                 ttATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATcag    >  1:1670111/1‑76 (MQ=255)
                                                                  tATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATcagg   >  1:506973/1‑76 (MQ=255)
                                                                   aTGGGGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATcaggc  <  2:608501/76‑1 (MQ=255)
                                                                       |                                                                       
AAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGC  >  NC_000913/502582‑502724

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