Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,483,242 G→C 6.8% S14T (AGC→ACC) ‡ yheT → UPF0017 family putative hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,483,2420GC 6.8% 203.2 / 4.7 59S14T (AGC→ACC) ‡yheTUPF0017 family putative hydrolase
Reads supporting (aligned to +/- strand):  ref base G (28/27);  new base C (2/2);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

AAGCCAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCTTGCGTG  >  NC_000913/3483105‑3483386
                                                                                                                                         |                                                                                                                                                
aaGCCAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGaa                                                                                                                                        <  4:1076745/148‑1 (MQ=255)
   ccAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCACTGCTGCTGAATTc                                                                                                                                     <  4:689396/148‑1 (MQ=255)
      gcgcCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCACTGCTGCTGAAt                                                                                                                                       <  4:624967/143‑1 (MQ=255)
      gcgcCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCACTGCTGCTGAAt                                                                                                                                       >  3:624967/1‑143 (MQ=255)
        gcCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGAAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATccc                                                                                                                                >  2:948757/1‑148 (MQ=255)
              tCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTcagcagcag                                                                                                                                               <  1:235838/127‑1 (MQ=255)
              tCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTcagcagcag                                                                                                                                               >  2:235838/1‑127 (MQ=255)
                      ggAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTa                                                                                                                  >  2:3666/1‑148 (MQ=255)
                                          ggcTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCa                                                                                                                              <  3:442971/116‑1 (MQ=255)
                                          ggcTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCa                                                                                                                              >  4:442971/1‑116 (MQ=255)
                                                aaGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCAtgctg                                                                                        <  1:3666/148‑1 (MQ=255)
                                                     cAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCGGTGCTGAATTCATCCCCATGCGCgg                                                                                                                       >  2:84038/1‑112 (MQ=255)
                                                     cAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCgg                                                                                                                       <  1:84038/112‑1 (MQ=255)
                                                     cAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATc                                                                                                      <  2:487609/129‑1 (MQ=255)
                                                     cAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATc                                                                                                      >  1:487609/1‑129 (MQ=255)
                                                      aGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATc                                                                                                                                  >  3:446454/1‑100 (MQ=255)
                                                      aGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATc                                                                                                                                  <  4:446454/100‑1 (MQ=255)
                                                        gAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAAt                                                                                                             >  4:402720/1‑119 (MQ=255)
                                                        gAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAAt                                                                                                             <  3:402720/119‑1 (MQ=255)
                                                           cagcTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCttt                                                                                                                   <  3:1093788/110‑1 (MQ=255)
                                                           cagcTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCttt                                                                                                                   >  4:1093788/1‑110 (MQ=255)
                                                           cagcTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGt                                                                              <  3:180396/147‑1 (MQ=255)
                                                           cagcTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACAGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGt                                                                              >  4:180396/1‑147 (MQ=255)
                                                                       aTGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCtt                                                                                                                    <  3:819010/97‑1 (MQ=255)
                                                                       aTGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCtt                                                                                                                    >  4:819010/1‑97 (MQ=255)
                                                                       aTGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCcgcg                                                                                   >  3:137336/1‑130 (MQ=255)
                                                                       aTGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCcgcg                                                                                   <  4:137336/130‑1 (MQ=255)
                                                                            gctgGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGt                                                                              >  4:500936/1‑130 (MQ=255)
                                                                            gctgGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGt                                                                              <  3:500936/130‑1 (MQ=255)
                                                                               ggAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCa                                                                                             >  4:372588/1‑112 (MQ=255)
                                                                               ggAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCa                                                                                             <  3:372588/112‑1 (MQ=255)
                                                                                gAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCaa                                                                                                              <  4:326766/94‑1 (MQ=255)
                                                                                gAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCaa                                                                                                              >  3:326766/1‑94 (MQ=255)
                                                                                gAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCa                                                                                                        >  1:106189/1‑100 (MQ=255)
                                                                                gAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCa                                                                                                        <  2:106189/100‑1 (MQ=255)
                                                                                gAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCAc                                                        <  1:948757/148‑1 (MQ=255)
                                                                                    gCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAAt                                                                                                             >  3:260174/1‑91 (MQ=255)
                                                                                    gCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAAt                                                                                                             <  4:260174/91‑1 (MQ=255)
                                                                                      cAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGaa                                                              >  3:145002/1‑136 (MQ=255)
                                                                                      cAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGaa                                                              <  4:145002/136‑1 (MQ=255)
                                                                                       aaaGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTc                                                          <  4:477324/139‑1 (MQ=255)
                                                                                       aaaGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTc                                                          >  3:477324/1‑139 (MQ=255)
                                                                                        aaGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCACTGCTGCTGAAtt                                                                                                                                      >  4:744942/1‑62 (MQ=25)
                                                                                         aGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTg                                               >  4:40054/1‑148 (MQ=255)
                                                                                             aCTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCcg                                                                                     >  1:439681/1‑106 (MQ=255)
                                                                                             aCTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCcg                                                                                     <  2:439681/106‑1 (MQ=255)
                                                                                                     cgcAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAAcc                                                                                              >  1:44299/1‑89 (MQ=255)
                                                                                                     cgcAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAAcc                                                                                              <  2:44299/89‑1 (MQ=255)
                                                                                                                            aaTGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCa                                                                                                        <  2:146602/56‑1 (MQ=255)
                                                                                                                            aaTGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCa                                                                                                        >  1:146602/1‑56 (MQ=255)
                                                                                                                             aTGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGttt                                                                            <  1:1094003/83‑1 (MQ=255)
                                                                                                                             aTGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGttt                                                                            >  2:1094003/1‑83 (MQ=255)
                                                                                                                              tGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAg                                           >  1:427195/1‑115 (MQ=255)
                                                                                                                              tGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAg                                           <  2:427195/115‑1 (MQ=255)
                                                                                                                              tGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGAt          <  3:40054/148‑1 (MQ=255)
                                                                                                                                 aTTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGt                                                   <  3:536496/104‑1 (MQ=255)
                                                                                                                                 aTTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGt                                                   >  4:536496/1‑104 (MQ=255)
                                                                                                                                   tCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCt                                      <  1:459543/115‑1 (MQ=255)
                                                                                                                                   tCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCt                                      >  2:459543/1‑115 (MQ=255)
                                                                                                                                   tCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCTTGc     >  4:911450/1‑148 (MQ=255)
                                                                                                                                   tCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCTTGc     >  3:520753/1‑148 (MQ=255)
                                                                                                                                     agcagcagTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTgg                                              <  4:556461/105‑1 (MQ=255)
                                                                                                                                     agcagcagTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTgg                                              >  3:556461/1‑105 (MQ=255)
                                                                                                                                      gcagcagTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCt                                                                                <  4:720313/70‑1 (MQ=255)
                                                                                                                                      gcagcagTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCt                                                                                >  3:720313/1‑70 (MQ=255)
                                                                                                                                      gcagcagTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGATGGCGATTTTGTCGATCTTGCGTg  >  4:724482/1‑148 (MQ=255)
                                                                                                                                         gcagTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTtcgtcg                                                                       >  3:299698/1‑76 (MQ=255)
                                                                                                                                         gcagTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTtcgtcg                                                                       <  4:299698/76‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                
AAGCCAGCGCCAAATCCGGCCTGGAAGAGTGCGAAATGGCATGGCTGGAAGCCCAGGAGCAGCTTGAGCAGATGCTGCTGGAAGGCCAAAGCAACTGATGGCGCAGATAACGACGACCGATGCCAATGAATTCAGCAGCAGTGCTGAATTCATCCCCATGCGCGGCTTTAGCAATTGTCATCTGCAAACCATGCTGCCGCGTCTGTTTCGTCGTCAGGTGAAATTCACCCCGTACTGGCAGCGGCTGGAGTTGCCCGACGGCGATTTTGTCGATCTTGCGTG  >  NC_000913/3483105‑3483386

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