Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,037,956 A→C 9.2% W51G (TGG→GGG)  recJ ← ssDNA exonuclease, 5' ‑‑> 3'‑specific

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,037,9560AC 9.2% 137.9 / 4.2 54W51G (TGG→GGG) recJssDNA exonuclease, 5' ‑‑> 3'‑specific
Reads supporting (aligned to +/- strand):  ref base A (20/29);  new base C (3/2);  total (23/31)
Fisher's exact test for biased strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

ACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCTGTTGTTT  >  NC_000913/3037810‑3038103
                                                                                                                                                  |                                                                                                                                                   
aCCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAgg                                                                                                                                                   >  4:145429/1‑149 (MQ=255)
     gcgagcgCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGgcagc                                                                                                                                              <  6:29017/149‑1 (MQ=255)
       gcgcgcATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCGTTTTCGACGCCACTCAGTTGCTGCCAGGgcagca                                                                                                                                             >  6:134323/3‑148 (MQ=255)
       gagcgCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCCTGgcagcaa                                                                                                                                            <  7:166167/149‑2 (MQ=255)
                gCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTaa                                                                                                                                    <  2:14588/148‑1 (MQ=255)
                   aGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCCTGgcagcaa                                                                                                                                            >  7:323126/1‑136 (MQ=255)
                   aGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCCTGgcagcaa                                                                                                                                            <  8:323126/137‑2 (MQ=255)
                   aGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACt                                                                                                                                >  1:353683/1‑149 (MQ=255)
                   aGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACt                                                                                                                                <  4:129024/149‑1 (MQ=255)
                        gCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAaca                                                                                                                                  <  8:175001/142‑1 (MQ=255)
                        gCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAaca                                                                                                                                  >  7:175001/1‑142 (MQ=255)
                        gCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGtt                                                                                                                           >  2:13339/1‑149 (MQ=255)
                         cTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTc                                                                                                                          <  1:23272/149‑1 (MQ=255)
                         cTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGGTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTc                                                                                                                          <  3:41953/149‑1 (MQ=255)
                          ttAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGCCCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTc                                                                                                                          >  3:297632/1‑148 (MQ=255)
                             gagCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCCTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAg                                                                                                                       <  7:113948/148‑1 (MQ=255)
                               gCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACACTAATCCGCGTTCCTTCGCGAAAAGCGGTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTc                                                                                                                    <  2:348016/149‑1 (MQ=255)
                                ccGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCt                                                                                                                   >  8:212936/1‑149 (MQ=255)
                                 cGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCtt                                                                                                                  >  8:88026/1‑149 (MQ=255)
                                        gTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCa                                                                                                             <  5:63158/147‑1 (MQ=255)
                                        gTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACt                                                                                                           <  1:13339/149‑1 (MQ=255)
                                         tCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTg                                                                                                          >  3:253010/1‑149 (MQ=255)
                                                      gTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGc                                                                                             <  6:271449/149‑1 (MQ=255)
                                                      gTCGAAATCACAGACCAAAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGc                                                                                             <  8:357724/149‑1 (MQ=255)
                                                        cGAAATCACCGACCACAATAATCCGCGGTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGTCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTg                                                                                           <  4:30998/149‑1 (MQ=255)
                                                           aaTCACCGACCACAATAATCAGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCg                                                                                        >  8:111942/1‑149 (MQ=255)
                                                                 cGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAacac                                                                                                                                 >  7:241746/1‑102 (MQ=255)
                                                                 cGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAacac                                                                                                                                 <  8:241746/102‑1 (MQ=255)
                                                                   aacacaATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGTATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACg                                                                                <  3:277730/147‑1 (MQ=255)
                                                                       caATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGgcgc                                                                            <  8:270843/149‑1 (MQ=255)
                                                                        aataatCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGgcg                                                                             >  5:315666/1‑147 (MQ=255)
                                                                           aatCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCgcagc                                                                         >  1:130956/1‑148 (MQ=255)
                                                                            atCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCaa                                                                       <  4:297632/149‑1 (MQ=255)
                                                                               cgcgTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAggg                                                                    <  2:353683/149‑1 (MQ=255)
                                                                                  gTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAgg                                                                 <  7:212936/149‑1 (MQ=255)
                                                                                       ttCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAAtt                                                            >  1:381464/1‑149 (MQ=255)
                                                                                         cgcgAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTc                                                           <  4:253010/148‑1 (MQ=255)
                                                                                                              tctcAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCa                                     >  6:82513/1‑149 (MQ=255)
                                                                                                                  aaCGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACTGCGCAGCAAGGGAGGAAATTCAGCGAGCAAGTCTGCCGTTTCa                                     <  2:53600/145‑1 (MQ=255)
                                                                                                                  aaCGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATc                                   <  6:315666/147‑1 (MQ=255)
                                                                                                                  aaCGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGa                                 <  8:60741/149‑1 (MQ=255)
                                                                                                                  aaCGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGa                                 >  5:9449/1‑149 (MQ=255)
                                                                                                                  aaCGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGa                                 >  3:121787/1‑149 (MQ=255)
                                                                                                                  aaCGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGAGCAAGTCTGCCGTTTCATc                                   <  7:111942/147‑1 (MQ=255)
                                                                                                                    cGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATc                                   <  6:38293/145‑1 (MQ=255)
                                                                                                                    cGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATc                                   >  5:38293/1‑145 (MQ=255)
                                                                                                                                      aggtatgctgcCCTGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTg              >  6:293612/6‑148 (MQ=255)
                                                                                                                                        aGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTAt           <  7:88026/149‑1 (MQ=255)
                                                                                                                                         tatgctgcCCTGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGCCTTCACGGCGACGAAGTTGTATc          >  8:142723/3‑149 (MQ=255)
                                                                                                                                            gctgcCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAgg                                                                 <  1:40306/91‑1 (MQ=255)
                                                                                                                                            gctgcCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAgg                                                                 >  2:40306/1‑91 (MQ=255)
                                                                                                                                             ctgcCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCtgt       <  6:233089/148‑1 (MQ=255)
                                                                                                                                              tgcCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGaa                  >  1:283954/1‑136 (MQ=255)
                                                                                                                                              tgcCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGaa                  <  2:283954/136‑1 (MQ=255)
                                                                                                                                                  aGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCTGTTGttt  >  4:350555/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCTGTTGTTT  >  NC_000913/3037810‑3038103

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