Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,173,361 Δ2 bp 100% pseudogene (917‑918/1358 nt) gatC ← galactitol‑specific PTS enzyme IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,173,3610C.81.8% 22.1 / 2.9 11pseudogene (918/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base . (5/4);  total (5/6)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.
*NC_0009132,173,3620C.81.8% 24.2 / 0.8 11pseudogene (917/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (5/4);  minor base G (0/1);  total (5/6)
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
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCA  >  NC_000913/2173225‑2173491
                                                                                                                                        ||                                                                                                                                 
ccACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCTAAACGGTAGCATCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgccgt                                                                                                                              <  1:237824/143‑6 (MQ=255)
                          cACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCaa                                                                                                  <  2:123050/143‑1 (MQ=255)
                           gcGAGGAAGGCGATGGTGGCGAGATCGCCAAACGGCGGCACCTGATTACCGGGCGCACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAGCAGGCTTGCCGTTACCA‑‑CCGCCGTATGTCCTAGCAGCAACGCCGGATCAAg                                                                                                 <  2:125188/142‑1 (MQ=255)
                                             gCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACACACAGCAATTAAAATAGTGAGAGGAATAAAAAACAGGATTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGACGGAGCAAGGCCAATCAGGAACTCCTg                                                                               >  1:40774/1‑143 (MQ=255)
                                                                     tGATTCCCCGGCACACAAACAGCAATTAAAATGGTGAGTGTGATAAAATTCAGTTTTTCTTTTAC‑‑CTCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACg                                                       <  1:410030/143‑1 (MQ=255)
                                                                     tGATTACCCGGCACACAAACAGCAATTCAAATGGTGAGCGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGACGAATAACCCAAAAGCAACGCCGGAACAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACg                                                       >  1:286225/1‑143 (MQ=255)
                                                                     tGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGTTCCCA‑‑CCTCCGTATGTCCCAGCAGCAACGGCGGATCAAGGCAAATAAGGAACTCCTAACCGCCGAACTTCGACTGTAAAca                                                       >  1:299439/1‑142 (MQ=255)
                                                                           cccGGCACACAAACAGCAATTAAAATGGGGAGTGGGATAAAAATCAGGAGAGCAGAAACCA‑‑CCGAAGTATGGCCCAGCAGCAACGCGGGATCAAGGCGACTCAAGCAATGCTCGCCGCCGAACCGACCCTGTAAAcacgtacg                                                 >  2:347517/1‑136 (MQ=255)
                                                                                cccacaAACTGCAATTAAAATGTTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCt                                            <  2:163795/141‑1 (MQ=255)
                                                                                                        tcgggggATAAAAAGCAGGCTTGGCGATAC‑‑CGCCGCCGTATGTCCCAGCAGCAAGGACGGCGCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccat                    <  2:49751/139‑1 (MQ=255)
                                                                                                                    tgtCAGGCGTGCCGATGCCA‑‑CCGCCGTAGGTCCCAGCAGCAACGGCGGATCGATGTGAATCAGGAGCTCCGGACCGCCGAAGTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGGTAAACCATCCATGATGGGtt        <  2:312975/141‑1 (MQ=255)
                                                                                                                          gCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGATAAACCATCCATGATGGGTTTAATCa  >  2:166584/1‑143 (MQ=255)
                                                                                                                                        ||                                                                                                                                 
CCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCA  >  NC_000913/2173225‑2173491

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