Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,903,497 T→C I139T (ATC→ACC)  manY → mannose‑specific enzyme IIC component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,903,4970TC97.9% 181.0 / NA 48I139T (ATC→ACC) manYmannose‑specific enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (19/28);  minor base A (1/0);  total (20/28)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GCAGGTATCGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGC  >  NC_000913/1903361‑1903642
                                                                                                                                        |                                                                                                                                                 
gCAGGTATCGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTtcttct                                                                                                                                      <  2:501744/150‑1 (MQ=255)
  aGGTATCGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTg                                                                                                                                   <  1:232780/151‑1 (MQ=255)
    gTATCGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGtt                                                                                                                                 <  1:324475/151‑1 (MQ=255)
        cGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCt                                                                                                                              <  1:840094/150‑1 (MQ=255)
          cACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTgca                                                                                                                           <  2:88466/151‑1 (MQ=255)
            cTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTgcaa                                                                                                                          <  2:810742/150‑1 (MQ=255)
               gCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTgcaagcaa                                                                                                                      <  1:70674/151‑1 (MQ=255)
               gCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTgcaagcaa                                                                                                                      <  1:565494/151‑1 (MQ=255)
                 aaTCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATg                                                                                                                    >  2:282288/1‑151 (MQ=255)
                 aaTCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATg                                                                                                                    >  1:686621/1‑151 (MQ=255)
                        ctGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTGCCGTTTCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGtgtg                                                                                                              <  2:845693/150‑1 (MQ=255)
                          ggCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTgg                                                                                                             <  2:672471/149‑1 (MQ=255)
                             cGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTAtt                                                                                                        <  1:806441/151‑1 (MQ=255)
                                    ggCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCg                                                                                                 <  1:255022/151‑1 (MQ=255)
                                       cAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTc                                                                                               >  1:143296/1‑150 (MQ=255)
                                               gACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTgcg                                                                                      >  2:766649/1‑151 (MQ=255)
                                                   atcatcGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGctgt                                                                                  >  2:476897/1‑151 (MQ=255)
                                                     catcGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGctgtct                                                                                <  1:283997/151‑1 (MQ=255)
                                                      atcGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATAGTTGCGctgtctg                                                                               >  2:883273/1‑151 (MQ=255)
                                                        cgttcgtACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGctgtctgt                                                                              <  1:336950/150‑1 (MQ=255)
                                                        cgttcgtACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGtt                                                                             <  1:45683/151‑1 (MQ=255)
                                                          ttcgtACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTgg                                                                           >  2:1009115/1‑151 (MQ=255)
                                                           tcgtACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTgg                                                                           <  2:71197/150‑1 (MQ=255)
                                                               aCTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTAcc                                                                       >  2:250444/1‑150 (MQ=255)
                                                                cTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTAcc                                                                       >  1:237714/1‑149 (MQ=255)
                                                                  aTTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGc                                                                    <  2:697968/150‑1 (MQ=255)
                                                                  aTTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCg                                                                   <  1:73751/151‑1 (MQ=255)
                                                                       cGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGt                                                               >  2:623497/1‑150 (MQ=255)
                                                                       cGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGt                                                               >  1:407862/1‑150 (MQ=255)
                                                                        gTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTAc                                                             >  2:242169/1‑151 (MQ=255)
                                                                             tttCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAAc                                                        >  2:740269/1‑151 (MQ=255)
                                                                                ccAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATg                                                     >  1:553118/1‑151 (MQ=255)
                                                                                  aGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCt                                                   <  1:253341/151‑1 (MQ=255)
                                                                                   gCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTg                                                  >  1:822462/1‑151 (MQ=255)
                                                                                    cACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGa                                                 <  2:162048/151‑1 (MQ=255)
                                                                                     aCGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGaa                                                <  1:399475/151‑1 (MQ=255)
                                                                                      cGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAAt                                               >  2:140335/1‑151 (MQ=255)
                                                                                          gCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGc                                             <  1:1012426/149‑1 (MQ=255)
                                                                                                      gctgATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAgtg                                <  1:563597/150‑1 (MQ=255)
                                                                                                       ctgATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAgtgg                               >  2:544070/1‑150 (MQ=255)
                                                                                                         gATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAgtggtg                             >  2:337079/1‑150 (MQ=255)
                                                                                                              cGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAAt                       <  2:26087/151‑1 (MQ=255)
                                                                                                                       gACAGCGATTTCCTGGACCTACGTTTCTTCTCTGTGCCTGCAATCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTTTTGGTACCAGCGAAGTACAGATCATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAat              <  2:240828/151‑1 (MQ=255)
                                                                                                                       gACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAat              <  1:605192/151‑1 (MQ=255)
                                                                                                                       gACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAat              <  1:1009222/151‑1 (MQ=255)
                                                                                                                           gCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCg          <  1:98121/151‑1 (MQ=255)
                                                                                                                              aTTTCCAGGAAACACGTTTCTTCTCTGTTCCTGCAAGCAATGCGCGTGGCTATTCCGGCCGTCATCGTTGCGCTGACTGTTGGTACCAGCTAAGGACAGAACATGCTGAATGCGATTCCGGATGTGGTGAGAAATGGTCTGAATATc           >  2:618730/1‑147 (MQ=255)
                                                                                                                                   cTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGc  >  2:555558/1‑151 (MQ=255)
                                                                                                                                        |                                                                                                                                                 
GCAGGTATCGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGC  >  NC_000913/1903361‑1903642

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑