Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,302,151 G→A 6.9% G223S (GGC→AGC)  yhbU → U32 peptidase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,302,1510GA 6.9% 187.0 / 2.9 59G223S (GGC→AGC) yhbUU32 peptidase family protein
Reads supporting (aligned to +/- strand):  ref base G (24/30);  new base A (2/2);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.75e-01

TCTGTCGTCGTATCTGACGGGTGAGTCGCCCAACACCATAGGCGCGTGTTCTCCGGCCCGTTTCGTGCGCTGGCAACAAACGCCGCAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTCTCAATACCCTGGAACTGCTGCCGGAGTTAATGGCGGCGAAT  >  NC_000913/3302015‑3302279
                                                                                                                                        |                                                                                                                                
tCTGTCGTCGTATCTGACGGGTGAGTCGCCCAACACCATAGGCGCGTGTTCTCCGGCCCGTTTCGTGCGCTGGCAACAAACGCCGCAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACg                                                                                                                         >  8:188225/1‑146 (MQ=255)
      gtcgtATCTGACGGGTGAGTCGCCCAACACCATAGGCGCGTGTTCTCCGGCCCGTTTCGTGCGCTGGCAACAAACGCCGCAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTAt                                                                                                                 <  1:363053/148‑1 (MQ=255)
               gACGGGTGAGTCGCCCAACACCATAGGCGCGTGTTCTCCGGCCCGTTTCGTGCGCTGGCAACAAACGCCGCAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTa                                                                                                        >  5:384702/1‑148 (MQ=255)
                aCGGGTGAGTCGCCCAACACCATAGGCGCGTGTTCTCCGGCCCGTTTCGTGCGCTGGCAACAAACGCCGCAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTa                                                                                                        <  2:101356/147‑1 (MQ=255)
                    gTGAGTCGCCCAACACCATAGGCGCGTGTTCTCCGGCCCGTTTCGTGCGCTGGCAACAAACGCCGCAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTaa                                                                                                   <  7:178457/148‑1 (MQ=255)
                       aGTCGCCCAACACCATAGGCGCGTGTTCTCCGGCCCGTTTCGTGCGCTGGCAACAAACGCCGCAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAgg                                                                                                >  5:96756/1‑148 (MQ=255)
                           gCCCAACACCATAGGCGCGTGTTCTCCGGCCCGTTTCGTGCGCTGGCAACAAACGCCGCAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGt                                                                                            <  8:125518/148‑1 (MQ=255)
                                 cacCATAGGCGCGTGTTCTCCGGCCCGTTTCGTGCGCTGGCAACAAACGCCGCAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTaa                                                                                                   >  7:75986/1‑135 (MQ=255)
                                 cacCATAGGCGCGTGTTCTCCGGCCCGTTTCGTGCGCTGGCAACAAACGCCGCAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTaa                                                                                                   <  8:75986/135‑1 (MQ=255)
                                    cATAGGCGCGTGTTCTCCGGCCCGTTTCGTGCGCTGGCAACAAACGCCGCAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCtggtg                                                                                   >  7:180587/1‑148 (MQ=255)
                                       aGGCGCGTGTTCTCCGGCCCGTTTCGTGCGCTGGCAACAAACGCCGCAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGAc                                                                                >  4:287701/1‑148 (MQ=255)
                                        ggCGCGTGTTCTCCGGCCCGTTTCGTGCGCTGGCAACAAACGCCGCAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACg                                                                               <  6:384702/148‑1 (MQ=255)
                                             gtgtTCTCCGGCCCGTTTCGTGCGCTGGCAACAAACGCCGCAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTaa                                                                                                   >  8:86980/1‑123 (MQ=255)
                                             gtgtTCTCCGGCCCGTTTCGTGCGCTGGCAACAAACGCCGCAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTaa                                                                                                   <  7:86980/123‑1 (MQ=255)
                                             gtgtTCTCCGGCCCGTTTCGTGCGCTGGCAACAAACGCCGCAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAg                                                                          >  7:298282/1‑148 (MQ=255)
                                                  ctcCGGCCCGTTTCGTGCGCTGGCAACAAACGCCGCAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCt                                                                      <  6:96756/147‑1 (MQ=255)
                                                    ccGGCCCGTTTCGTGCGCTGGCAACAAACGCCGCAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATc                                                                   >  7:166111/1‑148 (MQ=255)
                                                     cGGCCCGTTTCGTGCGCTGGCAACAAACGCCGCAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATc                                                                   <  3:287701/147‑1 (MQ=255)
                                                      gcccccTTTCGTGCGCTGGCAACAAACGCCTCAGGCGCTGGAATCCCGCCTGAACGAAGTGCTGAACGACCGTTATCACGACCGCGAAAACGCAAGTTAAACGACGCAATGAAAAGAGCTTTATCTGGAGCACGGCGAGCGATATCaa                                                                 >  2:180511/3‑147 (MQ=255)
                                                             ttCGTGCGCTGGCAACAAACGCCGCAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGCCCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCTTTATCTGGTGGACGGCGAGCGCTATCACGCGCTg                                                           >  1:117797/1‑147 (MQ=255)
                                                             ttCGTGCGCTGGCAACAAACGCCGCAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTaa                                                                                                   >  8:252296/1‑107 (MQ=255)
                                                             ttCGTGCGCTGGCAACAAACGCCGCAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTaa                                                                                                   <  7:252296/107‑1 (MQ=255)
                                                               cGTGCGCTGGCAACAAACGCCGCAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTgg                                                          >  7:266226/1‑146 (MQ=255)
                                                                   cgcTGGCAACAAACGCCGCAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGaac                                                    >  7:363190/1‑148 (MQ=255)
                                                                           acaaACGCCGCAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGAc                                                                                                            <  4:73814/84‑1 (MQ=255)
                                                                           acaaACGCCGCAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGAc                                                                                                            >  3:73814/1‑84 (MQ=255)
                                                                            caaACGCCGCAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGtc                                           <  8:76393/148‑1 (MQ=255)
                                                                             aaaCGCCGCAGGTGCTGGAATCCCGCCTTAACCAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCTACGCTAGGTAAAGGGCGTTATCTGGTGGACGCCGAGCGCTATCACGCGCTGGAAGAACCAACCAGtct                                          <  2:117797/148‑1 (MQ=255)
                                                                              aaCGCCGCAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATc                                                                   <  2:59583/122‑1 (MQ=255)
                                                                              aaCGCCGCAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATc                                                                   >  1:59583/1‑122 (MQ=255)
                                                                              aaCGCCGCAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGaaccaacc                                               <  5:39952/142‑1 (MQ=255)
                                                                              aaCGCCGCAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGaaccaacc                                               >  6:39952/1‑142 (MQ=255)
                                                                              aaCGCCGCAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACAGCGAAAACGCAGGTTATCCGACGCTATGTa                                                                                                    >  4:203477/1‑89 (MQ=255)
                                                                              aaCGCCGCAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACAGCGAAAACGCAGGTTATCCGACGCTATGTa                                                                                                    <  3:203477/89‑1 (MQ=255)
                                                                              aaCGCCGCAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACAGCGAAAACGCAGGTTATCCGACGCTATGTa                                                                                                    <  3:208237/89‑1 (MQ=255)
                                                                              aaCGCCGCAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACAGCGAAAACGCAGGTTATCCGACGCTATGTa                                                                                                    >  4:208237/1‑89 (MQ=255)
                                                                               aCGCCGCAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTCTCa                                        <  6:96541/148‑1 (MQ=255)
                                                                                 gccgcAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTCTCAAt                                      <  8:266226/148‑1 (MQ=255)
                                                                                     cAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTCTCAATAccc                                  <  8:180587/148‑1 (MQ=255)
                                                                                     cAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTCTCAATAccc                                  >  1:66953/1‑148 (MQ=255)
                                                                                        gggCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGtctc                                         <  6:170142/138‑1 (MQ=255)
                                                                                        gggCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGtctc                                         >  5:170142/1‑138 (MQ=255)
                                                                                         ggCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTCTCAATACCCTGGa                              <  8:363190/148‑1 (MQ=255)
                                                                                         ggCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTCTCAATACCCTGGa                              <  7:188225/148‑1 (MQ=255)
                                                                                                 tCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTCTCAATACCCTGGAACtgctgc                      <  2:66953/148‑1 (MQ=255)
                                                                                                  cccGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGaacca                                                  <  6:383364/119‑1 (MQ=255)
                                                                                                  cccGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGaacca                                                  >  5:383364/1‑119 (MQ=255)
                                                                                                       cTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTCTCAATACCCTGGAACTGCTGCCGGAGt                <  8:298282/148‑1 (MQ=255)
                                                                                                         gAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTCTCAATACCCTGGAACTGCTGCCGGAGtt               >  3:237944/1‑147 (MQ=255)
                                                                                                          aaCGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTCTCAATACCCTGGAACTGCTGCCGGAGTTa              <  8:166111/147‑1 (MQ=255)
                                                                                                          aaCGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTCTCAATACCCTGGAACTGCTGCCGGAGTTa              <  4:237944/147‑1 (MQ=255)
                                                                                                              aaGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAg                                             >  2:313180/1‑112 (MQ=255)
                                                                                                              aaGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAg                                             <  1:313180/112‑1 (MQ=255)
                                                                                                                  gCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGtctc                                         >  8:179368/1‑112 (MQ=255)
                                                                                                                  gCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGtctc                                         <  7:179368/112‑1 (MQ=255)
                                                                                                                   cTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTCTCAATACCCTGGAACTGCTGCCGGAGTTAATGgcggcg     >  4:45150/1‑147 (MQ=255)
                                                                                                                      aTCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTCTCAATACCCTGGAACTGCTGCCGGAGTTAATGGCGGCGAat  <  3:45150/147‑1 (MQ=255)
                                                                                                                           ccGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATc                                                                   >  7:222167/1‑77 (MQ=255)
                                                                                                                           ccGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATc                                                                   <  8:222167/77‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                
TCTGTCGTCGTATCTGACGGGTGAGTCGCCCAACACCATAGGCGCGTGTTCTCCGGCCCGTTTCGTGCGCTGGCAACAAACGCCGCAGGGGCTGGAATCCCGCCTGAACGAAGTGCTGATCGACCGTTATCAGGACGGCGAAAACGCAGGTTATCCGACGCTATGTAAAGGGCGTTATCTGGTGGACGGCGAGCGCTATCACGCGCTGGAAGAACCAACCAGTCTCAATACCCTGGAACTGCTGCCGGAGTTAATGGCGGCGAAT  >  NC_000913/3302015‑3302279

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