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.83.3% 41.9 / 9.0 18intergenic (+235/‑58)yccE/agpYccE family protein/bifunctional glucose‑1‑phosphatase/inositol phosphatase
Reads supporting (aligned to +/- strand):  ref base A (0/3);  new base . (8/7);  total (8/10)
Fisher's exact test for biased strand distribution p-value = 2.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

GGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAGAGGTGCCAGGAATGAACAAAACGCTAATCGCCGCAG  >  NZ_CP009273/1060894‑1061065
                                                                                         |                                                                                  
ggCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGa                                                                                    >  1:230132/1‑90 (MQ=255)
 gCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGaa                                                                                   >  1:973254/1‑90 (MQ=255)
      gTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGaaaaagt                                                                              >  2:144802/1‑88 (MQ=255)
              cAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTc                                                                     <  2:1011076/90‑1 (MQ=255)
                aaCCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTg                                                                   >  2:409700/1‑90 (MQ=255)
                 aCCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGt                                                                  <  1:307494/90‑1 (MQ=255)
                   cTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTca                                                                <  2:524802/90‑1 (MQ=255)
                      gCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCacac                                                             <  1:760910/90‑1 (MQ=255)
                            ttttGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAAGTCATATTTCTGTCACACTCttt                                                        <  2:973254/90‑1 (MQ=255)
                                   gATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGtgattg                                                <  1:409700/90‑1 (MQ=255)
                                   gATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGtgattg                                                <  2:812192/90‑1 (MQ=255)
                                    aTGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGtgattga                                               >  1:908335/1‑90 (MQ=255)
                                       tatataGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATa                                             <  2:538537/90‑1 (MQ=255)
                                             gAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCata                                                                         >  1:1054626/1‑55 (MQ=255)
                                             gAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCata                                                                         <  2:1054626/55‑1 (MQ=255)
                                                aaaCAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAGAgg                                   >  1:496360/1‑90 (MQ=255)
                                                    aGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAGAGGTGcc                               >  1:1032700/1‑90 (MQ=255)
                                                     ggCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAGAGGTGcc                               <  2:984714/90‑1 (MQ=255)
                                                       cGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAGAGGTGCCAgg                            >  2:1031298/1‑90 (MQ=255)
                                                                       ttGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAGAGGTGCCAGGAATGAACAAAACGCTa            >  2:801418/1‑90 (MQ=255)
                                                                                 gTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAGAGGTGCCAGGAATGAACAAAACGCTAATCACCGCAg  >  1:969347/1‑90 (MQ=255)
                                                                                         |                                                                                  
GGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAGAGGTGCCAGGAATGAACAAAACGCTAATCGCCGCAG  >  NZ_CP009273/1060894‑1061065

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