Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,585,505 C→T L419L (CTG→CTA hsdR ← type I restriction enzyme EcoKI endonuclease component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,585,5050CT98.3% 196.3 / NA 59L419L (CTG→CTAhsdRtype I restriction enzyme EcoKI endonuclease component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (28/30);  minor base G (1/0);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGG  >  NC_000913/4585373‑4585638
                                                                                                                                    |                                                                                                                                     
ccccTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCtt                                                                                                                                >  1:646530/1‑140 (MQ=255)
 cccTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTc                                                                                                                               >  1:595220/1‑140 (MQ=255)
     tGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCggg                                                                                                                           <  1:178575/140‑1 (MQ=255)
          aTTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTg                                                                                                                      <  1:575297/140‑1 (MQ=255)
              ccTTTTCAACCGCGCGGTCGGCTTCTTCCTGCTAAAAACCCAGGCCCAACTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCAc                                                                                                                  <  2:279170/140‑1 (MQ=255)
              ccTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCAc                                                                                                                  <  2:846559/140‑1 (MQ=255)
              ccTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCAc                                                                                                                  >  2:837938/1‑140 (MQ=255)
              ccTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCACTGCCGAGCATTTCCACTAGCACTTCCGGGCGGTGCCAc                                                                                                                  >  2:430068/1‑140 (MQ=255)
                     aaCCGCGCGGACGGCGTCTTCCTGACAATGACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTa                                                                                                           <  1:430068/140‑1 (MQ=255)
                      aCCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCTTTTCCAGGAGCGCTGCCGGGCGGTGCCACTCGGGTaa                                                                                                          >  1:279170/1‑140 (MQ=255)
                       ccGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAg                                                                                                         <  2:595220/140‑1 (MQ=255)
                        cgcgcgGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAgg                                                                                                        >  2:966994/1‑140 (MQ=255)
                         gcgcgGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACt                                                                                                                 <  1:868369/130‑1 (MQ=255)
                         gcgcgGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACt                                                                                                                 >  2:868369/1‑130 (MQ=255)
                         gcgcgGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGc                                                                                                       >  2:708196/1‑140 (MQ=255)
                         gcgcgGACGGCATCTACCTGATCATCACGCAGGCCCAACTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGc                                                                                                       <  1:269072/140‑1 (MQ=255)
                           gcgGACGGCATCTTCCTTATAATAACGCAGGCCCACCTCTCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCtt                                                                                                     <  2:296811/140‑1 (MQ=255)
                           gcgGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCtt                                                                                                     <  1:966994/140‑1 (MQ=255)
                                cGGCATCTACCTGATAATAACGAAGGCCCAGTTTGATCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGAc                                                                                                <  2:498456/140‑1 (MQ=255)
                                  gCATCTTCCTGATAAAAACGCAGGCACAGTTCGCTCAGGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACat                                                                                              <  1:331960/140‑1 (MQ=255)
                                             ataataACGCAGGCCTAGCTCGCTCCTGCCAGGGTTACCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTAt                                                                                   <  1:263693/140‑1 (MQ=255)
                                                 taaCGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATcacg                                                                               <  2:91744/140‑1 (MQ=255)
                                                 taaCGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATcacg                                                                               <  1:228482/140‑1 (MQ=255)
                                                 taaCGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTGTCGACAAATTGCGGGTATcacg                                                                               >  1:319161/1‑140 (MQ=255)
                                                      cAGGCCCAACTCGCTCTTGCCATGGTTATCGTCAAACCAATGATTCTGATTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGt                                                                          <  2:465091/140‑1 (MQ=255)
                                                         gCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATcacg                                                                               <  2:76822/132‑1 (MQ=255)
                                                         gCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATcacg                                                                               >  1:76822/1‑132 (MQ=255)
                                                         gCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTcgc                                                                       >  2:766176/1‑140 (MQ=255)
                                                         gCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGCTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTcgc                                                                       >  2:240650/1‑140 (MQ=255)
                                                                          cAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCtt                                                      <  2:788309/140‑1 (MQ=255)
                                                                             ggTTATCGGCAAACCACCGATTCTGTTTTTGTGGTTCGCAGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTg                                                   <  2:319161/140‑1 (MQ=255)
                                                                                  tCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCtt                                              >  2:894295/1‑140 (MQ=255)
                                                                                  tCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCtt                                              >  2:586765/1‑140 (MQ=255)
                                                                                  tCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCtt                                              >  2:902616/1‑140 (MQ=255)
                                                                                  tCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACACATTGCGGGGATCGCGCACGTCGCGACACCAAATGCCGCTTTTGGTCtt                                              >  2:295432/1‑140 (MQ=255)
                                                                                    ggCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCa                                            <  1:240650/140‑1 (MQ=255)
                                                                                    ggCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATAATGCGGGTATCACCCACATCGCGATACCAGACGCCGCTTTTGGTCTTCa                                            >  2:418481/1‑140 (MQ=255)
                                                                                       aaaCCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCAttg                                         >  2:702809/1‑140 (MQ=255)
                                                                                       aaaCCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCAttg                                         <  2:709507/140‑1 (MQ=255)
                                                                                        aaCCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCtt                                                      <  1:783810/126‑1 (MQ=255)
                                                                                        aaCCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCtt                                                      >  2:783810/1‑126 (MQ=255)
                                                                                        aaCCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCAttgt                                        >  2:543830/1‑140 (MQ=255)
                                                                                        aaCCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCAttgt                                        >  2:612650/1‑140 (MQ=255)
                                                                                         aCCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCAttgtt                                       <  1:702809/140‑1 (MQ=255)
                                                                                             cTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTgcgc                                   <  1:766176/140‑1 (MQ=255)
                                                                                                  tCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTAt                              <  1:213520/140‑1 (MQ=255)
                                                                                                     gTTTTTGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTcg                           >  1:842226/1‑140 (MQ=255)
                                                                                                           gCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCg                     >  2:582194/1‑140 (MQ=255)
                                                                                                            cGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGt                    <  1:586765/140‑1 (MQ=255)
                                                                                                              gTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGGTGCGAGGTATTCGCGCCCGTTg                  >  2:407115/1‑140 (MQ=255)
                                                                                                              gTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATAACAGATGCCGCTTATGGGCTTCATTGTTGCGCGGTATTCGCGCCCGTTg                  >  2:431000/1‑140 (MQ=255)
                                                                                                                tCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGt                <  2:646530/140‑1 (MQ=255)
                                                                                                                tCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGt                >  2:17213/1‑140 (MQ=255)
                                                                                                                tCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGt                <  1:113906/140‑1 (MQ=255)
                                                                                                                         aGCATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCa       >  2:886236/1‑140 (MQ=255)
                                                                                                                           cATTTCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGa     >  2:535818/1‑140 (MQ=255)
                                                                                                                              ttCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAgg  <  2:1002075/140‑1 (MQ=255)
                                                                                                                              ttCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAgg  <  1:582194/140‑1 (MQ=255)
                                                                                                                              ttCCAGTAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGATCTTCATTGTTGGGCGTTATTCGCGACCGTAGGTCGAGTAGCAGAAgg  >  2:291384/1‑140 (MQ=255)
                                                                                                                                    |                                                                                                                                     
CCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGG  >  NC_000913/4585373‑4585638

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