Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,223,217 C→T V138I (GTA→ATA)  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,2170CT84.6% 16.1 / ‑2.7 14V138I (GTA→ATA) iclRDNA‑binding transcriptional repressor IclR
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (7/4);  total (8/6)
Fisher's exact test for biased strand distribution p-value = 1.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

GTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTACCCGCACCGGAAGCGTGCATCGGCAATTTACCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAGACTCTTCCATTAGATTGCGCAGGATAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTG  >  NC_000913/4222993‑4223461
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
ntCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTACCCGCACCGGAAGCGTGCATCGGCAATTTACCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTATCTGGTCGATAATAATCGCTTCGTGAt                                                                                                                                                                                                                            >  1:201579/2‑251 (MQ=255)
          ttcttccgatctcgcgTTTGGGCGAGATCTTTTTTTAAATCCACATGAGACTCCAGCGTTGCGTGGGTAAAGCCAGGTACCCCTTTGCGTTGCGGCAGCTTCTTCCCTTTTTCTTCCCTAAGTGGGGTAAAAATGTCTTTCCCCGCACCGGAGGCGTGCATCGACAATTTACCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTATCTGGTCGATAATAATCGCTTCGTGATCGCTTTGAt                                                                                                                                                                                                                   <  1:144876/238‑1 (MQ=255)
                 gTTTGCGCGTTTGGGCGCGATCTTCTTTTAAATGCACCGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACACTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTAACCGCCCCGGAAGCGTGCATCGGCAATTTACCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTATCTGGCCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCAcc                                                                                                                                                                                                           >  1:240979/1‑251 (MQ=255)
                                                                                                              tcggtcccttttcttccccgtggtGGGCAAAAAAAGTTTTACCCGCACCGGAAGGCGTGACTGGCAATTTACCCGCCGATAGGCGCGGACATTCGCATCAGGTGCTTACACTTTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATAAAGCCCCGCCATATTGACCGTTTCGCCAGACTCTTCCATTAGATGGCGCAGGATAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAgctgc                                                                                                                <  1:27813/227‑1 (MQ=255)
                                                                                                                                 cNGTTGGGCTACAACGGCTTTACCCGCACCGGAAGCGTGCATCGGCAATTAACCGCCGATAGGCGCGGACATTCGCCTCAGATGCGTACACTGTATCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAGACTCTTCCATTAGATTGCGCAGGATAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGACGACCATAAAGGCAtgtg                                                                                           >  1:255013/1‑251 (MQ=255)
                                                                                                                                   gTTGGGATAAAAAGGCTTTACCCGCACCGGAAGCGTGCATCGGCAATTTACCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTATCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAGACTCTTCCATTAGATTGCGCAGGATAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTgcg                                                                                         >  1:72545/1‑251 (MQ=255)
                                                                                                                                                                    tttcccaatttaaatttccccgccgAGAGGCGCGGACATTCGCAGCAGATGCGTACACTGTATCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAGACTCTTCCATTAGATTGCGCAGGATAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAAcc                                                           <  1:91387/231‑1 (MQ=255)
                                                                                                                                                                     cggtccttcgtaattttcccgccgATAGGCGCGGACATGCGCATCAGATGCGTACACTGTACCTGGTGAATAATAATCGCTTCGTGATCTCTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAGACTCGTCCATTAGATTGCGCAGGATAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTg                                                         <  1:182067/233‑1 (MQ=255)
                                                                                                                                                                               gCAATTTACCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTAACTGGTCGATAATACTCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAGACTCTTCCATTAGATTGCGCAGGATAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGaa                                                >  1:126301/1‑250 (MQ=255)
                                                                                                                                                                                           ccgATAGGCGCGGACATTCGCATCAGATGCGTACACTGTATCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAGACTCTTCCATTAGATTGCGCCGGATCGGGTGAACAATCGCTAACAAATTACGGctc                                                                                                                             >  1:45084/1‑162 (MQ=255)
                                                                                                                                                                                                              gCATCAGATGCGTCCACTGTATCTGGTCGATAAAAATCTCTTCGTGACCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAGACTCTTCCATTAGATTGCGCAGGATAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCgn                                                                <  1:194052/204‑2 (MQ=255)
                                                                                                                                                                                                              gCATCAGATGCGTACACTGTATCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAGACTCTTCCATTAGATTGCGCAGGATAGGGTGAACAATCGCTAACAANTTANNNnnctnn                                                                                                                          >  2:194052/1‑144 (MQ=255)
                                                                                                                                                                                                                          tACACTGTATCTGGTCGATCATACTAGCTTCGTGATCGCTTTGATCAAGCACCGCCAAATTGACCGTTTCGCCAGACTCTTCCATTAGATTGCGCAGGATAGGGTGAACAATCGCTCACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGa      >  1:99074/1‑250 (MQ=255)
                                                                                                                                                                                                                             ccTGTATCTGGTCGATAATAATCGCTTCGTGATCGTTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAGACTCTTCCATTAGATTGCGCAGGATAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTg  <  1:298677/250‑1 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
GTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTACCCGCACCGGAAGCGTGCATCGGCAATTTACCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAGACTCTTCCATTAGATTGCGCAGGATAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTG  >  NC_000913/4222993‑4223461

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: