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.94.4% 67.2 / ‑2.5 18pseudogene (918/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (7/10);  total (7/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01
Rejected as polymorphism: E-value score below prediction cutoff.
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.94.4% 65.9 / ‑2.7 18pseudogene (917/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (7/10);  minor base T (0/1);  total (7/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGT  >  NC_000913/2173238‑2173496
                                                                                                                           ||                                                                                                                                      
cggcgggcATCGCTACGAAGACGCCAATGGCGGCAAGAACGCCAAACGGCGGCACCTGAGTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCagcagc                                                                                                                    <  1:361322/136‑1 (MQ=255)
  aCCGCCATCGCCACGAAGAATCCGATGGTGGCAAGATCGCCAAACGGAAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCaa                                                                                                                  <  1:273507/143‑1 (MQ=255)
             cACGAAGAAGCCGATGGTGGCACGATCGCCTACCGGCTGCACCTGATTACCCAGCACACAAACAACAATAAAAATGGAGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCac                                                                                                       >  1:12544/1‑142 (MQ=255)
             cACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCaa                                                                                                       >  2:187095/1‑143 (MQ=255)
                  agaagCCGATGGTGGCAAGATCGCTAAACGGCAGTACTTTTTTACTTGGCTCTCATACATCATTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCa                                                                                                  <  1:168703/143‑1 (MQ=255)
                          aTGGTGGCAAGCTCGCCGAATGGCAGCGCCTGATTACCCAGCACACAAACATCAATTAAAATGGTGAGTGGTATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGaa                                                                                          <  2:324751/143‑1 (MQ=255)
                                  aaGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATAACAAAGCAACAAAACCGAATCAAGGCCAATCAGGAACTCCTGAc                                                                                  >  1:490018/1‑143 (MQ=255)
                                    taTCGCAAAACGGCAGCACCTGATTACACGGCACACAAACAGCAATGAAAAGGGGGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccg                                                                                <  1:116356/142‑1 (MQ=255)
                                            agcGGCATCACCTGGTTGCCCGGCACACAAGCAGGAGTTAAAATGGTGATTGGGATAAAAATCTGTCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCACCGCCGCGGCAATGCCAATCATGCACTCCTGACCGCCGAACtt                                                                        <  2:118293/141‑1 (MQ=255)
                                                    catcTGATTACCCTGCACACAAACAGCAAGTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTa                                                                <  1:408938/140‑1 (MQ=255)
                                                                              cAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGCACTGCATCCTTCTCCTGTAAACAAATAAGAGCATACTTaaagaaa                                      >  1:99651/1‑137 (MQ=255)
                                                                               cagtAAAATGGGGAGTGGTATAAAAATCATGCTTGCTGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCCAGGCAAATCAGGAACTCCTGACCGCCGAACTTCGTCTGTAAACTATTACGAGCCTGCTTAGCGATgg                                     <  2:210298/140‑1 (MQ=255)
                                                                               aaTTAAAATGGTGAGTGGGATAAAAATCAGGATTGTCTATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCGGCTGGCCGATgg                                     >  2:65948/1‑143 (MQ=255)
                                                                                    aaaTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACAAGCCTGCTGAGAGATGGGCGtt                                >  2:160376/1‑143 (MQ=255)
                                                                                      aTGGTGAGTGTGATAAAAATCAGGCTTGCCGATAC‑‑CTCCGCCGTATGTCCTAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCGGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaa                              <  1:360663/143‑1 (MQ=255)
                                                                                                  aTAAAAAGTAGGCTTGCCGATACCA‑‑CCGACGTATGTGCCAGCAGCAACGCCGGAGCAAGGCCAATGAGGGACGCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgat                  <  1:259804/143‑1 (MQ=255)
                                                                                                          cAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTaa          <  1:187095/143‑1 (MQ=255)
                                                                                                                  ccGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTCACCGACTCCTAGCCTGCTTAGCGATGGGCGTAAAACCAAACATGATAAGTTTAATCACACGt  >  1:13194/1‑143 (MQ=255)
                                                                                                                           ||                                                                                                                                      
CGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGT  >  NC_000913/2173238‑2173496

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