Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,926,760 A→G Y82H (TAC→CAC)  asnC ← DNA‑binding transcriptional dual regulator AsnC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,926,7600AG75.3% 301.5 / 113.5 154Y82H (TAC→CAC) asnCDNA‑binding transcriptional dual regulator AsnC
Reads supporting (aligned to +/- strand):  ref base A (17/21);  new base G (49/67);  total (66/88)
Fisher's exact test for biased strand distribution p-value = 8.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.48e-01

Only 100 of 159 total aligned reads displayed.

GCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCAT  >  NC_000913/3926711‑3926807
                                                 |                                               
gCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTa                                                 <  1:2237507/50‑1 (MQ=255)
  ttCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGggtggt                                               >  1:253927/1‑50 (MQ=255)
    cAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCt                                             >  2:2979625/1‑50 (MQ=255)
    cAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAACGCGGACGGGTAGt                                               <  1:2642864/48‑1 (MQ=255)
     aGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCtt                                            >  2:3645788/1‑50 (MQ=255)
     aGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCtt                                            >  2:4915971/1‑50 (MQ=255)
     aGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCtt                                            <  1:2993692/50‑1 (MQ=255)
     aGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCtt                                            <  1:4701323/50‑1 (MQ=255)
     aGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCtt                                            <  1:5662481/50‑1 (MQ=255)
     aGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCtt                                            >  2:5552619/1‑50 (MQ=255)
     aGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCtt                                            <  2:1236181/50‑1 (MQ=255)
     aGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCt                                             <  2:3428755/49‑1 (MQ=255)
     aGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCt                                             <  2:503161/49‑1 (MQ=255)
      gTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCttt                                           <  1:1602197/50‑1 (MQ=255)
      gTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCttt                                           <  1:849154/50‑1 (MQ=255)
      gTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCttt                                           <  2:4617932/50‑1 (MQ=255)
       tGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCTTTg                                          <  1:248529/50‑1 (MQ=255)
       tGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCTTTg                                          <  2:4571323/50‑1 (MQ=255)
       tGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCTTTg                                          >  2:1894398/1‑50 (MQ=255)
       tGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTg                                          <  1:5320429/50‑1 (MQ=255)
        gACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCTTTgg                                         >  1:4999770/1‑50 (MQ=255)
         aCTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCTTTGgc                                        <  1:1954935/50‑1 (MQ=255)
          cTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCTTTGgcg                                       <  1:4086062/50‑1 (MQ=255)
          cTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCTTTGgcg                                       >  1:3765399/1‑50 (MQ=255)
          cTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCTTTGgcg                                       >  1:2652558/1‑50 (MQ=255)
           ttCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCTTTGgcgc                                      >  2:390051/1‑50 (MQ=255)
           ttCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCTTTGgcgc                                      >  2:1116916/1‑50 (MQ=255)
           ttCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCTTTGgcgc                                      >  1:4435779/1‑50 (MQ=255)
             catcaAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCTTTGgcgc                                      <  1:2299373/48‑1 (MQ=255)
             catcaAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCTTTGGCGctc                                    >  1:1810631/1‑50 (MQ=255)
             catcaAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCTTTGGCGctc                                    >  1:829254/1‑50 (MQ=255)
             catcaAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCTTTGGCGctc                                    >  2:1988972/1‑50 (MQ=255)
             catcaAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCTTTGGCGct                                     <  2:5168906/49‑1 (MQ=255)
             catcaAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCTTTGGCGct                                     <  2:4249448/49‑1 (MQ=255)
             catcaAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCTTTGGCGct                                     <  2:5405370/49‑1 (MQ=255)
             catcaAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGctc                                    >  2:2522438/1‑50 (MQ=255)
              atcaAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCTTTGGCGctct                                   <  1:4496999/50‑1 (MQ=255)
              atcaAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCTTTGGCGctc                                    <  2:4324336/49‑1 (MQ=255)
              atcaAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCTTTGGCGctc                                    <  2:587465/49‑1 (MQ=255)
              atcaAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGctct                                   <  2:4071443/50‑1 (MQ=255)
              atcaAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGctct                                   <  1:64008/50‑1 (MQ=255)
              atcaAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGctc                                    <  1:1055190/49‑1 (MQ=255)
               tcaAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCTTTGGCGCTCtt                                  >  2:1384482/1‑50 (MQ=255)
               tcaAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCtt                                  >  2:2059234/1‑50 (MQ=255)
               tcaAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCtt                                  <  1:22770/50‑1 (MQ=255)
                caAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCTTTGGCGCTCtt                                  >  2:838210/1‑49 (MQ=255)
                caAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCTTTGGCGCTCtt                                  >  2:933672/1‑49 (MQ=255)
                 aaGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCTTTGGCGCTCttt                                 <  2:2410586/49‑1 (MQ=255)
                 aaGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCTTTGGCGCTCTTTa                                <  1:3855137/50‑1 (MQ=255)
                  aGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCTTTGGCGCTCTTTaa                               >  1:4872862/1‑50 (MQ=255)
                  aGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCTTTGGCGCTCTTTaa                               <  1:3848534/50‑1 (MQ=255)
                  aGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTaa                               <  1:2507354/50‑1 (MQ=255)
                   ggCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCTTTGGCGCTCTTTaa                               >  2:2312915/1‑49 (MQ=255)
                   ggCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCTTTGGCGCTCTTTAat                              <  1:3565045/50‑1 (MQ=255)
                    gCTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCTTTGGCGCTCTTTAata                             <  2:3761070/50‑1 (MQ=255)
                     cTTTCCAGCTTTGCCAGCGCGGAAGGGTGGTCTTTGGCGCTCTTTAatat                            <  1:1995588/50‑1 (MQ=255)
                       ttCCAGCTTTGCCAGCGCGGAAGGGTGGTCTTTGGCGCTCTTTAATATaa                          >  1:5525478/1‑50 (MQ=255)
                        tCCAGCTTTGCCAGCGCGGAAGGGTGGTCTTTGGCGCTCTTTAATATaa                          <  1:3714117/49‑1 (MQ=255)
                         ccAGCTTTGCCAGCGCGGAAGGGTGGTCTTTGGCGCTCTTTAATATAATg                        >  1:1141449/1‑50 (MQ=255)
                         ccAGCTTTGCCAGCGCGGAAGGGTGGTCTTTGGCGCTCTTTAATATAATg                        >  1:5210873/1‑50 (MQ=255)
                         ccAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAAt                         <  1:2843053/49‑1 (MQ=255)
                          cAGCTTTGCCAGCGCGGAAGGGTGGTCTTTGGCGCTCTTTAATATAATg                        <  1:4234564/49‑1 (MQ=255)
                          cAGCTTTGCCAGCGCGGAAGGGTGGTCTTTGGCGCTCTTTAATATAATGc                       >  1:5058941/1‑50 (MQ=255)
                           aGCTTTGCCAGCGCGGAAGGGTGGTCTTTGGCGCTCTTTAATATAATGcc                      >  2:2752763/1‑50 (MQ=255)
                           aGCTTTGCCAGCGCGGAAGGGTGGTCTTTGGCGCTCTTTAATATAATGcc                      <  1:4389508/50‑1 (MQ=255)
                            gCTTTGCCAGCGCGGAAGGGTGGTCTTTGGCGCTCTTTAATATAATGCCg                     <  2:4093498/50‑1 (MQ=255)
                               ttGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGAt                   >  1:672636/1‑49 (MQ=255)
                                tGCCAGCGCGGAAGGGTGGTCTTTGGCGCTCTTTAATATAATGCCGATa                  <  2:1659914/49‑1 (MQ=255)
                                tGCCAGCGCGGAAGGGTGGTCTTTGGCGCTCTTTAATATAATGCCGATa                  <  2:1147931/49‑1 (MQ=255)
                                tGCCAGCGCGGAAGGGTGGTCTTTGGCGCTCTTTAATATAATGCCGATa                  <  2:2435814/49‑1 (MQ=255)
                                tGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATa                  <  2:5157024/49‑1 (MQ=255)
                                 gCCAGCGCGGAAGGGTGGTCTTTGGCGCTCTTTAATATAATGCCGATaa                 <  1:4712407/49‑1 (MQ=255)
                                  ccAGCGCGGAAGGGTGGTCTTTGGCGCTCTTTAATATAATGCCGATaaa                <  2:4768368/49‑1 (MQ=255)
                                  ccAGCGCGGAAGGGTGGTCTTTGGCGCTCTTTAATATAATGCCGATaaa                <  1:1553565/49‑1 (MQ=255)
                                  ccAGCGCGGAAGGGTGGTCTTTGGCGCTCTTTAATATAATGCCGATAAag               <  2:5563168/50‑1 (MQ=255)
                                   cAGCGCGGAAGGGTGGTCTTTGGCGCTCTTTAATATAATGCCGATAAag               >  2:1755880/1‑49 (MQ=255)
                                   cAGCGCGGAAGGGTGGTCTTTGGCGCTCTTTAATATAATGCCGATAAag               >  2:866534/1‑49 (MQ=255)
                                   cAGCGCGGAAGGGTGGTCTTTGGCGCTCTTTAATATAATGCCGATAAag               >  2:3733772/1‑49 (MQ=255)
                                   cAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAag               >  2:3547266/1‑49 (MQ=255)
                                    aGCGCGGAAGGGTGGTCTTTGGCGCTCTTTAATATAATGCCGATAAag               >  2:5694497/1‑48 (MQ=255)
                                    aGCGCGGAAGGGTGGTCTTTGGCGCTCTTTAATATAATGCCGATAAag               >  2:5380628/1‑48 (MQ=255)
                                     gcgcGGAAGGGTGGTCTTTGGCGCTCTTTAATATAATGCCGATAAagcag            <  2:3339463/50‑1 (MQ=255)
                                     gcgcGGAAGGGTGGTCTTTGGCGCTCTTTAATATAATGCCGATAAagcag            <  1:3649151/50‑1 (MQ=255)
                                     gcgcGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAagcag            <  1:1405173/50‑1 (MQ=255)
                                      cgcgGAAGGGTGGTCTTTGGCGCTCTTTAATATAATGCCGATAAagcagc           <  2:5588957/50‑1 (MQ=255)
                                       gcgGAAGGGTGGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGcc          <  1:1319315/50‑1 (MQ=255)
                                         ggAAGGGTGGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTa        >  1:3184863/1‑50 (MQ=255)
                                         ggAAGGGTGGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTa        >  1:1762392/1‑50 (MQ=255)
                                         ggAAGGGTGGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTa        >  1:1714171/1‑50 (MQ=255)
                                         ggAAGGGTGGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTa        <  1:5263008/50‑1 (MQ=255)
                                         ggAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTa        <  1:5228315/50‑1 (MQ=255)
                                          gAAGGGTGGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTAc       >  1:5628001/1‑50 (MQ=255)
                                          gAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTAc       <  2:1609193/50‑1 (MQ=255)
                                            cgggTGGTCTTTGGCGCTCTTTAATATAATGCAGATAAAGCAGCCTACgg     >  2:2053574/2‑49 (MQ=255)
                                            aGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGt     >  2:1605714/1‑50 (MQ=255)
                                             gggTGGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGt     >  2:5291142/1‑49 (MQ=255)
                                             gggTGGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTc    <  2:1403541/50‑1 (MQ=255)
                                             gggTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGt     >  2:1695170/1‑49 (MQ=255)
                                               gTAGTCTTTGGCGCTCTTTAATGTAATGCCGATAAAGCAGCCTACGTCAt  <  1:3928155/50‑1 (MQ=255)
                                                tAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCAt  >  1:954697/1‑49 (MQ=255)
                                                 |                                               
GCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCAT  >  NC_000913/3926711‑3926807

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