Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,604,517 C→T L120L (CTG→TTG)  bglJ → DNA‑binding transcriptional regulator BglJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,5170CT96.7% 99.8 / ‑6.6 30L120L (CTG→TTG) bglJDNA‑binding transcriptional regulator BglJ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (12/17);  minor base G (1/0);  total (13/17)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGT  >  NC_000913/4604286‑4604749
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                        
aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGAGTTTTCATCAGGa                                                                                                                                                                                                                       <  2:173298/251‑1 (MQ=255)
 gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGAATTTTCATCAGGaa                                                                                                                                                                                                                      <  1:201194/251‑1 (MQ=255)
          cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTg                                                                                                                                                                                                             <  1:332004/251‑1 (MQ=255)
           aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGGATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGt                                                                                                                                                                                                            <  1:157961/251‑1 (MQ=255)
                aGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTAATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGGAGAGAGGGAGTATCTTGCCTGACTGAACTGGGGGTTAAGTTTTCGCGTACCCGGGGGTTAGTTATTTCGGGTGGTGATATTTAAGCTTGGCTGATTGGGTTATTTGCGGCATCACCGCTGGTAGGTGTGTTAAGTAAAACGTCACAGTTGGGGAGTTTTCTTTAGGAACTTTTTTTTTCAttt                                                                                                                                                                                                       >  1:63914/1‑250 (MQ=255)
                  gATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCCTGGGGAGGGGGCGCAGAGGGGGGTTATCTTGCCTTTCTGTACTGGGGGTTAAGGTTTCGCGTACCCCGCGTTTTGTTTTTGCGGGGGAGGGTTTTTTAGCTCGGCTGGTTGGTTTTTTGTCGCCCTCACCGCTGGGAGGGGTATTAAGTAAAGCGGCAAAGGTGGGGGTTTTTCTTCCGGAACCCtttttttctttta                                                                                                                                                                                                      >  2:200442/1‑242 (MQ=255)
                              cgGATCTCGTGTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTGTGTCATTAAATGGTGTACGTCa                                                                                                                                                                                         <  2:40750/251‑1 (MQ=255)
                                             gcggTCATTTTGTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACt                                                                                                                                                                          <  2:97793/251‑1 (MQ=255)
                                             gaggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATGATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATAGGTTCATTGTCGCCATCACCGCTTGACGGTGTATTAAGTAAAGCGTCAACTTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCTACt                                                                                                                                                                          <  1:121141/249‑1 (MQ=255)
                                                          ctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGACCTGGCGATTAAGTTTCCGAGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGAGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTgg                                                                                                                                                             >  1:34710/1‑251 (MQ=255)
                                                            ttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGt                                                                                                                                                            >  1:283512/1‑250 (MQ=255)
                                                               tctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACa                                                                                                                                                         <  1:138281/250‑1 (MQ=255)
                                                                     aTGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAAcc                                                                                                                                                   >  1:272346/1‑250 (MQ=255)
                                                                                     gagagGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCACAGCCGGACGTTAAGc                                                                                                                                  >  1:232719/1‑251 (MQ=255)
                                                                                         gggATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGa                                                                                                                              <  2:93163/251‑1 (MQ=255)
                                                                                              tATCATGCCTGACTGAACTGGCGATTAAGTTTCCGCGAACCCGGCGTTTAGTTATTGCGGATGATGATATTGGAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACgg                                                                                                                           >  1:241230/1‑249 (MQ=255)
                                                                                                                tGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGCTGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGCACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCttt                                                                                                       >  2:305173/1‑251 (MQ=255)
                                                                                                                                  gTACCCGGCGTTTAGTTATTGCGGATGATGAAATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAAGGAGAATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGAGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAACGACGGACGTTAAGCCAGACGGAGAGAGCAATATTGCGCTTTATGTCGCGTGGCTACTc                                                                                      >  1:311248/1‑250 (MQ=255)
                                                                                                                                    aCCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAAt                                                                                   <  1:341943/251‑1 (MQ=255)
                                                                                                                                         gCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGacac                                                                              >  2:119919/1‑251 (MQ=255)
                                                                                                                                               aGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGTCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTg                                                                        <  1:309930/251‑1 (MQ=255)
                                                                                                                                                   aTTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAAc                                                                                                                                                    >  1:344115/1‑171 (MQ=255)
                                                                                                                                                   aTTGCGGATGATGATATTGAAGCGCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAAc                                                                                                                                                    <  2:344115/171‑1 (MQ=255)
                                                                                                                                                                    tGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTCTTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGAACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAat                                                    >  1:129872/1‑250 (MQ=255)
                                                                                                                                                                           cGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAAc                                            <  1:118680/251‑1 (MQ=255)
                                                                                                                                                                            ggCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACg                                           <  1:200442/251‑1 (MQ=255)
                                                                                                                                                                             gCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACg                                           >  1:148423/1‑250 (MQ=255)
                                                                                                                                                                                         aTTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATa                              <  2:251719/251‑1 (MQ=255)
                                                                                                                                                                                              cgtactCACTGCTGGACGATGTAGTAAGTAAAGCGTCCACGTTGGAGATTTTTCATCAGGAACGCTTTTTTGCATTAAATGTTGTACGGCAGGCGACCGGCCGACTGAACAATCAGTGGGACATTAACCAAAGCCGGACGTTAAGACCGACGGAGAGAGAAAAATTTCACTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAAttt                         <  1:126567/246‑1 (MQ=255)
                                                                                                                                                                                                                      tAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGt  <  1:285123/250‑1 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                        
AGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGT  >  NC_000913/4604286‑4604749

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