Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,821,182 T→A V11V (GTA→GTT chbC ← N,N'‑diacetylchitobiose‑specific enzyme IIC component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,821,1820TA92.3% 162.7 / 2.4 52V11V (GTA→GTTchbCN,N'‑diacetylchitobiose‑specific enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base T (3/1);  new base A (23/25);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 6.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

GGCCCCCGCAAGGGTTAACGGCATTAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACG  >  NC_000913/1821086‑1821281
                                                                                                |                                                                                                   
ggCCCCCGCAAGGGTTAACGGCATTAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGAAcc                                                                                                  >  2:873930/1‑100 (MQ=255)
ggCCCCCGCAAGGGTTAACGGCATTAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGAAcc                                                                                                  >  1:3597/1‑100 (MQ=255)
    cccGCAAGGGTTAACGGCATTAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGAACCtttt                                                                                              <  1:1359947/100‑1 (MQ=255)
      cGCAAGGGTTAACGGCATTAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGAACCTTTTCa                                                                                            >  2:5785/1‑100 (MQ=255)
            ggTTAACGGCATTAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGAACCTTTTCAAGCGAt                                                                                      <  1:379528/100‑1 (MQ=255)
            ggTTAACGGCATTAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGAACCTTTTCAAGCGAt                                                                                      <  2:1313926/100‑1 (MQ=255)
             gTTAACGGCATTAAGCGAATAAAGCCATTTTTTTTTTTTTTTACGTGTGGCTGCTTTCCTATTTTAACTGCTTAAGGGAGGAGAACCTTTTCAAGc                                                                                         >  2:202468/1‑96 (MQ=255)
             gTTAACGGCATTAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGAACCTTTTCAAGCGATg                                                                                     <  1:1266091/100‑1 (MQ=255)
             gTTAACGGCATTAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAAGAGAACCTTTTCAAGCGATg                                                                                     <  2:966015/100‑1 (MQ=255)
              ttAACGGCATTAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGAACCTTTTCAAGCGATGc                                                                                    <  2:1173344/100‑1 (MQ=255)
                 aCGGCATTAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGAACCTTTTCAAGCGATGCAAt                                                                                 >  2:209754/1‑100 (MQ=255)
                    gCATTAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGAACCTTTTCAAGCGATGCAATAAc                                                                              <  2:1488422/100‑1 (MQ=255)
                          aGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGAACCTTTTCAAGCGATGCAATAACATTACt                                                                        >  1:837444/1‑100 (MQ=255)
                          aGCGAATAAAGCCATTTTTGATTGCACTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGAACCTTTTCAAGCGATGCAATAACATTACt                                                                        >  1:960848/1‑100 (MQ=255)
                            cGAATAAAGCCATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGAACCTTTTCAAGCGATGCAATAACATTACTCa                                                                      <  2:1401805/100‑1 (MQ=255)
                            cGAATAAAGCCATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGAACCTTTTCAAGCGATGCAATAACATTACTCa                                                                      <  2:1400564/100‑1 (MQ=255)
                                  aaGCCATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGAACCTTTTCAAGCGATGCAATAACATTACTCATAGaaa                                                                >  2:568292/1‑100 (MQ=255)
                                      cATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGAACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATa                                                            >  1:1024866/1‑100 (MQ=255)
                                      cATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGAACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATa                                                            >  2:1098403/1‑100 (MQ=255)
                                       aTTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGAACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATAc                                                           <  1:608544/100‑1 (MQ=255)
                                         ttttGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGAACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATAccc                                                         >  1:9391/1‑100 (MQ=255)
                                          tttGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGAACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCt                                                        <  1:920069/100‑1 (MQ=255)
                                                 tctttAATGTGTGGCGGCTTTCCTATTTTAACTGCAAAAGGGAGGAGAACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAAc                                                 >  2:314386/4‑100 (MQ=255)
                                                   aTTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGAACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCg                                               >  2:672445/1‑100 (MQ=255)
                                                    ttAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGAACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGc                                              >  1:674829/1‑100 (MQ=255)
                                                     tAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACACTTAAAAACCGCa                                             >  2:570600/1‑100 (MQ=255)
                                                     tAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGAACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCa                                             >  2:41065/1‑100 (MQ=255)
                                                      aaCGTGTGGCTGCTTTCCTATTTTAACTTCAAAAGGGAGGAGAACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCa                                             <  2:393339/99‑1 (MQ=255)
                                                        cGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGAACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAAtt                                          <  1:587043/100‑1 (MQ=255)
                                                         gtgtGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGAACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAAttt                                         >  1:35263/1‑100 (MQ=255)
                                                              gctgctTTCCTATTTTAACTGCAAAAGGGAGGAGAACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAat                                     <  2:1290202/99‑1 (MQ=255)
                                                               ctgctTTCCTATTTTAACTGCAAAAGGGAGGAGAACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAatat                                   >  1:820383/1‑100 (MQ=255)
                                                                 gctTTCCTATTTTAACTGCAAAAGGGAGGAGAACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTg                                 <  1:914614/100‑1 (MQ=255)
                                                                  ctTTCATATTTTAACTGCAAAAGGCAGGAGTAACTTTTCAAGCGATGCAATAACATTACACATAGAAAAATACCCTTAAAAACCGCAATTTAAATAGTGc                                <  2:234055/100‑1 (MQ=255)
                                                                    ttCCTATTTTAACTGCAAAAGGGAGGAGAACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCgg                              <  2:455258/100‑1 (MQ=255)
                                                                      ccTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTa                            >  2:1456115/1‑100 (MQ=255)
                                                                      ccTATTTTAACTGCAAAAGGGAGGAGAACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTa                            >  2:120961/1‑100 (MQ=255)
                                                                         aTTTTAACTGCAAAAGGGAGGAGAACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTg                         <  1:209754/100‑1 (MQ=255)
                                                                               aCTGCAAAAGGGAGGAGAACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTAt                   <  1:744373/100‑1 (MQ=255)
                                                                               aCTGCAAAAGGGAGGAGAACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTAt                   >  1:742262/1‑100 (MQ=255)
                                                                                cTGCAAAAGGGAGGAGAACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATg                  <  1:1337310/100‑1 (MQ=255)
                                                                                 tGCAAAAGGGAGGAGAACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGa                 <  1:1167371/100‑1 (MQ=255)
                                                                                  gCAAGAGGGAGGAGAACCTGTTCAAGCGATGCAATAGCATTACGCATAGAGAAATACCCTTAAAAACCGCAATTTAAAGATTGCGGTATTGAttt                     <  2:63364/95‑1 (MQ=255)
                                                                                  gCAAAAGGGAGGAGAACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGaa                >  2:696286/1‑100 (MQ=255)
                                                                                   cAAAAGGGAGGAGAACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGaaa               <  1:630403/100‑1 (MQ=255)
                                                                                    agaaGGGAGGAGAACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAAt              <  2:960843/98‑1 (MQ=255)
                                                                                     aaaGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATa             >  2:463658/1‑100 (MQ=255)
                                                                                     aaaGGGAGGAGAACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATa             >  2:1303138/1‑100 (MQ=255)
                                                                                     aaaGGGAGGAGAACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATa             >  2:893738/1‑100 (MQ=255)
                                                                                        gggAGGAGAACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAAct          <  1:364260/100‑1 (MQ=255)
                                                                                            ggagAACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCttt      <  2:35263/100‑1 (MQ=255)
                                                                                            ggagAACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCttt      >  2:1130097/1‑100 (MQ=255)
                                                                                                tACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACg  <  1:280108/100‑1 (MQ=255)
                                                                                                |                                                                                                   
GGCCCCCGCAAGGGTTAACGGCATTAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACG  >  NC_000913/1821086‑1821281

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