Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,223,596 Δ1 bp coding (33/825 nt) iclR ← DNA‑binding transcriptional repressor IclR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,223,5930T.95.8% 82.6 / 1.5 25coding (36/825 nt)iclRDNA‑binding transcriptional repressor IclR
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (0/23);  total (2/23)
Fisher's exact test for biased strand distribution p-value = 4.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGT  >  NC_000913/4223461‑4223713
                                                                                                                                    |                                                                                                                        
gggTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGt                                                                                                                          >  2:320900/1‑133 (MQ=255)
gggTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGt                                                                                                                          >  1:279932/1‑133 (MQ=255)
      cccGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGtttctgc                                                                                                                    >  1:409243/1‑129 (MQ=255)
           cTTGTTGCGGTAGTTTCGTGAGTGGCGCTCTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTCTGCCGCGt                                                                                                              <  2:289019/133‑1 (MQ=255)
              gttgCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTCTGCCGCGTTTc                                                                                                           <  2:522855/133‑1 (MQ=255)
                  cgcCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTCTGCCGCGTTTCGCgg                                                                                                       <  2:279932/133‑1 (MQ=255)
                          ccGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTCTGCCGCGTTTCGCGGGAATGGGt                                                                                               <  2:409243/133‑1 (MQ=255)
                                 tGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTCTGCCGCGGTTCGCGGGAAAGGGGGCGGCCa                                                                                        >  1:556693/1‑133 (MQ=255)
                                 tGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTATGCCGCGGTTCCCGGGAATGGTGGCGGCCa                                                                                        >  1:211630/1‑133 (MQ=255)
                                          gCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGtct                                                                               <  1:556951/133‑1 (MQ=255)
                                          gCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGtct                                                                               <  2:428931/133‑1 (MQ=255)
                                           ccATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGtctc                                                                              <  2:418319/133‑1 (MQ=255)
                                            cATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTctct                                                                             <  1:392854/133‑1 (MQ=255)
                                            cATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTctct                                                                             <  2:50301/133‑1 (MQ=255)
                                                 gATTCGGCAAGCCACGCCAGTAATCTCACGCCGCGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCtttttt                                                                        <  2:211630/133‑1 (MQ=255)
                                                   ttCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑GTTCTGCCCCGTTTTGCGGGGATGGGGGCGGCCACGAGAGGCTCTTTTTTCt                                                                      >  2:444207/1‑133 (MQ=255)
                                                                        tttCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTT‑TGCCGCGTTTCGCGGGGAAGGGTGGGACCCTGACAGGCTCTTTTTTTTGTATCGTGGGAAtcattttta                                                 >  2:157067/1‑131 (MQ=255)
                                                                               ccACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTT‑TGCCGCGTTTTGCGGGGATGGGGGGGGCCATGACAGTCTCTTTTTTTTGTTTCGTGGAAATCATTTTCAttttttt                                          >  1:88992/1‑131 (MQ=255)
                                                                               ccACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTCTGTCCCGCTTCGCGGGGAAGGGGGCGACCACCGCAGTCTCTTTTTTTTGTATCGTGGGAATTAATTTCATTTTTAt                                          >  2:82643/1‑133 (MQ=255)
                                                                                cACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTAtt                                         <  1:148362/133‑1 (MQ=255)
                                                                                    cgTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑TGCCGCGGTTCGCGGGGATGGGGGCGGCCCTGACAGTCTCTTTTTTCTGTATCGTGGGAATCATTTTTATTTTTTTTTTTa                                     >  2:521090/1‑133 (MQ=255)
                                                                                      ttAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGc                                   <  2:645294/133‑1 (MQ=255)
                                                                                           gACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATg                              <  1:82643/133‑1 (MQ=255)
                                                                                               gAACCTGTCCCGTCGCTGGTGCGGCGGCAACGTCGGG‑TTTCCGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAAt                          <  2:168076/133‑1 (MQ=255)
                                                                                               gAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAAt                          <  2:82765/133‑1 (MQ=255)
                                                                                               gAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAAt                          <  1:548433/133‑1 (MQ=255)
                                                                                               gAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAAt                          <  1:320900/133‑1 (MQ=255)
                                                                                                    tGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTa                     <  1:444207/133‑1 (MQ=255)
                                                                                                     gTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑TGCCGCGGTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTTTGTATCGTGGGAATCATTTTTATTTTTTGTTTATGTTAATGCGATAGTTac                    >  1:168076/1‑133 (MQ=255)
                                                                                                         aGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactgaa                <  1:157067/133‑1 (MQ=255)
                                                                                                            cGCTGGTGCGGTGGCAACGGCGGGGTTT‑TGCCGCGGTTCGCGGGGATGGGGGCGACCATGACAGTCTCTTTTTTTTGTATCGTGGGAATCATTTTTATTTTTTTTTTTTGCTAATGCAATAGTTactgaactg             >  2:614714/1‑133 (MQ=255)
                                                                                                             gCTGGTGCGGTGGCAACGGCGGGTTTT‑TGCCGCGGTTCCCGGGGGTGGGGGGGGCCCTTGCAGTCTCTTTTTTTTTTTTCGCGGAAATGATTTTTCTTTTTTTTTTTTTTTAATGCCAAATTTtctgcactga            >  2:398917/1‑133 (MQ=255)
                                                                                                                  tGCGGTGGCAACGGCGGG‑TTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGa       <  2:117001/133‑1 (MQ=255)
                                                                                                                  tGCGGTGGCAACGGCGGG‑TTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGa       <  1:398917/133‑1 (MQ=255)
                                                                                                                  tGCGGTGGCAACGGCGGG‑TTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGa       <  1:34450/133‑1 (MQ=255)
                                                                                                                       tGGCAACGGCGGG‑TTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGt  <  1:614714/133‑1 (MQ=255)
                                                                                                                                    |                                                                                                                        
GGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGT  >  NC_000913/4223461‑4223713

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