Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,586,255 C→A 8.8% A170S (GCG→TCG)  hsdR ← endonuclease R Type I restriction enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,586,2550CA 8.8% 123.2 / 2.6 46A170S (GCG→TCG) hsdRendonuclease R Type I restriction enzyme
Reads supporting (aligned to +/- strand):  ref base C (18/23);  new base A (2/2);  total (20/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 = 8.09e-01

TCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGAC  >  NC_000913/4586111‑4586399
                                                                                                                                                |                                                                                                                                                
tCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCg                                                                                                                                              >  1:24220/1‑149 (MQ=255)
 gtGTTTTGCCTGTTCCGCGTTCTTCTCGGCGGGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGa                                                                                                                                             <  2:239019/148‑1 (MQ=255)
  tGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGAc                                                                                                                                            >  2:516146/1‑149 (MQ=255)
         ccAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGcc                                                                                                                                     >  1:248259/1‑149 (MQ=255)
          cAGTTCCGCGTTCTGCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCt                                                                                                                                    <  2:388354/149‑1 (MQ=255)
                 gcgTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGGTGTTGCGCTTCGCCTTTTGCCTGAGTCtc                                                                                                                             >  1:204188/1‑148 (MQ=255)
                    ttcttcTCGGCGAGCTGTGCTTCTAGTGCCGCAAGGCGCGGGTGGTGTGACGCGCCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTgcgc                                                                                                                          <  1:329743/149‑1 (MQ=255)
                      cttctCGGCGGGCTGTCCTTCCAGTGCCCCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCtt                                                                                                                        <  1:220158/149‑1 (MQ=255)
                      cttctCGGCGAGCTGGGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCATAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCtt                                                                                                                        <  1:256987/149‑1 (MQ=255)
                        tctcGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCGGTTGCGCTTCGACTTCTGCCGGAGTCGGGCCtttt                                                                                                                      >  1:232485/1‑149 (MQ=255)
                                     tgCTTCCCGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACctgct                                                                                                         <  1:120537/149‑1 (MQ=255)
                                            aGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGGGTGTTCGCGCACCTGCTGTTCAAg                                                                                                  >  1:568945/1‑149 (MQ=255)
                                              tGCCGCAAGGCGGGCCGGGGTTTGGCCTTCGGTTGCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCt                                                                                                <  1:3881/149‑1 (MQ=255)
                                                ccGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCtgt                                                                                              <  1:264545/149‑1 (MQ=255)
                                                  gCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTACTGCCTGAGTCTGCGCTTTTTCTCGCCCCTGCTGTTCAGGCtgtgg                                                                                            >  1:172290/1‑147 (MQ=255)
                                                      ggCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTGCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGtttt                                                                                        <  2:134798/149‑1 (MQ=255)
                                                       gcgcgcCTGGGTTTCCGCCTCGGTTTCCGGCTGTTTGCTTTCCAGGATGGTGATATAGCCGTTCATGGCCTCCAGCTTCTGCTGTTGCGCTTCGTCTTCTGCCTGAGTCCGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTa                                                                                       <  2:627890/149‑1 (MQ=255)
                                                              tGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAg                                                                                >  2:589179/1‑149 (MQ=255)
                                                                 gTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGATTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCAc                                                                             >  2:445134/1‑149 (MQ=255)
                                                                 gTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGATTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCAc                                                                             >  2:445077/1‑149 (MQ=255)
                                                                       gcTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTg                                                                       <  1:689935/149‑1 (MQ=255)
                                                                       gcTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTg                                                                       <  1:618709/149‑1 (MQ=255)
                                                                        cTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTgg                                                                      <  1:66382/149‑1 (MQ=255)
                                                                                tCCTGCGGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTATTAGCGTCAGCACTTCCTGGTGATagag                                                              <  1:383667/149‑1 (MQ=255)
                                                                                  ctgctgTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAggg                                                            >  2:687138/1‑148 (MQ=255)
                                                                                   tgctgTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGtt                                                           <  2:248259/149‑1 (MQ=255)
                                                                                             cTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTCTTCACCAAGt                                                 >  2:680053/1‑149 (MQ=255)
                                                                                              ttCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGGCTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTTTTCAGGCTGTTGTTTTAGCGTCCGCACTTCCGGGTGATAGAGGGTTTCACCACGtt                                                >  1:57007/1‑149 (MQ=255)
                                                                                                    aaTGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCa                                          <  1:150481/149‑1 (MQ=255)
                                                                                                       ggCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAaca                                       >  1:620171/1‑149 (MQ=255)
                                                                                                              taGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCGCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAacac                                <  1:227870/149‑1 (MQ=255)
                                                                                                              taGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTCTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAacac                                <  1:227880/149‑1 (MQ=255)
                                                                                                                    ttCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAAcaccggcac                          >  1:387894/1‑149 (MQ=255)
                                                                                                                          gCACCCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGaa                    <  2:204188/149‑1 (MQ=255)
                                                                                                                          gCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGCGGTTATCACCACGTTCCGGCAACACAAACACCGGCCCCGGGaa                    >  2:380234/1‑149 (MQ=255)
                                                                                                                            aaCCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACCCCGGCACCGGGAAGt                  >  2:219008/1‑149 (MQ=255)
                                                                                                                            aaCCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAgc                  >  2:47212/1‑148 (MQ=255)
                                                                                                                               cAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCAt               <  1:252580/149‑1 (MQ=255)
                                                                                                                                      tgctgTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTTTCTTTTTCTGTTCAAACAGTGGATTTAGCGTCAACACATCCTGTTGATAGAGGTTTTCACCACGTTCCGGCACCACAAACACCGGCACCGGGGCGCCATAATCgtt        >  2:591862/1‑149 (MQ=255)
                                                                                                                                      atctgTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCttt        <  2:172290/147‑1 (MQ=255)
                                                                                                                                       gctgTTGCGCTTCGCCTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCTTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTa       <  1:319106/149‑1 (MQ=255)
                                                                                                                                       gctgTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACAGCGGCACCGGGAAGTCATAATCTTTa       <  1:413765/149‑1 (MQ=255)
                                                                                                                                         tgTTGCGTTTCGCCTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGt     <  2:571954/149‑1 (MQ=255)
                                                                                                                                           ttGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGa   >  1:104768/1‑149 (MQ=255)
                                                                                                                                            tGCGATTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGAc  <  1:445077/149‑1 (MQ=255)
                                                                                                                                            tGCGATTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGAc  <  1:445134/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
TCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGAC  >  NC_000913/4586111‑4586399

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