Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,580,834 T→C 8.0% T210T (ACA→ACG)  hsdS ← specificity determinant for hsdM and hsdR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,580,8340TC 8.0% 132.7 / 2.5 50T210T (ACA→ACG) hsdSspecificity determinant for hsdM and hsdR
Reads supporting (aligned to +/- strand):  ref base T (24/22);  new base C (2/2);  total (26/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 = 7.54e-01

CACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGG  >  NC_000913/4580686‑4580976
                                                                                                                                                    |                                                                                                                                              
cACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCt                                                                                                                                                >  2:697726/1‑149 (MQ=255)
cACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCt                                                                                                                                                >  2:792877/1‑149 (MQ=255)
 aCGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTg                                                                                                                                               >  2:541246/1‑149 (MQ=255)
   gTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTc                                                                                                                                             >  1:404305/1‑149 (MQ=255)
                 cGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAAc                                                                                                                               <  2:706928/149‑1 (MQ=255)
                 cGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAAc                                                                                                                               <  1:792877/149‑1 (MQ=255)
                 cGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAAc                                                                                                                               <  1:697726/149‑1 (MQ=255)
                          ggATGCCCAACACCACTTTCTTTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGccccc                                                                                                                      <  2:123900/149‑1 (MQ=255)
                           gATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGcccccc                                                                                                                     >  1:114804/1‑149 (MQ=255)
                           gATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTACATTAACTGCGcccccc                                                                                                                     >  1:634/1‑149 (MQ=255)
                             tGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCaa                                                                                                                   <  1:541246/149‑1 (MQ=255)
                             tGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCaa                                                                                                                   >  1:725725/1‑149 (MQ=255)
                             tGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACAGCGCCCCCCaa                                                                                                                   >  2:127784/1‑149 (MQ=255)
                                 cAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATAcc                                                                                                               >  2:471322/1‑149 (MQ=255)
                                    caccacTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGAGAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCt                                                                                                            <  2:130474/149‑1 (MQ=255)
                                         cTTTCATTTGGCTTTGATGAAAGACCATTACGTAATCCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACg                                                                                                       <  1:189224/149‑1 (MQ=255)
                                          tttCATTTGGCTTTGATGAAAGAGCATTACGTAATTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTATGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGa                                                                                                      >  1:429234/1‑149 (MQ=255)
                                          tttCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCCGTCAATTTTCCATTAACTGCGCCCCCCACTACCGCTTGACgg                                                                                                      >  2:556520/1‑148 (MQ=255)
                                          tttCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCCGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGa                                                                                                      >  2:90473/1‑149 (MQ=255)
                                                ttGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACgg                                                                                                >  1:613537/1‑148 (MQ=255)
                                                    ccttGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACATTACa                                                                                            <  2:495958/147‑1 (MQ=255)
                                                     tttGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAg                                                                                           >  1:43587/1‑149 (MQ=255)
                                                     tttGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAg                                                                                           <  2:196669/149‑1 (MQ=255)
                                                              aGACCATTACGTAATTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTgg                                                                                  >  1:83374/1‑149 (MQ=255)
                                                                   aTTACGTAATTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGACTTGTGGGAttt                                                                             >  2:124653/1‑149 (MQ=255)
                                                                          aaTTCCGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCaaa                                                                      <  2:100922/149‑1 (MQ=255)
                                                                          aaTTCAGTTAAGATAGATTCAAAATTTACCTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCaaa                                                                      <  2:424305/149‑1 (MQ=255)
                                                                                gTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGc                                                                <  1:594704/149‑1 (MQ=255)
                                                                                   aaGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTTTTTCTTTTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGGGGGATTTGCTCAAAAAGTGCttt                                                             >  2:535851/1‑149 (MQ=255)
                                                                                    agatagatTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTg                                                            >  2:236507/1‑149 (MQ=255)
                                                                                     gatagatTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTgg                                                           <  2:535987/149‑1 (MQ=255)
                                                                                         gatTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCt                                                       <  2:375659/149‑1 (MQ=255)
                                                                                              aaaaTTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGGTTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCtc                                                  >  1:291925/1‑148 (MQ=255)
                                                                                              aaaaTTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAACCGTGCTTTGGTGCGTTCTa                                                  >  1:56069/1‑149 (MQ=255)
                                                                                              aaaaTTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCGGTCTa                                                  >  1:291851/1‑149 (MQ=255)
                                                                                                     aaCTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTgcgc                                           <  1:629592/149‑1 (MQ=255)
                                                                                                     aaCTTCTGAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTgcgc                                           <  2:720771/149‑1 (MQ=255)
                                                                                                        ttcttAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATAACGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCcag                                        <  2:404305/149‑1 (MQ=255)
                                                                                                           ttAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGGTTGGCTCAAACCGTGCTTTGGTTCTCTCTACCTGCGCcagcag                                     >  1:399167/1‑149 (MQ=255)
                                                                                                           ttAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCcagcag                                     >  2:456493/1‑149 (MQ=255)
                                                                                                               aTACAGAATTATTCGGGTAAAACCTACTCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTa                                 <  1:182443/149‑1 (MQ=255)
                                                                                                               aTACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTa                                 <  1:79779/149‑1 (MQ=255)
                                                                                                                tACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGGTTTGGTGTTGTCTACCTGCGCCagcagcagat                                >  2:148843/1‑145 (MQ=255)
                                                                                                                  cAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCg                              >  1:247293/1‑149 (MQ=255)
                                                                                                                        gTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGtttt                        >  1:99119/1‑149 (MQ=255)
                                                                                                                        gTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGtttt                        >  1:34718/1‑149 (MQ=255)
                                                                                                                               cTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCga                 <  2:521637/149‑1 (MQ=255)
                                                                                                                                   aaaTTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCgatgat             <  2:725725/149‑1 (MQ=255)
                                                                                                                                        aCGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGAGTGGTGGGATTTGCTCAAACCGTGATTTGTTTCTGTCTACTTGCGCCAGCAGCGTATCAGGCTTTTCAGCGATAGtttttt        >  2:114291/1‑149 (MQ=255)
                                                                                                                                              tttttCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCgg  <  2:429234/149‑1 (MQ=255)
                                                                                                                                              tttttCCGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCgg  <  1:90473/149‑1 (MQ=255)
                                                                                                                                              tttttCCGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGCGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCgg  <  1:556520/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
CACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAAATTTAACTTCTTAAATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGG  >  NC_000913/4580686‑4580976

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑