Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,906,503 T→A 8.0% S265T (TCC→ACC)  ygcG → TPM domain protein, putative phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,906,5030TA 8.0% 171.1 / 4.8 50S265T (TCC→ACC) ygcGTPM domain protein, putative phosphatase
Reads supporting (aligned to +/- strand):  ref base T (22/24);  new base A (2/2);  total (24/26)
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.86e-01

ATTTATTCCTGTTCTCCGATAAAATTTTTTCCATTGCTGTATTTTCCTTCTTTTTTACTTTCACCATATTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTA  >  NC_000913/2906360‑2906645
                                                                                                                                               |                                                                                                                                              
aTTTATTCCTGTTCTCCGATAAAATTTTTTCCATTGCTGTATTTTCCTTCTTTTTTACTTTCACCATATTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGAACCGCt                                                                                                                                           >  3:784172/1‑149 (MQ=255)
  ttATTCCTGTTCTCCGATAAAATTTTTTCCATTGCTGTATTTTCCTTCTTTTTTACTTTCACCATATTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTg                                                                                                                                          >  2:474013/1‑148 (MQ=255)
  ttATTCCTGTTCTCCGATAAAATTTTTTCCATTGCTGTATTTTCCTTCTTTTTTACTTTCACCATATTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTg                                                                                                                                          >  1:669697/1‑148 (MQ=255)
   tATTCCTGTTCTCCGATAAAATTTTTTCCATTGCTGTATTTTCCTTCTTTTTTACTTTCACCATATTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGa                                                                                                                                        <  3:628366/149‑1 (MQ=255)
   tATTCCTGTTCTCCGATAAAATTTTTTCCATTGCTGTATTTTCCTTCTTTTTTACTTTCACCATATTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGAACCGCTgtt                                                                                                                                        <  4:784172/149‑3 (MQ=255)
     ttCCTGTTCTCCGATAAAATTTTTTCCATTGCTGTATTTTCCTTCTTTTTTACTTTCACCATATTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAgg                                                                                                                                      >  3:419197/1‑149 (MQ=255)
        cTGTTCTCCGATAAAATTTTTTCCATTGCTGTATTTTCCTTCTTTTTTACTTTCACCATATTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCt                                                                                                                                    <  2:637491/148‑1 (MQ=255)
        cTGTTCTCCGATAAAATTTTTTCCATTGCTGTATTTTCCTTCTTTTTTACTTTCACCATATTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCt                                                                                                                                    <  2:250880/148‑1 (MQ=255)
        cTGTTCTCCGATAAAATTTTTTCCATTGCTGTATTTTCCTTCTTTTTTACTTTCACCATATTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTc                                                                                                                                   <  3:416881/149‑1 (MQ=255)
                  aTAAAATTTTTTCCATTGCTGTATTTTCCTTCTTTTTTACTTTCACCATATTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAgcgg                                                                                                                          <  2:118933/148‑1 (MQ=255)
                       aTTTTTTCCATTGCTGTATTTTCCTTCTTTTTTACTTTCACCATATTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGAt                                                                                                                     >  1:635106/1‑148 (MQ=255)
                        ttttttCCATTGCTGTATTTTCCTTCTTTTTTACTTTCACCATATTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGAtt                                                                                                                    >  1:564403/1‑148 (MQ=255)
                         tttttCCATTGCTGTATTTTCCTTCTTTTTTACTTTCACCATATTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGAttt                                                                                                                   <  2:669697/148‑1 (MQ=255)
                              ccATTGCTGTATTTTCCTTCTTTTTTACTTTCACCATATTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCt                                                                                                                                           <  4:86466/119‑1 (MQ=255)
                              ccATTGCTGTATTTTCCTTCTTTTTTACTTTCACCATATTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCt                                                                                                                                           >  3:86466/1‑119 (MQ=255)
                                aTTGCTGTATTTTCCTTCTTTTTTACTTTCACCATATTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAgcggcg                                                                                                            <  1:392389/148‑1 (MQ=255)
                                aTTGCTGTATTTTCCTTCTTTTTTACTTTCACCATATTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAgcggcg                                                                                                            <  2:635106/148‑1 (MQ=255)
                                        aTTTTCCTTCTTTTTTACTTTCACCATATTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAg                                                                                                                             <  1:44506/123‑1 (MQ=255)
                                        aTTTTCCTTCTTTTTTACTTTCACCATATTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAg                                                                                                                             >  2:44506/1‑123 (MQ=255)
                                        aTTTTCCTTCTTTTTTACTTTCACCATATTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGt                                                                                                    >  1:656972/1‑148 (MQ=255)
                                              cttcttTTTTACTTTCACCATATTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTtcttct                                                                                              <  1:481570/148‑1 (MQ=255)
                                                tcttTTTTACTTTCACCATATTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGa                                                                                           >  4:833256/1‑149 (MQ=255)
                                                          tttCACCATATTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGgcggcg                                                                                 >  3:33100/1‑149 (MQ=255)
                                                                  tatTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCgg                                                                          >  1:183820/1‑148 (MQ=255)
                                                                                aCCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCa                                                            >  1:399753/1‑148 (MQ=255)
                                                                                 ccTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCaa                                                           <  2:564403/148‑1 (MQ=255)
                                                                                  cTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGttcttc                                                                                               >  4:729036/1‑111 (MQ=255)
                                                                                  cTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGttcttc                                                                                               <  3:729036/111‑1 (MQ=255)
                                                                                  cTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAAcc                                                         >  4:548770/1‑149 (MQ=255)
                                                                                           tgtACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAAt                                                <  4:352690/149‑1 (MQ=255)
                                                                                             tACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATg                                               <  1:239611/148‑1 (MQ=255)
                                                                                               cTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCa                                            <  4:419197/149‑1 (MQ=255)
                                                                                                ttCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAg                                                                                                                 >  4:883730/1‑79 (MQ=255)
                                                                                                ttCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAg                                                                                                                 <  3:883730/79‑1 (MQ=255)
                                                                                                     aaaaGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCag                                                                              >  2:464687/1‑108 (MQ=255)
                                                                                                     aaaaGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCAt                                                                              <  1:464687/109‑1 (MQ=255)
                                                                                                      aaaGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCgg                                      <  1:474013/148‑1 (MQ=255)
                                                                                                           gCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGGTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTg                                >  4:321888/1‑149 (MQ=255)
                                                                                                               cTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCAt                            <  3:548770/149‑1 (MQ=255)
                                                                                                                tGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCAt                            <  2:656972/148‑1 (MQ=255)
                                                                                                                 ggTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCAttt                          <  3:321888/149‑1 (MQ=255)
                                                                                                                        gTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATc                   >  3:214826/1‑149 (MQ=255)
                                                                                                                               ttCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATaaa            >  3:553558/1‑149 (MQ=255)
                                                                                                                                tCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGGGGTTCTTCTGGGGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATaaa            >  2:531992/1‑148 (MQ=255)
                                                                                                                                tCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGGTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATaaa            >  2:238933/1‑148 (MQ=255)
                                                                                                                                 cAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAag           <  2:183820/148‑1 (MQ=255)
                                                                                                                                  aGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAagtc         <  4:214826/149‑1 (MQ=255)
                                                                                                                                       ccaGCGGAACCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCAtt                           <  1:100234/125‑1 (MQ=255)
                                                                                                                                       ccaGCGGAACCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCCGAGGCTGGCAtt                           >  2:100234/2‑126 (MQ=255)
                                                                                                                                          gCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTa  >  2:143063/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                              
ATTTATTCCTGTTCTCCGATAAAATTTTTTCCATTGCTGTATTTTCCTTCTTTTTTACTTTCACCATATTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTA  >  NC_000913/2906360‑2906645

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