Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,173,363 Δ2 bp intergenic (‑1/+1) gatC ← / ← gatC pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC/pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,173,3610C.88.6% 121.3 / 7.7 35pseudogene (2/442 nt)gatCpseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC
Reads supporting (aligned to +/- strand):  ref base C (0/4);  new base . (14/17);  total (14/21)
Fisher's exact test for biased strand distribution p-value = 1.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009132,173,3620C.91.2% 123.0 / 5.1 35pseudogene (1/442 nt)gatCpseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base . (14/17);  total (14/21)
Fisher's exact test for biased strand distribution p-value = 2.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAG‑CAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGG  >  NC_000913/2173222‑2173498
                                                                                                                                           ||                                                                                                                                         
tttCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGCTTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgc                                                                                                                                         >  1:243376/1‑141 (MQ=255)
    cACGACACACGGCACCCGCCCTCGCCACGCAGACGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTCCCTGGCACACAAACAGCAATTAAAATGGTGAGTGGTATAAAAATCAGGCTTGCCGATaccaccgccgta                                                                                                                                     <  1:525622/143‑7 (MQ=255)
          gcccGGCGACCGCCACCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTccc                                                                                                                             <  2:433939/140‑1 (MQ=255)
             cGGCGCCCGCCATCGCCACGACGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagc                                                                                                                          <  2:733876/143‑1 (MQ=255)
              ggtgACCGCCTTCTCTTCGGAGGAGGCGATGGTGGCAAGGTCGCCAAACGGCAGCTCCGGATTGCGCGGTACACATGCAGCAGTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagca                                                                                                                         <  2:465119/140‑1 (MQ=255)
                   ccGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCGCACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAG‑GAAc                                                                                                                   >  1:724853/1‑143 (MQ=255)
                   ccGCCACAGCGACGAAGAAGCAGATGGTGGCAAGACCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAG‑CAAc                                                                                                                   <  1:934027/143‑1 (MQ=255)
                    cGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAG‑CAACg                                                                                                                  >  1:680086/1‑143 (MQ=255)
                            ccACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATCTCCCAGCAG‑CAACGCCGGATCa                                                                                                          >  1:22114/1‑143 (MQ=255)
                             cACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGGCCCAGCAC‑CAACGCCGGATcca                                                                                                         >  2:19704/1‑141 (MQ=255)
                             cACGAAGAAGCCGATGGCGGGAAGAGCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAG‑CAACGCCGGATCaa                                                                                                         <  1:636433/143‑1 (MQ=255)
                              aCGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCGGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAG‑CAACGCCGGATCAAg                                                                                                        <  1:322960/143‑1 (MQ=255)
                               cGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGACGATACCA‑‑CCGCCGTATGTCCCAGCCA‑AAAGATCGGATCAAgg                                                                                                       >  2:707337/1‑143 (MQ=255)
                                   gaagCCGAGGGTGGCAAGACCGCGAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGGGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAG‑CAACGCCGGATCAAGGCCaa                                                                                                   <  2:478532/143‑1 (MQ=255)
                                       ccGATCGGGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAG‑CAACGCCGGATCAAGGCCAATCAg                                                                                               <  2:165532/143‑1 (MQ=255)
                                         gATGGTGGCAAGATCGCCAAAAGGCAGCACCTGATTACCCGGCACACAAAAAGCAACAGTGATCGTACACAGAAAAAAAAACAGACATGCCGATACCA‑‑CCGCCGTATGGCCCAGCAG‑CAACGGCGGATCAAGGCCCATCAGGa                                                                                             >  2:679271/1‑143 (MQ=255)
                                                    gATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAG‑CAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccg                                                                                  <  1:758080/143‑1 (MQ=255)
                                                                          aTTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAG‑CAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGAc                                                            <  1:634782/143‑1 (MQ=255)
                                                                               ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAG‑CAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGACTGTAAAAGACTACGa                                                       >  2:484944/1‑143 (MQ=255)
                                                                                     cacaATCAGCAATTAAAATGGTGAGTGGGATAAAGATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAG‑CAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGc                                                 <  2:751350/143‑1 (MQ=255)
                                                                                        aaaCAGCAATTAAAATGGTGAGTGGGATAACAATCAGGCCTGCCGATAC‑‑CGCCGCCGTATGTCCCAGCAG‑CAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTa                                              <  2:569289/143‑1 (MQ=255)
                                                                                            aGCAATTAAAATGGGGAGTGGGATAAAAATCAGGCTTGCCGATAC‑‑CCGCG‑CGGATACGGTAGCCGTCTACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGa                                          <  1:441247/143‑1 (MQ=255)
                                                                                             gCAATCAAAATCGTGAGTGGGATAAACATCAGCCTTGCCGATAC‑‑CCCCGCCGGATGTCCCAGCAG‑CAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAt                                         <  1:428616/143‑1 (MQ=255)
                                                                                                 ttAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAG‑CAACGCCGGATCAAGGCCAATCAGGAACTCCGGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGc                                     <  2:680086/143‑1 (MQ=255)
                                                                                                   aaaaTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAG‑CAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGt                                   <  2:521576/143‑1 (MQ=255)
                                                                                                   aaaaGGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAG‑CAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGt                                   <  2:368605/143‑1 (MQ=255)
                                                                                                    aaaTGGGGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAG‑CAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGtt                                  <  2:504763/143‑1 (MQ=255)
                                                                                                          ggAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAG‑CAACGCCTGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAcca                            <  2:470803/142‑1 (MQ=255)
                                                                                                           gAGTGGGATAAAAATGAGGCTTGCCGATACCA‑‑CCGCGGTATGTCCCAGCAG‑CAACGCCGGATCAAGGCCAATGAGGAACTCCTGACCGCCGAACTTCGCCTGGAACAGACTACGAGCCTGCTTAGCGATCGGCGCTAAAcaaa                           >  2:118985/1‑140 (MQ=255)
                                                                                                                  aTAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAG‑CAACGCCGGATCAAGGCCGATCAGGACCTCATGACCGCCGAACTTCGCCTGTAAACGACTACGAGACTGCTTAGCGATGGGCGTTAAACCATCCatgat                    >  1:202866/1‑143 (MQ=255)
                                                                                                                   tAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAG‑CAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGCACTTCGCCTGCAAACGACTACGAGCCTGCTTAGCCATGGGCGTTAAACCATCCacgacg                   >  2:319455/1‑141 (MQ=255)
                                                                                                                    aaaaaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAG‑CAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATgg                  >  1:862767/1‑143 (MQ=255)
                                                                                                                       aagcAGGCTGGCCGATAC‑‑CCCCGCCGTATGGCCCAGCAG‑CAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGtt               <  2:653097/140‑1 (MQ=255)
                                                                                                                           aGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAG‑CAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGCGCGTTAAACCATCCATGCTGGGTCTAAt           >  1:300808/1‑143 (MQ=255)
                                                                                                                               ttGCCGATACCA‑‑CCGCCGTATGTCCCAGCAG‑CAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTCAACCATCCATGATGGGTTTCATCAcc       >  1:287117/1‑143 (MQ=255)
                                                                                                                               ttGCCGATACCA‑‑CCGCCGTATGTCCCAGCAG‑CAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAcc       <  2:724853/143‑1 (MQ=255)
                                                                                                                               ttGCCGATACCA‑‑CCGCCGTATGTCCCAGCAG‑CAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAcc       >  1:866781/1‑143 (MQ=255)
                                                                                                                                      gatacccccGCCGTATGTCCCAGCAG‑CAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTgg  <  1:769287/138‑1 (MQ=255)
                                                                                                                                           ||                                                                                                                                         
TTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAG‑CAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGG  >  NC_000913/2173222‑2173498

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: ‑