Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,360,900 Δ1 bp 10.0% coding (1036/1539 nt) cadC ← cadBA operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,360,9000A.10.0% 118.7 / 7.4 40coding (1036/1539 nt)cadCcadBA operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base A (21/15);  new base . (2/2);  total (23/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGATGTTGAGAATGGCGCACGATATCAACTAATGCTTTTTCTGCTCTCGCGTAGGTAAATTCTGGGGATGATTGAACAATCTCACCTAATAATTCACTGGCACGGTTCAATGATTTATCATCGCCATGCAGTAAATAATCATGTGCCTGATAAAAATTAGTTAATAACGCACCACGATGCGGCAAAATTTTCTGGAGCGTCTCCTGCATTCGTTGTGGCCACGGTTGGTTTAACGCTTTTGATAAACTCTCCAGTAAATCATTTTGAATCGC  >  NC_000913/4360763‑4361035
                                                                                                                                         |                                                                                                                                       
aGGATGTTGAGAATGGCGCACGATATCAACTAATGCTTTTTCTGCTCTCGCGTAGGTAAATTCTGGGGATGATTGAACAATCTCACCTAATAATTCACTGGCACGGTTCAATGATTTATCATCGCCATGCAGTAAATAATCATGTGCCt                                                                                                                              <  1:618526/149‑1 (MQ=255)
         agaATGGCGCACGATATCAACTAATGCTTTTTCTGCTCTCGCGTAGGTAAATTCTGGGGATGATTGAACAATCTCACCTAATAATTCACTGGCACGGTTCAATGATTTATCATCGCCATGCAGTAAATAATCATGTGCCTGATAAAAAt                                                                                                                     <  1:124164/149‑1 (MQ=255)
         agaATGGCGCACGATATCAACTAATGCTTTTTCTGCTCTCGCGTAGGTAAATTCTGGGGATGATTGAACAATCTCACCTAATAATTCACTGGCACGGTTCAATGATTTATCATCGCCATGCAGTAAATAATCATGTGCCTGATAAAAAt                                                                                                                     >  2:704790/1‑149 (MQ=255)
         agaATGGCGCACGATATCAACTAATGCTTTTTCTGCTCTCGAGTAGGTAAATTCTGGGGATGATTGAACAATCTCACCTAATAATTCACTGGCACGGTTCAATGATTTATCATCGCCATGCAGTAAATAATCATGTGCCTGATAAAAAt                                                                                                                     <  1:341311/149‑1 (MQ=255)
             tGGCGCACGATATCAACTAATGCTTTTTCTGCTCTCGCGTAGGTAAATTCTGGGGATGATTGAACAATCTCACCTAATAATTCACTGGCACGGTTCAATGATTTATCATCGCCATGCAGTAAATAATCATGTGCCTGATAAAAATTAGt                                                                                                                 <  1:31214/149‑1 (MQ=255)
             tGGCGAACGATATCAACTAATGCTTTTTCTGCTCTCGCGTAGGTAAATTCTGGGGATGATTGAACAATCTCACCTAATAATTCACTGGCACGGTTCAATGATTTATCATCGCCATGCAGTAAATAATCATGTGCCTGATAAAAATTAGt                                                                                                                 <  1:119253/149‑1 (MQ=255)
              ggCGCACGATATCAACTAATGCTTTTTCTGCTCTCGCGTAGGTAAATGCTGGGGATGATTGAACAATCTCACCTAATAATTCACTGGCACGGTTCAATGATTTATCATCGCCATGCAGTAAATAATCATGTGCCTGATAAAAAATAGtt                                                                                                                >  2:715373/1‑149 (MQ=255)
                  cACGATATCAACTAATGCTTTTTCTGCTCTCGCGTAGGTAAATTCTGGGGATGATTGAACAATCTCACCTAATAATTCACTGGCACGGTTCAATGATTTATCATCGCCATGCAGTAAATAATCATGTGCCTGATAAAAATTAGTtaata                                                                                                            >  1:120274/1‑149 (MQ=255)
                  cACGATATCAACTAATGCTTTTTCTGCTCTCGCGTAGGTAAATTCTGGGGATGATTGAACAATCTCACCTAATAATTCACTGGCACGGTTCAATGATTTATCATCGCCATGCAGTAAATAATCATGTGCCTGATAAAAATTAGTtaata                                                                                                            >  2:447755/1‑149 (MQ=255)
                      atatCAACTAATGCTCTTTCTGCTCTCGCGTAGGTAAATTCTGGGGATGATTGAACAATCTCACCTAATAATTCACTGGCACGGTTCAATGATTTATCATCGCCATGCAGTAAATAATCATGTGCCTGATAAAAATTAGTTAATAACGc                                                                                                        <  2:527031/149‑1 (MQ=255)
                            aCTAATGCTTTTTCTGCTCTCGCGTAGGTAAATTCTGGGGATGATTGAACAATCTCACCTAATAATTCACTGGCACGGTTCAATGATTTATCATCGCCATGCAGTAAATAATCATGTGCCTGATAAAAATTAGTTAATAACGCACCACg                                                                                                  >  1:759221/1‑149 (MQ=255)
                                 tGCTTTTTCTGCTCTCGCGTAGGTAAATTCTGGGGATGATTGAACAATCTCACCTAATAATTCACTGGCACGGTTCAATGATTTATCATCGCCATGCAGTAAATAATCATGTGCCTGATAAAAATTAGTTAATAACGCACCACGATGCg                                                                                             <  1:610799/149‑1 (MQ=255)
                                 tGCTTTTTCTGCTCTCGCGTAGGTAAAGTCTGGGGATGATTGAACAATCTCACCTAATAATTCACTGGCACGGTTCAATGATTTATCATCGCCACGCAGTAAATAATCATGTGCCTGATAAAAATTAGTTAATAACGCACCACGATGCg                                                                                             <  1:610806/149‑1 (MQ=255)
                                   cTTTTTCTGCTCTCGCGTAGGTAAATTCTGGGGATGATTGAACAATCTCACCTAATAATTCACTGGCACGGTTCAATGATTTATCATCGCCATGCAGTAAATAATCATGTGCCTGATAAAAATTAGTTAATAACGCACCACGATGCGGc                                                                                           >  1:657061/1‑149 (MQ=255)
                                   cTTTTTCTGCTCTCGCGTAGGTAAATTCTGGGGATGATTGAACAATCTCACCTAATAATTCACTGGCACGGTTCAATGATTTATCATCGCCATGCAGTAAATAATCATGTGCCTGATAAAAATTAGTTAATAACGCACCACGATGCGGc                                                                                           >  1:75078/1‑149 (MQ=255)
                                   cTTTTTCTGCTCTCGCGTAGGTAAATTCTGGGGATGATTGAACAATCTCACCTAATAATTCACTGGCACGGTTCAATGATTTATCATCGCCATGCAGTAAAT‑ATCATGTGCCTGATAAAAATTAGTTAATAACGCACCACGATGCGGCa                                                                                          >  2:724908/1‑149 (MQ=255)
                                   cTTTTTCTGCTCTCGCGTAGGTAAATTCTGGGGATGATTGAACAATCTCACCTAATAATTCACTGGCACGGTTCAATGATTTATCATCGCCATGCAGTAAAT‑ATCATGTGCCTGATAAAAATTAGTTAATAACGCACCACGATGCGGCa                                                                                          >  2:650115/1‑149 (MQ=255)
                                           gCTCTCGCGTAGGTAAATTCTGGGGTTGATTGAACAATCTCACCTAATAATTCACTGGCACGGTTCAATGATTTATCATCGCCATGCAGTAAATAATCATGTGCCTGATAAAAATTAGTTAATAACGCACCACGATGCGGCAAAAtttt                                                                                   >  1:392796/1‑149 (MQ=255)
                                                       gTAAATTCTGGGGATGATTGAACAATCTCACCTAATAATTCACTGGCACGGTTCAATGATTTATCATCGCCATGCAGTAAATAATCATGTGCCTGATAAAAATTAGTTAATAACGCACCACGATGCGGCAAAATTTTCTGGAGCGtctc                                                                       <  2:545308/149‑1 (MQ=255)
                                                               tGGGGATGATTGAACAATCTCACCTAATCATTCACTGGCACGGTTCAATGATTTATCATCGCCATGCAGTAAATAATCATGTGCCTGATAATAATTAGTTAATAACGCACCACGATGCGGCAAAATTTTCTGGAGCGTCTCCTGCATTc                                                               >  1:419645/1‑149 (MQ=255)
                                                               tGGGGATGATTGAACAATCTCACCTAATAATTCACTGGCACGGTTCAATGATTTATCATCGCCATGCAGTAAATAATCATGTGCCTGATAAAAATTAGTTAATAACGCACCACGATGCGGCAAAATTTTCTGGAGCGTCTCCTGCATTc                                                               <  1:753689/149‑1 (MQ=255)
                                                                         tGAACAATCTCACCTAATAATTCACTGGCACGGTTCAATGATTTATCATCGCCATGCAGTAAATAATCATGTGCCTGATAAAAATTAGTTAATAACGCACCACGATGCGGCAAAATTTTCTGGAGCGTCTCCTGCATTCGTTGTGGCCa                                                     >  2:268499/1‑149 (MQ=255)
                                                                           caCAATCTCACCTAACAATTCACTGGCACGGGTCAATGATTTATCATCGCCATGCAGTAAATAATCATGTGCCTGATAAAAATTAGTTAATAACGCACCACGATGCGGCAAAATTTTCTGGAGCGTCTCCTGCATTCGTTGTGGCCACg                                                   >  2:573611/2‑149 (MQ=255)
                                                                           aaCAATCTCACCTAATAATTCACTGGCACGGTTCAATGATTTATCATCGCCATGCAGTAAATAATCATGTGCCTGATAAAAATTAGTTAATAACGCACCACGATGCGGCAAAATTTTCTGGAGCGTCTGCTGCATTAGTTGTGGCCACg                                                   >  2:45682/1‑149 (MQ=255)
                                                                           aaCAATCTCACCTAATAATTCACTGGCACGGTTCAATGATTTATCATCGCCATGCAGTAAATAATCATGTGCCTGATAAAAATTAGTTAATAACGCACCACGATGCGGCAAAATTTTCTGGAGCGTCTCCTGCATTCGTTGTGGCCACg                                                   >  1:114177/1‑149 (MQ=255)
                                                                           aaCAATCTCACCTAATAATTCACAGGCACGGTTCAATGATTTATCATCGCCATGCAGTAAAAAATCATGTGCCTGATAAAAATTAGTTAATAACGCACCACGATGCGGCAAAATTTTCTGAAGCGTCTCCTGCATTCGTTGTGGCCACg                                                   >  1:355859/1‑149 (MQ=255)
                                                                                         ataatTCACTGGCACGGTTCAATGATTTATCATCGCCATGCAGTAAATAATCATGTGCCTGATAAAAATTAGTTAATAACGCACCACGATGCGGCAAAATTTTCTGGCGCGTCTCCTCCATTCGTTGTGGCCAAGGGTGGGTTCACGcg                                     >  2:158360/1‑148 (MQ=255)
                                                                                              tCACTTGCCCGGTTCAATGATTTATCATCGCCATGCAGTAAATAATCATGTGCCTGATAAAAATTAGTTAATAACGCACCACGATGCGGCAAAATTTTCTGGAGCGTCTCCTGCATTCGTTGTGGCCACGGTTGGTTTAACGCTTTTGa                                <  2:117800/149‑1 (MQ=255)
                                                                                               cACTGGCACGGTTCAATGATTTATCATCGCCATGCAGTAAATAATCATGTGCCTGATAAAAATTAGTTAATAACGCACCACGATGCGGCAAAATTTTCTGGAGCGTCTCCTGCATTCGTTGTGGCCACGGTTGGTTTAACGCTTTTGAt                               <  1:469178/149‑1 (MQ=255)
                                                                                                  tGGCACGGTTCAATGATTTATCATCGCCATGCAGTAAATAATCATGTGCCTGATAAAAATTAGTTAATAACGCACCACGATGCGGCAAAATTTTCTGGAGCGTCTCCTGCATTCGTTGTGGCCACGGTTGGTTTAACGCTTTTGaatat                            >  2:582154/1‑145 (MQ=255)
                                                                                                      aCGGTTCAATGATTTATCATCGCCATGCAGTAAATAATCATGTGCCTGATAAAAATTAGTTAATAACGCACCACGATGCGGCAAAATTTTCTGGAGCGTCTCCTGCATTCGTTGTGGCCACGGTTGGTTTAACGCTTTTGATCAActgt                        >  2:617719/1‑147 (MQ=255)
                                                                                                      aCGGTTCAATGATTTATCATCGCCATGCAGTAAAT‑ATCATGTGCCTGATAAAAATTAGTTAATAACGCACCACGATGCGGCAAAATTTTCTGGAGCGTCTCCTGCATTCGTTGTGGCCACGGTTGGTTTAACGCTTTTGATAAACtctc                       <  1:724908/149‑1 (MQ=255)
                                                                                                      aCGGTTCAATGATTTATCATCGCCATGCAGTAAAT‑ATCATGTGCCTGATAAAAATTAGTTAATAACGCACCACGATGCGGCAAAATTTTCTGGAGCGTCTCCTGCATTCGTTGTGGCCACGGTTGGTTTAACGCTTTTGATAAACtctc                       <  1:650115/149‑1 (MQ=255)
                                                                                                       cGGTTCAATGATTTATCATCGCCATGCAGTAAATAATCATGTGCCTGATAAAAATTAGTTAATAACGCACCACGATGCGGCAAAATTTTCTGGAGCGTCTCCTGCATTCGTTGTGGCCACGGTTGGTTTAACGCTTTTGATAAACtctc                       <  2:75078/149‑1 (MQ=255)
                                                                                                          ttCAATGATTTATCATCGCCATGCAGTAAATAATCATATGCCTGAGAAAAATTAGTTAATAACGCACCACGATGCGGCAAAATTTTCTGGAGCGTCTCCTGCATTCGTTGTGGCCACGGTTGGTTTCACGCTTTTGATAAACTCTCCAg                    >  1:497357/1‑149 (MQ=255)
                                                                                                                   ttATCATCGCCATGCAGTAAATAATCATGTGCCTGATAAAAATTAGTTAATAACGCACCACGATGCGGCAAAATTTTCTGGAGCGTCTCCTGCATTCGTTGTGGCCACGGTTGGTTTAACGCTTTTGATAAACTCTCCAGTAAATCAtt           <  1:784502/149‑1 (MQ=255)
                                                                                                                      tcatcGCCATGCAGTAAATAATCATGTGCCTGATAAAAATTAGTTAATAACGCACCACGATGCGGCAAAATTTTCTGGAGCGTCTCCTGCATTCGTTGTGGCCACGGTTGGTTTAACGCTTTTGATAAACTCTCCAGTAAATCATTTTg        >  2:121981/1‑149 (MQ=255)
                                                                                                                        atcGCCATGCAGTAAATAATCATGTGCCTGATAAAAATTAGTTAATAACGCACCACGATGCGGCAAAATTTTCTGGAGCGTCTCCTGCATTCGTTGTGGCCACGGTTGGTTTAACGCTTTTGATAAACTCTCCAGTAAATCATTTTGaa      <  2:455087/149‑1 (MQ=255)
                                                                                                                            ccATGCAGTAAATAGTCATGTGCCTGATAAAAATTAGTTAATAACGCACCACGATGCGGCAAAATTTTCTGGAGCGTCTCCTGCATTAGCTGTGGCGACGGTTGGTTTAACGCTTTTGATAAACTCTCCAGTAAATCATTTTGAATCGc  >  1:768104/1‑149 (MQ=255)
                                                                                                                            ccATGCAGTAAATAATCATGTGCCTGATACAAATTAGTTAATAACGCACCACGATGCGGCAAAATTTGCTGGAGCGTCTCCTGCATTCGTTGTGGCCACGGTTGGTTTAACGCTTTTGATAAACTCTCCAGTAAATCATTTTGAATCGc  >  1:475377/1‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                       
AGGATGTTGAGAATGGCGCACGATATCAACTAATGCTTTTTCTGCTCTCGCGTAGGTAAATTCTGGGGATGATTGAACAATCTCACCTAATAATTCACTGGCACGGTTCAATGATTTATCATCGCCATGCAGTAAATAATCATGTGCCTGATAAAAATTAGTTAATAACGCACCACGATGCGGCAAAATTTTCTGGAGCGTCTCCTGCATTCGTTGTGGCCACGGTTGGTTTAACGCTTTTGATAAACTCTCCAGTAAATCATTTTGAATCGC  >  NC_000913/4360763‑4361035

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑