Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,906,502 T→A 7.9% G264G (GGT→GGA)  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,5020TA 7.9% 172.5 / 2.3 50G264G (GGT→GGA) 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 = 1.00e+00

TTTAGATTTATTCCTGTTCTCCGATAAAATTTTTTCCATTGCTGTATTTTCCTTCTTTTTTACTTTCACCATATTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTA  >  NC_000913/2906355‑2906645
                                                                                                                                                   |                                                                                                                                               
tttAGATTTATTCCTGTTCTCCGATAAAATTTTTTCCATTGCTGTATTTTCCTTCTTTTTTACTTTCACCATATTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGt                                                                                                                                                 >  1:250880/1‑148 (MQ=255)
     aTTTATTCCTGTTCTCCGATAAAATTTTTTCCATTGCTGTATTTTCCTTCTTTTTTACTTTCACCATATTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGAACCGCt                                                                                                                                           >  3:784172/1‑149 (MQ=255)
       ttATTCCTGTTCTCCGATAAAATTTTTTCCATTGCTGTATTTTCCTTCTTTTTTACTTTCACCATATTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTg                                                                                                                                          >  1:669697/1‑148 (MQ=255)
       ttATTCCTGTTCTCCGATAAAATTTTTTCCATTGCTGTATTTTCCTTCTTTTTTACTTTCACCATATTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTg                                                                                                                                          >  2:474013/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:250880/148‑1 (MQ=255)
             cTGTTCTCCGATAAAATTTTTTCCATTGCTGTATTTTCCTTCTTTTTTACTTTCACCATATTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCt                                                                                                                                    <  2:637491/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                                                                                                                                           >  3:86466/1‑119 (MQ=255)
                                   ccATTGCTGTATTTTCCTTCTTTTTTACTTTCACCATATTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCt                                                                                                                                           <  4:86466/119‑1 (MQ=255)
                                     aTTGCTGTATTTTCCTTCTTTTTTACTTTCACCATATTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAgcggcg                                                                                                            <  2:635106/148‑1 (MQ=255)
                                     aTTGCTGTATTTTCCTTCTTTTTTACTTTCACCATATTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAgcggcg                                                                                                            <  1:392389/148‑1 (MQ=255)
                                             aTTTTCCTTCTTTTTTACTTTCACCATATTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAg                                                                                                                             >  2:44506/1‑123 (MQ=255)
                                             aTTTTCCTTCTTTTTTACTTTCACCATATTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAg                                                                                                                             <  1:44506/123‑1 (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)
                                                                                                                                                   |                                                                                                                                               
TTTAGATTTATTCCTGTTCTCCGATAAAATTTTTTCCATTGCTGTATTTTCCTTCTTTTTTACTTTCACCATATTTATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTA  >  NC_000913/2906355‑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: ‑