Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,432,994	(C)6→7	coding (567/1125 nt)	smf ←	DNA recombination‑mediator A family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,432,988	1	.	C	100.0%	162.0 / NA	46	G191G (GGC→GGG)	smf	DNA recombination‑mediator A family protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (26/20);  total (26/20)

TGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTA  >  NC_000913/3432850‑3433109                                                                                                                                            |                                                                                                                          tGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAgcg                                                                                                                            <  1:889262/139‑1 (MQ=255) tGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAgcg                                                                                                                            >  1:28794/1‑139 (MQ=255)      tCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAGATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑ccccc                                                                                                                      <  2:28794/139‑1 (MQ=255)         ccAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑ccccccct                                                                                                                   >  1:952499/1‑137 (MQ=255)                aCCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAAACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAg                                                                                                             >  2:801747/1‑139 (MQ=255)                            aGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGc                                                                                                 >  2:673126/1‑139 (MQ=255)                            aGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑ccccccc                                                                                                                    >  2:553038/1‑118 (MQ=255)                            aGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑ccccccc                                                                                                                    <  1:553038/119‑2 (MQ=255)                                              tACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcggcg                                                                               >  2:526116/1‑139 (MQ=255)                                                   cTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGAt                                                                          >  1:322226/1‑139 (MQ=255)                                                   cTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGAt                                                                          >  2:703021/1‑139 (MQ=255)                                                   cTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑ccccc                                                                                                                      >  1:845566/1‑94 (MQ=255)                                                   cTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑ccccc                                                                                                                      <  2:845566/94‑1 (MQ=255)                                                     ttACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGa                                                                        >  1:607419/1‑139 (MQ=255)                                                     ttACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGa                                                                        <  1:801747/139‑1 (MQ=255)                                                      tACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGaa                                                                       <  2:607419/139‑1 (MQ=255)                                                      tACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGaa                                                                       >  2:933625/1‑139 (MQ=255)                                                       aCTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                      <  2:368073/139‑1 (MQ=255)                                                        cTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATg                                                                     >  1:387095/1‑139 (MQ=255)                                                             aCCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCTGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                      <  2:559196/133‑1 (MQ=255)                                                             aCCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCTGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                      >  1:559196/1‑133 (MQ=255)                                                             aCCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTAtt                                                                >  2:60549/1‑139 (MQ=255)                                                             aCCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTAtt                                                                >  2:291546/1‑139 (MQ=255)                                                                  tGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                      >  1:882702/1‑128 (MQ=255)                                                                  tGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAAt                                                                      <  2:882702/128‑1 (MQ=255)                                                                  tGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTc                                                           >  1:494169/1‑139 (MQ=255)                                                                  tGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTc                                                           >  2:351534/1‑139 (MQ=255)                                                                         gCGATTTTTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAGTTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTccc                                                    <  2:952499/139‑1 (MQ=255)                                                                         gCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTccc                                                    <  1:933625/139‑1 (MQ=255)                                                                            aTTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCaa                                                 >  2:666328/1‑139 (MQ=255)                                                                                   cGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAACGGTATTAAGTCCATTCCCCAATACAGCa                                          <  1:703021/139‑1 (MQ=255)                                                                                          aaTTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGa                                   <  1:351534/139‑1 (MQ=255)                                                                                            ttGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACg                                 <  1:60549/139‑1 (MQ=255)                                                                                                       gggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGa                                   >  1:45585/1‑126 (MQ=255)                                                                                                       gggTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGa                                   <  2:45585/126‑1 (MQ=255)                                                                                                             aaCATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGgc                <  2:494169/139‑1 (MQ=255)                                                                                                                  cGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGa                                                                                    >  1:365816/1‑66 (MQ=255)                                                                                                                  cGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGa                                                                                    <  2:365816/66‑1 (MQ=255)                                                                                                                  cGAGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCtt           >  2:930744/1‑139 (MQ=255)                                                                                                                    aGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCa                                          >  1:142567/1‑106 (MQ=255)                                                                                                                    aGGGGAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCa                                          <  2:142567/106‑1 (MQ=255)                                                                                                                       ggAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCAt                                                                                      >  2:764319/1‑59 (MQ=255)                                                                                                                       ggAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCAt                                                                                      <  1:764319/59‑1 (MQ=255)                                                                                                                       ggAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACg                                 >  2:241422/1‑112 (MQ=255)                                                                                                                       ggAAATTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACg                                 <  1:241422/112‑1 (MQ=255)                                                                                                                           aTTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAg                                                                                                  <  2:246204/43‑1 (MQ=255)                                                                                                                           aTTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAg                                                                                                  >  1:246204/1‑43 (MQ=255)                                                                                                                           aTTCCGAGACGAGAGCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTa  >  2:689936/1‑139 (MQ=255)                                                                                                                                gagacgagaGCGCCCCCCCTGTTCAAGCAGACTGGCAGcc                                                                                                <  2:1001198/40‑1 (MQ=255)                                                                                                                                gagacgagaGCGCCCCCCCTGTTCAAGCAGACTGGCAGcc                                                                                                >  1:1001198/1‑40 (MQ=255)                                                                                                                                      agagCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAg                                                                                              <  2:768725/36‑1 (MQ=255)                                                                                                                                      agagCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAg                                                                                              >  1:768725/1‑36 (MQ=255)                                                                                                                                      agagCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCttt          >  2:20831/1‑120 (MQ=255)                                                                                                                                      agagCGCCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCttt          <  1:20831/120‑1 (MQ=255)                                                                                                                                            |                                                                                                                          TGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑CCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTA  >  NC_000913/3432850‑3433109

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