Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 2,131,466 (A)9→8 intergenic (‑125/+175) gatY ← / ← fbaB D‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,131,4580A.95.2% 84.7 / 0.2 22intergenic (‑117/+183)gatY/fbaBD‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (10/10);  minor base T (1/0);  total (12/10)
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

ATGATTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGC  >  REL606/2131385‑2131522
                                                                         |                                                                
atgaTTTGTTTTTAAATGTCGAGATAACGATGACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaa                                                                >  1:469763/1‑76 (MQ=255)
atgaTTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaa                                                                <  2:706144/76‑1 (MQ=255)
        tttttAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCaaaaaaaacgg                                                        >  1:24480/1‑73 (MQ=255)
            tAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCGAAAAAAAACGGCCt                                                    >  2:746957/1‑76 (MQ=255)
            tAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTc                                                   <  1:1078507/76‑1 (MQ=255)
            tAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTc                                                   >  1:482298/1‑76 (MQ=255)
            tAAATGTCGAGAGAACGATCACAAAAACTGCAATATGAAAATTAATAGAGGAATGAATGGCAAAAAAAACGGCCTc                                                   <  2:1074773/76‑1 (MQ=255)
             aaaTGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTcc                                                  >  1:414725/1‑76 (MQ=255)
              aaTGTCGAGATAACGATCGCAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTccc                                                 >  1:459874/1‑76 (MQ=255)
              aaTGTCGAGATAACGATCACAAAACCGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTccc                                                 >  2:220600/1‑76 (MQ=255)
                 gTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGa                                              >  2:403460/1‑76 (MQ=255)
                   cGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGANGAATGATAGGCAAAAAAAACGGCCTCCCGGAgg                                            >  1:958208/1‑76 (MQ=255)
                      naTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGaa                                         <  1:1042729/75‑1 (MQ=255)
                      gaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGaa                                         <  2:1114968/76‑1 (MQ=255)
                      gaTAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGaa                                         <  2:415522/76‑1 (MQ=255)
                          aCGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCg                                     >  2:22614/1‑76 (MQ=255)
                             aTCACAAAAACGACAATATGCAAATCAATAGAGGATTGAAAGGTAACAAAAACGTCCTCCCGGAGGGAGGCCGcca                                  >  1:129039/1‑74 (MQ=255)
                             aTCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGcaa                                  <  2:815047/76‑1 (MQ=255)
                                  aaaaaCGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAg                             >  1:106944/1‑76 (MQ=255)
                                   aaaaCGACAATATGAAAATAAATAGAGAAAGAAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGt                            <  2:638147/76‑1 (MQ=255)
                                         aCAATATAAAAATTAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCTGAGGGAAGCCGCAACAAAGTGCGCGt                      <  1:1055441/76‑1 (MQ=255)
                                          aaaTATAAAAATAAATAGNGAAATGAAAGGAAATAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTg                     <  1:937601/75‑1 (MQ=255)
                                                aaaaaTTAATAGAGGAATAAAAGGAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTa               <  2:449702/75‑1 (MQ=255)
                                                     ttAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAg          <  2:1106534/76‑1 (MQ=255)
                                                      tAATAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAgg         >  1:342449/1‑76 (MQ=255)
                                                        aTAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGAt       >  2:362181/1‑76 (MQ=255)
                                                        aTAGAGGAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGAt       <  2:32733/76‑1 (MQ=255)
                                                        aTAGAGGAATAAAAGGAAAAAAAACCGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGAt       <  2:1159098/76‑1 (MQ=255)
                                                             ggAATGAAAGGCAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATgcggc  <  2:1134432/76‑1 (MQ=255)
                                                                         |                                                                
ATGATTTGTTTTTAAATGTCGAGATAACGATCACAAAAACGACAATATGAAAATTAATAGAGGAATGAAAGGCAAAAAAAAACGGCCTCCCGGAGGGAAGCCGCAACAAAGTGCGCGTGGTTTTATGCAGGATGCGGC  >  REL606/2131385‑2131522

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