Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,030,732 T→C 6.8% T63A (ACC→GCC)  vsr ← DNA mismatch endonuclease of very short patch repair

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,030,7320TC 6.8% 177.1 / 3.3 59T63A (ACC→GCC) vsrDNA mismatch endonuclease of very short patch repair
Reads supporting (aligned to +/- strand):  ref base T (33/22);  new base C (2/2);  total (35/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

CCAGCCGAGTTCCTGCAAGCGACTGATATCGCGGCGATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTCCTGAACGCGAAATGCCAGGCCCTGCCCGGTTAACAGACTGGCGAGGCGCTTCTCTATCGCCGTATCACGCGTGGCAATCG  >  NC_000913/2030592‑2030872
                                                                                                                                            |                                                                                                                                            
ccAGCCGAGTTCCTGCAAGCGACTGATATCGCGGCGATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCa                                                                                                                                      <  2:184076/149‑1 (MQ=255)
 cAGCCGAGTTCCTGCAAGCGACTGATATCGCGGCGATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCAc                                                                                                                                     <  4:302141/149‑1 (MQ=255)
  aGCCGAGTTCCTGCAAGCGACTGATATCGCTGCGATCGCGCTAAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACg                                                                                                                                    <  2:89621/149‑1 (MQ=255)
   gCCGAGTTCCTGCAAGCGACTGATATCGCGGCGATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACg                                                                                                                                    >  3:36784/1‑148 (MQ=255)
     cGAGTTCCTGCAAGCGACTGATATCGCGGCGATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAg                                                                                                                                 >  6:26627/1‑149 (MQ=255)
        gTTCCTGCAAGCGACTGATATCGCGGCGATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGa                                                                                                                              >  5:191645/1‑149 (MQ=255)
            cTGCAAGCGACTGATATCGCGGCGATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATAtt                                                                                                                          >  8:56460/1‑149 (MQ=255)
               cAAGCGACTGATATCGCGGCGATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATtc                                                                                                                         <  4:36784/147‑1 (MQ=255)
                aaGCGACTGATATCGCGGCGATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATtcatc                                                                                                                      <  5:106048/149‑1 (MQ=255)
                   cGACTGATATCGCGGCGATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAAc                                                                                                                   <  3:230142/149‑1 (MQ=255)
                         atatCGCGGCGATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACTACaa                                                                                                              >  5:99752/1‑148 (MQ=255)
                         atatCGCGGCGATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGACACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACa                                                                                                               <  4:106468/147‑1 (MQ=255)
                            tCGCGGCGATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAAt                                                                                                           <  5:26627/148‑1 (MQ=255)
                                 gcgATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGAc                                                                                                     <  5:29251/149‑1 (MQ=255)
                                     tCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTcc                                                                                                 >  5:142295/1‑149 (MQ=255)
                                        cgcTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCCCGCAGCGATATTCATCACCGACAAAATCCGGACGTCCggg                                                                                              >  6:317198/1‑149 (MQ=255)
                                         gcTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCggg                                                                                              >  1:266897/1‑148 (MQ=255)
                                           tCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCggg                                                                                              >  2:290855/1‑146 (MQ=255)
                                             aaCATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGAc                                                                                         >  3:82909/1‑149 (MQ=255)
                                                  ttttACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCg                                                                                    >  8:41865/1‑149 (MQ=255)
                                                     tACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTc                                                                                  <  8:12066/148‑1 (MQ=255)
                                                       ccTATCTTCTCCAGCCAGAACTCGGTTCGATGCGCAGGCACTTTACACAGATAGCAGTGATGAGGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAACATCCGGACGTCCGGGCAGACTGGCGTCCt                                                                                <  8:124645/148‑1 (MQ=255)
                                                          aTCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTCCTGAAc                                                                            >  2:138493/1‑149 (MQ=255)
                                                          aTCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGCAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTCCTGAAc                                                                            >  2:319307/1‑149 (MQ=255)
                                                          aTCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGCAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTCCTGAAc                                                                            >  2:323885/1‑149 (MQ=255)
                                                           tcttctCCAGCCAGAACTCGGTTCGAGTCGCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTCCTGAAc                                                                            <  7:56460/148‑1 (MQ=255)
                                                             ttctCCAGCCAGAACTCGGTTCGAGTCGCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTCCTGAAcgc                                                                          >  5:33986/1‑148 (MQ=255)
                                                             ttctCCAGCCAGAACTCGGTTCGAGTCGCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTCCTGAAcgc                                                                          >  7:64234/1‑148 (MQ=255)
                                                                 ccagccagAACTCGGTTCGAGTCGCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTCCTGAACGCGAAAt                                                                     >  3:220574/1‑149 (MQ=255)
                                                                  cagccagAACTCGGTTCGAGTCGCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTCCTGAACGCGAAATg                                                                    >  2:194410/1‑149 (MQ=255)
                                                                            tCGGTTCGAGTCGCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTCCTGAACGCGAAATGCCAGGCCCTg                                                          >  3:74837/1‑149 (MQ=255)
                                                                                  cGAGTCGCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTCCTGAACGCGAAATGCCAGGCCCTGCCCgg                                                     <  1:194410/148‑1 (MQ=255)
                                                                                      tCGCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTCCTGAACGCGAAATGCCAGGCCCTGCCCGGTTAAc                                                <  7:41865/149‑1 (MQ=255)
                                                                                       cGCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTCCTGAAcg                                                                           >  8:322033/1‑121 (MQ=255)
                                                                                       cGCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTCCTGAAcg                                                                           <  7:322033/121‑1 (MQ=255)
                                                                                        gCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTCCTGAACGCGAAATGCCAGGCCCTGCCCGGTTAACa                                               <  4:220574/148‑1 (MQ=255)
                                                                                          aGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTCCTGAACGCGAAATGCCAGGCCCTGCCCGGTTAACAGAc                                            <  4:47271/149‑1 (MQ=255)
                                                                                             cACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTCCTGAACGCGAAATGCCAGGCCCTGCCCGGTTAACAGACTgg                                         <  6:33986/149‑1 (MQ=255)
                                                                                                tttAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTCCTGAACGCGAAATGCCAGGCCCTGCCCGGTTAACAGACTGGc                                        >  3:119835/1‑147 (MQ=255)
                                                                                                tttAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTCCTGAACGCGAAATGCCAGGCCCTGCCCGGTTAACAGACTGGc                                        >  2:224765/1‑147 (MQ=255)
                                                                                                tttAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTCCTGAACGCGAAATGCCAGGCCCTGCCCGGTTAACAGACTGGc                                        >  7:79996/1‑147 (MQ=255)
                                                                                                  tAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTCCTGAACGCGAAATGCCAGGCCCTGCCCGGTTAACAGACTGGCGAgg                                    >  5:16694/1‑149 (MQ=255)
                                                                                                   aaaCAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTCCTGAACGCGAAATGCCAGGCCCTGCCCGGTAAACAGACTGGCGAGgc                                   >  3:179338/1‑149 (MQ=255)
                                                                                                     aCAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGAAGTCCGGGCAGACTGGCGTCCTGAACGCGAACTGCCAGGCCCTGCCCGGTTAACAGACTGGCGAGgcgc                                 >  6:18717/1‑149 (MQ=255)
                                                                                                         aTAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTc                                                                                  >  7:9372/1‑96 (MQ=255)
                                                                                                         aTAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTc                                                                                  <  8:9372/96‑1 (MQ=255)
                                                                                                                       gatgCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTCCTGAACGCGAAATGCCAGGCCCTGCCCGGTTAACAGACTGGCGAGGCGCTTCTCTATCGCCGTATCa               >  2:216770/1‑149 (MQ=255)
                                                                                                                        atgCCAGAAGCAGCCATGGGTAAATATCCCGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTCCTGAACGCGAAATGCCAGGCCCTGCCCGGTTAACAGACTGGCGAGGCGCTTCTCTATCGCCGTATCAc              >  5:67338/1‑149 (MQ=255)
                                                                                                                        atgCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTCCTGAACGCGAAATGCCAGGCCCTGCCCGGtt                                                   >  7:25459/1‑112 (MQ=255)
                                                                                                                        atgCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTCCTGAACGCGAAATGCCAGGCCCTGCCCGGtt                                                   <  8:25459/112‑1 (MQ=255)
                                                                                                                        atgCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTCCTGAACGCGAAATGCCAGGCCCTGCCCGGTTAACAGACTGGCGAGGCGCTTCTCTATCGCCGTATCAc              >  8:293219/1‑149 (MQ=255)
                                                                                                                         tgCCAGAAGCAGCCATGGGCAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTCCTGAACGCGAAATGCCAGGCCCTGCCCGGTTAACAGACTGGCGAGGCGCTTCTCTATCGCCGTATCAcg             <  1:323885/149‑1 (MQ=255)
                                                                                                                         tgCCAGAAGCAGCCATGGGCAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTCCTGAACGCGAAATGCCAGGCCCTGCCCGGTTAACAGACTGGCGAGGCGCTTCTCTATCGCCGTATCAcg             <  1:319307/149‑1 (MQ=255)
                                                                                                                            cAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTCCTGAACGCGAAATGCCAGGCCCTGCCCGGTTAACAGACTGGCGAGGCGCTTCTCTATCGCCGTATCACGCGt          >  7:37417/1‑149 (MQ=255)
                                                                                                                             aGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTCCTGAACGCGAAATGCCAGGCCCTGCCCGGTTAACAGACTGGCGAGGCGCTTCTCTATCGCCGTATCACGCGTg         <  3:297740/149‑1 (MQ=255)
                                                                                                                             aGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTCCTGAACGCGAAATGCCAGGCCCTGCCCGGTTAACAGACTGGCGAGGCGCTTCTCTATCGCCGTATCACGCGTg         <  1:216770/149‑1 (MQ=255)
                                                                                                                               aaGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTCCTGAACGCGAAATGCCAGGCCCTGCCCGGTTAACAGACTGGCGAGGCGCTTCTCTATCGCCGTATCACGCGTgg        >  2:153357/1‑148 (MQ=255)
                                                                                                                                  cagcCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTCCTGAACGCGAAATGCCAGGCCCTGCCCGGTTAACAGACTGGCGAGGCGCTTCTCTATCGCCGTATCACGCGTGGCAAt    >  7:317260/1‑149 (MQ=255)
                                                                                                                                    gcCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTCCTGAACGCGAAATGCCAGGCCCTGCCCGGTTAACAGACTGGCGAGGCGCTTCTCTATCGCCGTATCACGCGTGGCAATcg  >  6:99890/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                            
CCAGCCGAGTTCCTGCAAGCGACTGATATCGCGGCGATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGTCCGGGCAGACTGGCGTCCTGAACGCGAAATGCCAGGCCCTGCCCGGTTAACAGACTGGCGAGGCGCTTCTCTATCGCCGTATCACGCGTGGCAATCG  >  NC_000913/2030592‑2030872

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: ‑