Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,060,988 (A)6→5 intergenic (+240/‑53) yccE → / → agp YccE family protein/bifunctional glucose‑1‑phosphatase/inositol phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,060,9830A.81.2% 32.3 / 10.2 16intergenic (+235/‑58)yccE/agpYccE family protein/bifunctional glucose‑1‑phosphatase/inositol phosphatase
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base . (7/6);  total (9/7)
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

GTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAGAGGTGCCAGGAATGAACAAAACGC  >  NZ_CP009273/1060900‑1061053
                                                                                   |                                                                      
gTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGaaaaagt                                                                  >  1:795174/1‑88 (MQ=255)
gTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAAg                                                                  >  1:464396/1‑90 (MQ=255)
  ttGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGaaaaagtca                                                                >  2:297977/1‑86 (MQ=255)
    gCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCata                                                             >  1:322937/1‑90 (MQ=255)
        cAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAAGTCATAttt                                                          >  2:296843/1‑90 (MQ=255)
              ttGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTcac                                                   <  1:718141/90‑1 (MQ=255)
              ttGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTcac                                                   >  2:287788/1‑90 (MQ=255)
                gCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCacac                                                 <  1:11612/90‑1 (MQ=255)
                gCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCacac                                                 <  1:664836/90‑1 (MQ=255)
                 ccTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAAGTCATATTTCTGTCacac                                                 <  1:171172/90‑1 (MQ=255)
                             gATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGtgattg                                    <  1:297977/90‑1 (MQ=255)
                              aTGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGtgattga                                   >  1:594377/1‑90 (MQ=255)
                                gTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATa                                 <  1:499741/90‑1 (MQ=255)
                                gTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATa                                 <  2:709433/90‑1 (MQ=255)
                                    ataGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACaa                             >  1:160587/1‑90 (MQ=255)
                                       gAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAag                          >  2:827017/1‑90 (MQ=255)
                                          aaaCAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAGAgg                       >  2:586537/1‑90 (MQ=255)
                                                               ttttGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAGAGGTGCCAGGAATGAACAAAACGc  >  2:542880/1‑90 (MQ=255)
                                                                                   |                                                                      
GTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAGAGGTGCCAGGAATGAACAAAACGC  >  NZ_CP009273/1060900‑1061053

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