Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,365,034 T→A 8.9% Q3H (CAA→CAT) 
N108I (AAC→ATC) 		 dsbD ←
cutA ←		 protein disulfide oxidoreductase ‑ DsbDreduced
copper binding protein CutA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,365,0340TA 8.9% 119.5 / 2.0 45Q3H (CAA→CAT) 
N108I (AAC→ATC) 		dsbD
cutA		protein disulfide oxidoreductase ‑ DsbDreduced
copper binding protein CutA
Reads supporting (aligned to +/- strand):  ref base T (22/19);  new base A (2/2);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.58e-01

ATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTC  >  NC_000913/4364896‑4365176
                                                                                                                                          |                                                                                                                                              
aTTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAg                                                                                                                                            >  1:369459/1‑143 (MQ=255)
aTTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGATGAg                                                                                                                                            >  1:349855/1‑143 (MQ=255)
aTTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGATGAg                                                                                                                                            >  1:205340/1‑143 (MQ=255)
    aGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCAt                                                                                                                                        <  1:639407/143‑1 (MQ=255)
      gTCATGTTGGTTTTGCTGAAAAGCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATga                                                                                                                                      <  1:416720/143‑1 (MQ=255)
        cATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATgaga                                                                                                                                    <  1:575998/143‑1 (MQ=255)
                tgttGCTGACAATCGATATCACAGGCTTGACCCGCGGGGACAAACTGGGAACATCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCt                                                                                                                            <  2:419303/141‑1 (MQ=255)
                      tGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGtct                                                                                                                      >  1:622551/1‑143 (MQ=255)
                      tGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTACCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAGAGATGCGTTGAGCAGAGAGAAGTTGTCtgtggcc                                                                                                                      >  2:43595/1‑140 (MQ=255)
                        aaaaTCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTcc                                                                                                                    >  2:329664/1‑143 (MQ=255)
                            tCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCgtgt                                                                                                                >  1:448311/1‑143 (MQ=255)
                            tCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGCCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTCAAGATGCGTTGAGCCATGAGAGGTAATCTGTGACTCCgggt                                                                                                                >  2:12237/1‑143 (MQ=255)
                                  gCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACa                                                                                                          >  1:535250/1‑143 (MQ=255)
                                          ttGATCCGCGGGGACAAATTGTGACCGTGCCGGCGAGTCGAATAATCCGGCAAAGAAGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTGGAGCCACGAGAGGTACTCTGTGTCTCCGTGTGTGACAGGTagagc                                                                                                  >  2:659353/1‑139 (MQ=255)
                                             aTCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGATAAAACCAg                                                                                               >  2:43626/1‑143 (MQ=255)
                                              tgcgcgAGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCACGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGa                                                                                              <  2:643899/141‑1 (MQ=255)
                                                 gcgGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGAACAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACAAGCAga                                                                                           >  2:436764/1‑142 (MQ=255)
                                                    gggTCAAAGTGTGAATGTCCCGGCGAGTCTAATTATCCGGCAAGAACGGTAGTGCTGCAAGGTAGCAGGATCAGCGTAAAGGTGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTcc                                                                                        <  2:124960/143‑1 (MQ=255)
                                                     ggACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCg                                                                                       >  1:435159/1‑143 (MQ=255)
                                                      taCAAATTGTGAACGCCCCGGCGCGTCGAAGAACCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCgg                                                                                      <  2:38763/142‑1 (MQ=255)
                                                      gACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGCGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCgg                                                                                      <  1:149835/143‑1 (MQ=255)
                                                       aCAAACTGTGAGCGTCCCGGCGCGTCGAAGAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCggg                                                                                     <  1:220367/143‑1 (MQ=255)
                                                          aaTTGTGAACCTCCCGCCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGtt                                                                                  <  1:607176/143‑1 (MQ=255)
                                                              gtgAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGGAGTTCCGGGGtgcgac                                                                              >  2:564979/1‑138 (MQ=255)
                                                               tgAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGata                                                                             <  1:204225/143‑1 (MQ=255)
                                                                      cccGGCGCGGCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGgatga                                                                      <  1:565377/143‑1 (MQ=255)
                                                                                    taatCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGATGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCa                                                        <  2:205340/143‑1 (MQ=255)
                                                                                    gaatCCGGCAAAAACGGGAGTGCGGCAAAGTAGCAGGATCAGCGTAAAGATGCGATGAGCCATGAGAGTTAATCTTTGTGTTAGTGTAAAACAGGTAAAACCAGAAGTTCCGGGCTTTGATATGGATGATGAGACTTCAGGCa                                                        <  2:349855/142‑1 (MQ=255)
                                                                                       tgcGGCAAGAACGGAAGGGCGGCAAAGTGGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTGCGTGTCTAACAGGTAAAACCAGAAGTTCCGTGGTTTGATATGTATGATGACACTTCAGGCATTc                                                     <  2:389614/141‑1 (MQ=255)
                                                                                         cGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCca                                                   >  2:601093/1‑143 (MQ=255)
                                                                                           gCAAAAGCGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAGAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTGCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCcagc                                                 <  1:213470/143‑1 (MQ=255)
                                                                                               aaaCGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTg                                             >  2:259370/1‑143 (MQ=255)
                                                                                                aaCGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTGCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCCGTGc                                            >  2:8378/1‑143 (MQ=255)
                                                                                                      aGTGCTGCACAGTAGCAGGATCAGCGTAAAGATGCGTCGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTCGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCctgctg                                      <  1:348391/143‑1 (MQ=255)
                                                                                                           tgcAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTACAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGGgag                                 >  2:668511/1‑143 (MQ=255)
                                                                                                           tgcAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTgag                                 >  1:646832/1‑143 (MQ=255)
                                                                                                              aaaGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATa                              <  2:435159/143‑1 (MQ=255)
                                                                                                              aaaGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGGTTGACATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGGTa                              >  2:221127/1‑143 (MQ=255)
                                                                                                              aaaGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGGGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACCTCAGGCATTCCAGCAGTGCATGCTGGGGAGATa                              >  2:495257/1‑143 (MQ=255)
                                                                                                                      aGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGtc                      <  1:15201/143‑2 (MQ=255)
                                                                                                                       ggATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGttt                     <  2:415478/143‑1 (MQ=255)
                                                                                                                       ggAGCAGGGTAAAGGTGCGTTGAGCCATGAGAGGGAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGttt                     <  2:415476/143‑1 (MQ=255)
                                                                                                                                   aGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCt         >  1:41522/1‑143 (MQ=255)
                                                                                                                                   aGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAAACAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATCCGGTAGTCTTTAAActcatgt         >  2:328005/1‑141 (MQ=255)
                                                                                                                                    gATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGAATTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTg        >  1:361622/1‑143 (MQ=255)
                                                                                                                                          ttGAGCCATGCGAGGTAATCTGGGTCTGCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTc  <  2:403600/143‑1 (MQ=255)
                                                                                                                                          ttGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTc  >  1:105840/1‑143 (MQ=255)
                                                                                                                                          |                                                                                                                                              
ATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGCGTTGAGCCATGAGAGGTAATCTGTGTCTCCGTGTGTAACAGGTAAAACCAGAAGTTCCGGGGTTTGATATGGATGATGAGACTTCAGGCATTCCAGCAGTGCCTGCTGGTGAGATACGGTAGTTTTTAAAATCATCTGCACTTC  >  NC_000913/4364896‑4365176

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