Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007779 502,653 C→T intergenic (‑191/+47) ybaL ← / ← fsr transporter with NAD(P)‑binding Rossmann‑fold domain/fosmidomycin efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_007779502,6530CT97.6% 161.2 / ‑5.8 42intergenic (‑191/+47)ybaL/fsrtransporter with NAD(P)‑binding Rossmann‑fold domain/fosmidomycin efflux system
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (26/15);  minor base A (1/0);  total (27/15)
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

AATAATTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATA  >  NC_007779/502559‑502749
                                                                                              |                                                                                                
aataatTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGt                                                                                            <  2:109632/101‑1 (MQ=255)
aataatTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGt                                                                                            <  1:1295127/101‑1 (MQ=255)
            aTTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAg                                                                                >  1:637104/1‑101 (MQ=255)
             ttttGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGc                                                                               >  1:691890/1‑101 (MQ=255)
              tttGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGcc                                                                              >  1:1154364/1‑101 (MQ=255)
              tttGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGcc                                                                              >  1:836129/1‑101 (MQ=255)
                 gCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAAGCGGATTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCaaa                                                                           >  1:22570/1‑101 (MQ=255)
                   aGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGt                                                                         >  1:970276/1‑101 (MQ=255)
                     aaaaaTGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGttt                                                                       <  1:271761/101‑1 (MQ=255)
                       aaaTGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGc                                                                     >  1:1098308/1‑101 (MQ=255)
                              aTTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGc                                                              >  4:57731/1‑101 (MQ=255)
                                 aGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTAAGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGc                                                           >  1:871498/1‑101 (MQ=255)
                                     gCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATg                                                       <  3:293120/101‑1 (MQ=255)
                                      cAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGa                                                      <  3:879260/101‑1 (MQ=255)
                                      cAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGa                                                      <  2:144037/101‑1 (MQ=255)
                                      cAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGa                                                      <  2:142673/101‑1 (MQ=255)
                                         aaaaaGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAAt                                                   >  1:1513310/1‑101 (MQ=255)
                                           aaaGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCa                                                 <  1:665929/101‑1 (MQ=255)
                                            aaGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAg                                                >  3:522369/1‑101 (MQ=255)
                                            aaGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAg                                                >  3:531389/1‑101 (MQ=255)
                                            aaGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAg                                                >  1:83432/1‑101 (MQ=255)
                                             aGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGt                                               >  3:973122/1‑101 (MQ=255)
                                                  ggATCGGCATCCAGCCGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTa                                          >  4:8877/1‑101 (MQ=255)
                                                  ggATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTa                                          <  3:643325/101‑1 (MQ=255)
                                                   gATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTAt                                         >  1:175139/1‑101 (MQ=255)
                                                         aaacacgAAGGAGGTTTACCAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGt                                   >  4:160395/7‑101 (MQ=255)
                                                         aaaaCAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGt                                   >  1:572335/1‑101 (MQ=255)
                                                                aaGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATcagg                            <  3:465555/101‑1 (MQ=255)
                                                                  ggCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATcaggca                          <  1:754281/101‑1 (MQ=255)
                                                                     ggTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGa                       >  3:355667/1‑101 (MQ=255)
                                                                      gTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGaa                      >  3:1558254/1‑101 (MQ=255)
                                                                        ttAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAata                    <  3:1108708/101‑1 (MQ=255)
                                                                             agagGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTc               <  1:558479/101‑1 (MQ=255)
                                                                              gagGAAAGGCTTATGGAGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGCAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCa              >  1:517461/1‑101 (MQ=255)
                                                                               agGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCaa             <  3:698805/101‑1 (MQ=255)
                                                                               agGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCaa             >  2:134983/1‑101 (MQ=255)
                                                                                   aaGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATc         >  1:764584/1‑101 (MQ=255)
                                                                                         tATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAAt   >  3:647513/1‑101 (MQ=255)
                                                                                         tATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAAt   >  3:129328/1‑101 (MQ=255)
                                                                                          aTGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATa  <  3:82319/101‑1 (MQ=255)
                                                                                          aTGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATa  >  3:971134/1‑101 (MQ=255)
                                                                                          aTGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATa  >  1:334168/1‑101 (MQ=255)
                                                                                              |                                                                                                
AATAATTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATA  >  NC_007779/502559‑502749

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