Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,971,442 C→A 10.1% A61A (GCC→GCA)  wecC → UDP‑N‑acetyl‑D‑mannosaminuronic acid dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,971,4420CA10.1% 147.1 / 3.6 46A61A (GCC→GCA) wecCUDP‑N‑acetyl‑D‑mannosaminuronic acid dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (18/23);  new base A (3/2);  total (21/25)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCTGCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGC  >  NC_000913/3971303‑3971576
                                                                                                                                           |                                                                                                                                      
ggCTGCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGTTTTGGCGAGTGTAGTAAAAACTGcc                                                                                                                                        >  1:1079024/1‑140 (MQ=255)
ggCTGCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGcc                                                                                                                                        >  1:888301/1‑140 (MQ=255)
ggCTGCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGcc                                                                                                                                        >  2:995690/1‑140 (MQ=255)
ggCTGCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGcc                                                                                                                                        >  2:905995/1‑140 (MQ=255)
ggCTGCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGcc                                                                                                                                        >  2:1536019/1‑140 (MQ=255)
ggCTGCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGcc                                                                                                                                        >  1:557375/1‑140 (MQ=255)
ggCTGCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGcc                                                                                                                                        >  2:250988/1‑140 (MQ=255)
ggCTGCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGcc                                                                                                                                        >  2:251023/1‑140 (MQ=255)
ggCTGCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGcc                                                                                                                                        >  1:257786/1‑140 (MQ=255)
    gCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCAGTAg                                                                                                                                    >  1:848857/1‑140 (MQ=255)
    gCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCAGTAg                                                                                                                                    >  1:174583/1‑140 (MQ=255)
     ccAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGa                                                                                                                                   >  2:1487905/1‑140 (MQ=255)
                    tttGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTAcg                                                                                                                    >  1:686731/1‑140 (MQ=255)
                                 aaaaaCAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCacgacg                                                                                                       <  1:167250/140‑1 (MQ=255)
                                 aaaaaCAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGACGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCacgacg                                                                                                       <  1:264112/140‑1 (MQ=255)
                                  aaaaCAGGTAATTGGTGTCGATATCAGCCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGc                                                                                                      <  1:620750/140‑1 (MQ=255)
                                                   tCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGcc                                                                                     >  1:178241/1‑140 (MQ=255)
                                                      atatCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTgg                                                                                  >  1:69705/1‑140 (MQ=255)
                                                        atCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAATAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCt                                                                                <  1:1582353/140‑1 (MQ=255)
                                                           aaccaacATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTATAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGAt                                                                             >  2:1575784/1‑140 (MQ=255)
                                                                 cATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGt                                                                       >  1:211229/1‑140 (MQ=255)
                                                                       gTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCAc                                                                 <  2:786238/140‑1 (MQ=255)
                                                                       gTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCAc                                                                 >  2:690310/1‑140 (MQ=255)
                                                                       gTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCAc                                                                 <  2:1258631/140‑1 (MQ=255)
                                                                       gTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCAc                                                                 >  1:1101125/1‑140 (MQ=255)
                                                                       gTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGAAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCAc                                                                 >  2:690203/1‑140 (MQ=255)
                                                                        ttGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGTGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACg                                                                <  2:887974/140‑1 (MQ=255)
                                                                        ttGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACg                                                                <  2:1592451/140‑1 (MQ=255)
                                                                        ttGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACg                                                                <  1:1085148/140‑1 (MQ=255)
                                                                        ttGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACg                                                                <  1:660771/140‑1 (MQ=255)
                                                                        ttGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACg                                                                >  1:606867/1‑140 (MQ=255)
                                                                        ttGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCAGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACg                                                                <  2:848857/140‑1 (MQ=255)
                                                                        ttGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCAGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACg                                                                <  2:174583/140‑1 (MQ=255)
                                                                                  cAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAgg                                                      >  2:1414125/1‑140 (MQ=255)
                                                                                  cAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAgg                                                      >  1:290154/1‑140 (MQ=255)
                                                                                    aTCGTGGCGTAATCCTTTTCGTCGATCGTGATTTGGCGAGTGTTGTAAATACTGCCGTAGAAGGCGGTTTTTTTCGAGCGAGCTCGTCGCCAGTTGAAGCGGATGCCTGGCTGTTTGCTGTACCCACGCCGTTTAAGGGc                                                    <  1:1428258/140‑1 (MQ=255)
                                                                                         ggCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCa                                               >  1:1625545/1‑140 (MQ=255)
                                                                                         ggCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCa                                               >  1:751097/1‑140 (MQ=255)
                                                                                              aaTCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGc                                          <  2:501834/140‑1 (MQ=255)
                                                                                              aaTCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGc                                          <  2:501164/140‑1 (MQ=255)
                                                                                               aTCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGcc                                         >  2:1456266/1‑140 (MQ=255)
                                                                                                 ccATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGG‑‑TTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGat                                     <  1:1221804/140‑1 (MQ=255)
                                                                                                            aaCCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTAc                            <  1:33600/140‑1 (MQ=255)
                                                                                                            aaCCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTAc                            <  1:263029/140‑1 (MQ=255)
                                                                                                                 gATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGa                       <  1:546153/140‑1 (MQ=255)
                                                                                                                 gATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGa                       <  1:546753/140‑1 (MQ=255)
                                                                                                                  aTTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTgta                      >  1:156649/1‑138 (MQ=255)
                                                                                                                        cGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGgct                <  2:1065216/140‑1 (MQ=255)
                                                                                                                        cGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGgct                <  2:1065151/140‑1 (MQ=255)
                                                                                                                        cGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGgct                <  2:1064652/140‑1 (MQ=255)
                                                                                                                           gtgtAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGgctgct             <  1:234108/140‑1 (MQ=255)
                                                                                                                            tgtAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTgctc            >  1:872836/1‑140 (MQ=255)
                                                                                                                              tagtaAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTgctcgc          <  2:47032/140‑1 (MQ=255)
                                                                                                                                      aCTGCCGTAGAAGGCGGTTTTTTACGTGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGTCCTACGTTGAATCGGCTGCGCGCTCCATTgc  >  1:779482/1‑140 (MQ=255)
                                                                                                                                      aCTGCAGTAGAAGGCTTTTTTTTACGATCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGATCACACGCCGTTTAAGTGCTAACAAGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTgc  >  1:779481/1‑140 (MQ=255)
                                                                                                                                           |                                                                                                                                      
GGCTGCCAACGGCAGCAGCGTTTGCCTCACGGCAAAAACAGGTAATTGGTGTCGATATCAACCAACATGCGGTTGATACCATCAATCGTGGCGAAATCCATATCGTCGAACCTGATTTGGCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCGCTCCATTGC  >  NC_000913/3971303‑3971576

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