Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,173,363 Δ2 bp pseudogene (915‑916/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.86.4% 71.2 / 5.4 22pseudogene (918/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base . (8/11);  total (8/14)
Fisher's exact test for biased strand distribution p-value = 2.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009132,173,3620C.90.5% 71.2 / 2.8 22pseudogene (917/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base . (8/11);  minor base G (0/2);  total (8/14)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAA  >  NC_000913/2173234‑2173488
                                                                                                                               ||                                                                                                                              
aCGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTTCCCGGCCCACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcag                                                                                                                <  1:232595/143‑1 (MQ=255)
aCGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGACAAACGGCAGCACCTGATTACCCGGCACACCAACAGCAAATAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcag                                                                                                                >  1:345151/1‑143 (MQ=255)
                ccACGCCGAAGCCGATGGGGGCAAGATCTCGAAGCGGGAGCACCTGATTACCCGGCACACAATCAGCAATTAAAATGGGGAGTTGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCa                                                                                                <  2:157214/143‑1 (MQ=255)
                ccACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCa                                                                                                >  1:552321/1‑143 (MQ=255)
                            cGATGGTGGCAAGATCGCCAGACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAgg                                                                                    <  2:466756/143‑1 (MQ=255)
                              aTGGTGTCAAGATCGCCAAACGGCAGCACCGGATGACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGaa                                                                                  <  1:206523/143‑1 (MQ=255)
                                     gaaGATCGCCAACCGGGAGCACCTGATTACCCGGCACACAAGCAGCAGTTCAAATGGGGAGTCGGATAAAACTCAGGCTTGCCGATAC‑‑CGCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGa                                                                           <  2:267652/142‑1 (MQ=255)
                                           cGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgccg                                                                     >  1:37302/1‑143 (MQ=255)
                                             gcgagcGGCAGCACCTGATGACCCGGGACACAAACAGCAATTCAAATGGTGAGTGGGATAAAAAGCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGaa                                                                   <  2:668020/138‑1 (MQ=255)
                                                    gcagcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAGCTTCGcc                                                            <  2:598304/143‑1 (MQ=255)
                                                       gcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAACAACGCCGGATCAAGGCCAATCAGGAACCCCTGACCGCCGCACCTCGCCTgc                                                         >  2:1543/1‑142 (MQ=255)
                                                           cTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATAAGGCTTGCGCATAGAA‑‑ACTCCGCATTGCCCAGCAGCAACGCCGGATCAAGGCCAAGCAGGAACTCCTGACCGCCGCCCTCCGCCTGTAGAc                                                     >  2:628433/1‑143 (MQ=255)
                                                                      gcccacaAACAGCAACTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGcc                                          <  1:377303/140‑1 (MQ=255)
                                                                          acagacagCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGTGCGTGCt                                      >  1:472672/1‑143 (MQ=255)
                                                                            aaaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTCCGCGCCTGCTTa                                    >  1:150180/1‑143 (MQ=255)
                                                                                   aaTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATgg                             <  1:280524/143‑1 (MQ=255)
                                                                                       aagaTGGTGAGTGGGATAAAAATCAGGCTTGCCGATAC‑‑CGCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGt                         <  2:649251/140‑1 (MQ=255)
                                                                                         aaagggtaGTGGGATAAAAAGCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAAGTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTa                       <  1:598260/136‑1 (MQ=255)
                                                                                         aaTGGTGAGCGGCATAACAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTa                       <  1:155889/143‑1 (MQ=255)
                                                                                           tGGTGAGCGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGTCCAATCTGGAACACTTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaaa                     <  1:683622/143‑1 (MQ=255)
                                                                                                           aaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATAGGCGTTAAACCATCCATGATGGGtt     >  1:35043/1‑143 (MQ=255)
                                                                                                              gtggCTTGCCGATAC‑‑CCCCGCCGTATGGCCCGGGAGCAGCGCCGGATCAAGGCCCAGCCGCAACTCTGGGCCGCCGAACTGCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTaa  <  2:252309/141‑1 (MQ=255)
                                                                                                                               ||                                                                                                                              
ACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAA  >  NC_000913/2173234‑2173488

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