Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,555,047 2 bp→TG coding (1963‑1964/2706 nt) malT → DNA‑binding transcriptional activator MalT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,555,0470GT100.0% 35.5 / NA 14V655C (GTC→TGC) malTDNA‑binding transcriptional activator MalT
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (6/8);  total (6/8)
*NC_0009133,555,0480TG92.9% 33.3 / NA 15V655C (GTC→TGC) malTDNA‑binding transcriptional activator MalT
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (5/8);  total (6/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

GGGCGCGGCTGGATGAAGCCGAAGCGTCGGCGCGTAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCGGAGTCTGCGGTTGATG  >  NC_000913/3554832‑3555256
                                                                                                                                                                                                                       ||                                                                                                                                                                                                                
ccttccgatctgATGAAGCCGAAGCGCCGGCGCGTAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCTATGTTGATTCAATGCTCGCGGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTGGGGAAAACCTGCTGGGGAATGGCAAATATCACTGCGACTGGATCTCTAACGCCATAAAATGCCGGGTGCTTTACTGGCAAATGCCCGGCGATaa                                                                                                                                                                                 <  1:138178/239‑1 (MQ=255)
      ggCTGGATGAAGCCGAAGCGTCGGCGCGTAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGCTAAAg                                                                                                                                                                               >  1:9974/1‑246 (MQ=255)
                   cGAAGCGTCGGCGCGTAGCGAGATTGAAGCCTTGTCGTCTTAGCAGCCACAGCAACAGCTTCAGTGCTTGTCAAGGGTGATTCAATGCCCGCTGGCCCGTGTTGATTTAGATAACGCCCGTATCCAGCTGAACCATCTGGAAACCCGGCTGGGGAATGGCAAATATCACAGCCACTGTATCTCTATCGCCAACAAATTCCGTGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGt                                                                                                                                                             <  1:163456/251‑1 (MQ=255)
                                                                                     ccTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAgc                                                                                                                                                                              >  1:79436/1‑168 (MQ=255)
                                                                                     ccTGGCAATGTTGANNCAATGCTCGNTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACNTGCTGGGGANNGGNAAATANCNNNGNNNCTNGATCTCTAACGCCAACAAANGCNNGGTGATNTACTGGCAAATGACCGGCGATAAAgc                                                                                                                                                                              <  2:79436/168‑1 (MQ=255)
                                                                                           aaTGGTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGTATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCTTCATACGGCTAAACCATAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTg                                                                                     <  1:226754/251‑1 (MQ=255)
                                                                                            aTGTTGATTCAATGCTCGATGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGc                                                                                    <  1:25179/251‑1 (MQ=255)
                                                                                                                                gATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTCACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAAc                                                                                                                                           >  1:143490/1‑160 (MQ=255)
                                                                                                                                                               aaaaCCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATNg                                              <  1:275871/222‑1 (MQ=255)
                                                                                                                                                                  aCCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATACCGGGTGCTTTACTGGCAAATGACCGGCGATAAAGCCGCCCCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTcc                                                                                                                   >  1:145504/1‑150 (MQ=255)
                                                                                                                                                                       ctgGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATACCCGGAGTCTGCg         >  1:112859/1‑251 (MQ=255)
                                                                                                                                                                              aTGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCGGAGTCTGCGGTtgatn  <  1:90387/251‑2 (MQ=255)
                                                                                                                                                                              aTGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAAAAATTGTGCTCGAAGAACTCAATGAAAATGCCCGGAGTCTGCGGTtgatn  <  1:63193/251‑2 (MQ=255)
                                                                                                                                                                                 gCAAATATCACAGCGACTGGATCTCTAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAANNGGTTGCNNc                                                                                                                                                       >  2:21617/1‑99 (MQ=255)
                                                                                                                                                                                                        tctAACGCCAACAAATGCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTGTGCGAACAACCACTTCCTGCAAGGTCAAGGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTtc                                           <  1:227401/184‑1 (MQ=255)
                                                                                                                                                                                                                       ||                                                                                                                                                                                                                
GGGCGCGGCTGGATGAAGCCGAAGCGTCGGCGCGTAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTCCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCGGAGTCTGCGGTTGATG  >  NC_000913/3554832‑3555256

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