Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,252,901 A→C intergenic (‑126/+297) psuK ← / ← fruA pseudouridine kinase/PTS fructose transporter subunit IIBC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,252,9010AC94.1% 59.1 / ‑3.0 17intergenic (‑126/+297)psuK/fruApseudouridine kinase/PTS fructose transporter subunit IIBC
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (13/3);  total (14/3)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.47e-01

ACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATT  >  NZ_CP009273/2252816‑2252970
                                                                                     |                                                                     
acaAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAAtt                                                                   >  2:317835/1‑90 (MQ=255)
           gATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCagaga                                                        >  1:124729/1‑90 (MQ=255)
             tACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAgg                                                      >  2:303575/1‑90 (MQ=255)
              aCTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGc                                                     >  1:415426/1‑90 (MQ=255)
              aCTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGc                                                     >  1:422168/1‑90 (MQ=255)
                    aaGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATg                                               >  1:149964/1‑90 (MQ=255)
                                 cATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGc                                  <  2:202100/90‑1 (MQ=255)
                                    gtgtTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAgg                               >  2:438070/1‑90 (MQ=255)
                                    gtgtTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAgg                               >  2:393674/1‑90 (MQ=255)
                                     tgtTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGc                              <  2:28688/90‑1 (MQ=255)
                                           tatTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACaa                        >  1:137413/1‑90 (MQ=255)
                                                     tgtgCTCATAGTATGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTa              >  2:264200/1‑90 (MQ=255)
                                                           cATAGTATGCATGTTTTACGTAAAGTCAAATACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGcc        >  2:305801/1‑90 (MQ=255)
                                                                 aTGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGGGCAGGCTCACAAAATGCGATTAAATGCCTGAAtt  >  2:225782/1‑90 (MQ=255)
                                                                 aTGCATGTTTTACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAAtt  >  1:261719/1‑90 (MQ=255)
                                                                           tACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGc                                  <  2:111184/48‑1 (MQ=255)
                                                                           tACGTAAAGTCAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGc                                  >  1:111184/1‑48 (MQ=255)
                                                                                     |                                                                     
ACAAGAATATAGATACTGCAAAGACGGTTTAAACATGTGTTCATATTTAAACATGTGCTCATAGTATGCATGTTTTACGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATT  >  NZ_CP009273/2252816‑2252970

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