Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 414,832 G→T 6.4% H308N (CAT→AAT)  sbcC ← exonuclease, dsDNA, ATP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913414,8320GT 6.4% 218.8 / 5.6 62H308N (CAT→AAT) sbcCexonuclease, dsDNA, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (29/29);  new base T (2/2);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.06e-01

TCCTGTAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGACGGCGGC  >  NC_000913/414693‑414975
                                                                                                                                           |                                                                                                                                               
tCCTGTAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATg                                                                                                                                                 >  4:263485/1‑140 (MQ=255)
tCCTGTAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATg                                                                                                                                                 <  3:263485/140‑1 (MQ=255)
   tGTAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGcgccg                                                                                                                                      <  3:178970/148‑1 (MQ=255)
         cAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTGAATCTGCTGGCGAATATGCGCCAGCGCCGCGCtgt                                                                                                                                >  4:741609/1‑148 (MQ=255)
              aTTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCt                                                                                                                           <  3:741609/148‑1 (MQ=255)
                         gttgCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTc                                                                                                                <  1:792642/148‑1 (MQ=255)
                                     gTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTc                                                                                                                >  4:118984/1‑136 (MQ=255)
                                     gTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTc                                                                                                                <  3:118984/136‑1 (MQ=255)
                                          tCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGAt                                                                                               >  3:475095/1‑148 (MQ=255)
                                                gACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTg                                                                                         >  4:138611/1‑148 (MQ=255)
                                                gACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTg                                                                                         <  3:138611/148‑1 (MQ=255)
                                                    gCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCgg                                                                                     >  4:965255/1‑148 (MQ=255)
                                                     cTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGt                                                                                    <  4:475095/148‑1 (MQ=255)
                                                         gccgcGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATGTACTTCTGCAATCTGCTGGCGAATATGCACCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTtgtg                                                                                <  3:965255/148‑1 (MQ=255)
                                                             cgTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTg                                                                                         >  2:354234/1‑135 (MQ=255)
                                                             cgTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTg                                                                                         <  1:354234/135‑1 (MQ=255)
                                                                  tggCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAg                                                                                                            >  3:593929/1‑111 (MQ=255)
                                                                  tggCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAg                                                                                                            <  4:593929/111‑1 (MQ=255)
                                                                       aaTGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGCATATTCGCCAgc                                                                                                                                          >  3:829113/1‑76 (MQ=255)
                                                                       aaTGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGCATATTCGCCAgc                                                                                                                                          <  4:829113/76‑1 (MQ=255)
                                                                                         cATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGt                                                                                                                  <  1:1650/82‑1 (MQ=255)
                                                                                         cATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGt                                                                                                                  >  2:1650/1‑82 (MQ=255)
                                                                                               gCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTa                                                                        <  3:106923/118‑1 (MQ=255)
                                                                                               gCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTa                                                                        >  4:106923/1‑118 (MQ=255)
                                                                                               gCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTg                                                    >  4:830854/1‑138 (MQ=255)
                                                                                               gCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTg                                                    <  3:830854/138‑1 (MQ=255)
                                                                                                 tctGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTtg                                                                                  <  4:575959/106‑1 (MQ=255)
                                                                                                 tctGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTtg                                                                                  >  3:575959/1‑106 (MQ=255)
                                                                                                    gTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCgcgc                                                                                                                                   <  2:627163/54‑1 (MQ=255)
                                                                                                    gTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCgcgc                                                                                                                                   >  1:627163/1‑54 (MQ=255)
                                                                                                     tAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAgtg                                                                                                          <  4:963917/78‑1 (MQ=255)
                                                                                                     tAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAgtg                                                                                                          >  3:963917/1‑78 (MQ=255)
                                                                                                         cGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGa                                                                                                                       >  2:314219/1‑61 (MQ=255)
                                                                                                         cGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGa                                                                                                                       <  1:314219/61‑1 (MQ=255)
                                                                                                         cGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCa                                                              <  4:192813/118‑1 (MQ=255)
                                                                                                         cGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCa                                                              >  3:192813/1‑118 (MQ=255)
                                                                                                              aTTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCTGACTAAGCGCCGCCAGTTGAGGTTGCGCttttt                                            >  2:282849/1‑131 (MQ=255)
                                                                                                              aTTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCttttt                                            <  1:282849/131‑1 (MQ=255)
                                                                                                                   cttcttCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTg                                                                                         >  2:790403/1‑81 (MQ=255)
                                                                                                                   cttcttCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTg                                                                                         <  1:790403/81‑1 (MQ=255)
                                                                                                                    ttcttcAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTg                                                    <  1:130400/117‑1 (MQ=255)
                                                                                                                    ttcttcAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTg                                                    >  2:130400/1‑117 (MQ=255)
                                                                                                                     tcttcAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCCAACGcc                           >  2:247607/1‑141 (MQ=255)
                                                                                                                     tcttcAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCCAACGcc                           <  1:247607/141‑1 (MQ=255)
                                                                                                                        tcAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAgg                 >  3:497061/1‑148 (MQ=255)
                                                                                                                           aTCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCa                                                              >  3:584459/1‑100 (MQ=255)
                                                                                                                           aTCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCa                                                              <  4:584459/100‑1 (MQ=255)
                                                                                                                               gctgGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCttt                                              <  1:501918/112‑1 (MQ=255)
                                                                                                                               gctgGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCttt                                              >  2:501918/1‑112 (MQ=255)
                                                                                                                               gctgGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCctgctg          >  1:624460/1‑148 (MQ=255)
                                                                                                                                ctgGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGa         <  4:497061/148‑1 (MQ=255)
                                                                                                                                 tgGCGCATATTCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGAc        >  3:765141/1‑148 (MQ=255)
                                                                                                                                  ggCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACg                             >  1:488202/1‑126 (MQ=255)
                                                                                                                                  ggCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACg                             <  2:488202/126‑1 (MQ=255)
                                                                                                                                  ggCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAgg                 <  2:458446/138‑1 (MQ=255)
                                                                                                                                  ggCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAgg                 >  1:458446/1‑138 (MQ=255)
                                                                                                                                   gCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTgcg                                                  >  2:300551/1‑104 (MQ=255)
                                                                                                                                   gCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTgcg                                                  <  1:300551/104‑1 (MQ=255)
                                                                                                                                      aaTATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGcc                           <  4:802361/124‑1 (MQ=255)
                                                                                                                                      aaTATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGcc                           >  3:802361/1‑124 (MQ=255)
                                                                                                                                       atatTCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGACggcggc  <  4:765141/148‑1 (MQ=255)
                                                                                                                                       atatGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGACggcggc  >  2:327795/1‑148 (MQ=255)
                                                                                                                                         atGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCgg                                   <  1:133904/113‑1 (MQ=255)
                                                                                                                                         atGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCgg                                   >  2:133904/1‑113 (MQ=255)
                                                                                                                                           |                                                                                                                                               
TCCTGTAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGACGGCGGC  >  NC_000913/414693‑414975

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