Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,433,146	T→C	T139A (ACA→GCA)	smf ←	protein Smf

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,433,146	0	T	C	100.0%	133.9 / NA	39	T139A (ACA→GCA)	smf	protein Smf
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/17);  total (22/17)

AGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTACACCATCGATTCCACGCGCCAGTCCACTCGTAATTGTCACTCCACGCGTCGCCAGAGTTTCGCAAAATAATCGTCCCCATCGCTCGCCATACCATGAATGCGCCCGACTCCCCACTACGGCAAGCTGAAATGAATGCAGCGCGTGCAGTTCTCCT  >  NC_000913/3433016‑3433264                                                                                                                                   |                                                                                                                       aGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTACACCATCGATTCCACGCGCCAGTCCACTCGTAATTGCCa                                                                                                                      >  1:64998/1‑133 (MQ=255)        cATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTACACCATCGATTCCACGCGCCAGTCCACTCGTAATTGCCACTCCAcg                                                                                                               <  1:494348/133‑1 (MQ=255)                          tGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTACACCATCGATTCCACGCGCCAGTCCACTCGTAATTGCCACTCCACGCGTCGCCAGAGTTTCGCa                                                                                             >  2:411378/1‑133 (MQ=255)                             tATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTACACCATCGATTCCACGCGCCAGTCCACTCGTAATTGCCACTCCACGCGTCGCCAGAGTTTCGCaaaa                                                                                          >  2:17964/1‑133 (MQ=255)                             tATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTACACCATCGATTCCACGCGCCAGTCCACTCGTAATTGCCACTCCACGCGTCGCCAGAGTTTCGCaaaa                                                                                          >  1:494463/1‑133 (MQ=255)                               ttAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTACACCATCGATTCCACGCGCCAGTCCACTCGTAATTGCCACTCCATGCGTCGCCAGAGTTTCGCAAaata                                                                                        >  2:393662/1‑133 (MQ=255)                               ttAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTACACCATCGATTCCACGCGCCAGTCCACTCGTAATTGCCACTCCACGCGTCGCCAGAGTTTCGCAAaata                                                                                        >  2:285721/1‑133 (MQ=255)                                  aGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTACACCATCGATTCCACGCGCCAGTCCACTCGTAATTGCCACTCCACGCGTCGCCAGAGTTTCGCAAAATAATc                                                                                     >  1:504968/1‑133 (MQ=255)                                       aTTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTACACCATCGATTCCACGCGCCAGTCCACTCGTAATTGCCACTCCACGCGTCGCCAGAGTTTCGCAAAATAATCGTccc                                                                                >  2:45074/1‑133 (MQ=255)                                        ttCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTACACCATCGATTCCACGCGCCAGTCCACTCGTAATTGCCACTCCACGCGTCGCCAGAGTTTCGCAAAATAATCGTcccc                                                                               >  2:392057/1‑133 (MQ=255)                                             cAATACAGTAATGCTGACGCCACTTACCTGTAAGGCTGCTTTATGCGCTACACCATCGATTCCACGCGCCAGTCCACTCGTAATTGCCACTCCACGCGTCGCCAGAGTTTCGCAAAATAATCGTCCCCAtcgc                                                                          <  1:185070/133‑1 (MQ=255)                                               aTACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTACACCATCGATTCCACGCGCCAGTCCACTCGTAATTGCCACTCCACGCGTCGCCAGAGTTTCGCAAAATAATCGTCCCCAtcgctc                                                                        <  2:504968/133‑1 (MQ=255)                                                   aGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTACACCATCGATTCCACGCGCCAGTCCACTCGTAATTGCCACTCCACGCGTCGCCAGAGTTTCGCAAAATAATCGTCCCCATCGCTCGCCa                                                                    >  1:384113/1‑133 (MQ=255)                                                    ggaaTGCTGTCGCCATTTACCTGTAAAGCTGCTTTATGCGCGACACCATCGATTCCACGCGCCAGTCCACTCGTAATTGCCACTCCACGCGTCGCCAGAGTTTCGCAAAATAATCGTCCCCATCGCTCGCCAt                                                                   <  2:249357/131‑1 (MQ=255)                                                     cAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTACACCATCGATTCCACGCGCCAGTCCACTCGTAATTGCCACTCCACGCGTGGCCAGAGTTTCGTAAAGTAATCGTCTCCATTGCTCGCCATa                                                                  >  2:148388/1‑133 (MQ=255)                                                     cAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTACACCATCGATTCCACGCGCCAGTCCACTCGTAATTGCCACTCCACGCGTCGCCAGAGTTTCGCAAAATAATCGTCCCCATCGCTCGCCATa                                                                  >  2:42884/1‑133 (MQ=255)                                                     cAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTACACCATCGATTCCACGCGCCAGTCCACTCGTAATTGCCACTCCACGCGTCGCCAGAGTTTCGCAAAATAATCGTCCCCATCGCTCGCCATa                                                                  >  1:79505/1‑133 (MQ=255)                                                     cAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTACACCATCGATTCCACGCGCCAGTCCACTCGTAATTGCCACTCCACGCGTCGCCAGAGTTTCGCAAAATAATCGTCCCCATCGCTCGCCATa                                                                  >  2:107879/1‑133 (MQ=255)                                                     cAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTACACCATCGATTCCACGCGCCAGTCCACTCGTAATTGCCACTCCACGCGTCGCCAGAGCTTCGCAAAATAATCGTCCCCATCGCTCGCCATa                                                                  >  2:492439/1‑133 (MQ=255)                                                             aCGCCATTTACCTGTAAGGCTGCTTTATGCGCTACACCATCGATTCCACGCGCCAGTCCACTCGTAATTGCCACTCCACGCGTCGCCAGAGTTTCGCAAAATAATCGTCCCCATCGCTCGCCATACCATGAAt                                                          >  2:1127/1‑133 (MQ=255)                                                                ccATTTACCTGTAAGGCTGCTTTATGCGCTACACCATCGATTCCACGCGCCAGTCCACTCGTAATTGCCACTCCATGCGTCGCCAGAGTTTCGCAAAATAATCGTCCCCATCGCTCGCCATACCATGAATgcg                                                       <  1:393662/133‑1 (MQ=255)                                                                  aTTTACCTATAAGGCTGCATTAAGCGCTACGCCATCGATTCCACGCGCCAGTCCACTCGTAATTGCCACTCCACGCGTCGCCAGAGTTTCGCAAAATAATCGTCCCCATCGCTCGCCATACCATGAATGCGcc                                                     <  1:217525/133‑1 (MQ=255)                                                                     tACCTGTAAGGCTGCTTTATGCGCTACACCATCGATTCCACGCGCCAGTCCACTCGTAATTGCCACTCCACGCGTCGCCAGAGTTTCGCAAAATAATCGTCCCCATCGCTCGCCATACCATGAATGCGCCCGa                                                  <  1:285721/133‑1 (MQ=255)                                                                     tACCTGTAAGGCTGCGTTTTGCTCTACACCATCGATTCCACGCGCCAGTCCACTCGTAATTGCCACTCCACGCGTCGCCAGAGTTTCGCAAAATAATCGTCCCCATCGCTCGCCATACCATGAATGCGCCCGa                                                  <  1:181866/133‑1 (MQ=255)                                                                        cTGTAAGGCTGCTTTATGCGCTACACCATCGATTCCACGCGCCAGTCCTCTCGTAATTGCCACTCCACGCGTCGCCAGAGTTTCGCAAAATAATCGTCCCCATCGCTCGCCATACCATGAATGCGCCCGACTc                                               <  1:493971/133‑1 (MQ=255)                                                                           tAAGGCTGCTTTATGCGCTACACCATCGATTCCACGCGCCAGTCCACTCGTAATTGCCACTCCACGCGTCGCCAGAGTTTCGCAAAATAATCGTCCCCATCGCTCGCCATACCATGAATGCGCCCGACTcccc                                            >  1:314546/1‑133 (MQ=255)                                                                                  gctTTATGCGCTACACCATCGATTCCACGCGCCAGTCCACTCGTAATTGCCACTCCACGCGTCGCCAGAGTTTCGCAAAATAATCGTCCCCATCGCTCGCCATACCATGAATGCGCCCGACTCCCCACTACgg                                     >  2:96349/1‑133 (MQ=255)                                                                                      tATGCGCTACACCATCGATTCCACGCGCCAGTCCACTCGTAATTGCCACTCCACGCGTCGCCAGAGTTTCGCAAAATAATCGTCCCCATCGCTCGCCATACCATGAATGCGCCCGACTCCCCACTACGGCAAg                                 <  1:107879/133‑1 (MQ=255)                                                                                             tACACCATCGATTCCACGCGCCAGTCCACTCGTAATTGCCACTCCACGCGTCGCCAGAGTTTCGCAAAATAATCGTCCCCATCGCTCGCCATACCATGAATGCGCCCGACTCCCCACTACGGCAAGCTGAaat                          >  2:409299/1‑133 (MQ=255)                                                                                                    tCGATTCCACGCGCCAGTCCACTCGTAATTGCCACTCCACGCGTCGCCAGAGTTTCGCAAAATAATCGTCCCCATCGCTCGCCATACCATGAATGCGCCCGACTCCCCACTACGGCAAGCTGAAATGAATGCa                   >  2:308341/1‑133 (MQ=255)                                                                                                      gATTCCACGCGCCAGTCCACTCGTAATTGCCACTCCACGCGTCGCCAGAGTTTCGCAAAATAATCGTCCCCATCGCTCGCCATACCATGAATGCGCCCGACTCCCCACTACGGCAAGCTGAAATGAATGCAgc                 <  1:308341/133‑1 (MQ=255)                                                                                                         tCCACGCGCCAGTCCACTCGTAATTGCCACTCCACGCGTCGTCAGAGTTTCGCAAAATAATCGTCCCCATCGCTCGCCATACCATGAATGCGCCCGACTCACCAGTACAGCAAGCTGGAATGGATGCAGcgcg              >  2:145034/1‑133 (MQ=255)                                                                                                         tCCACGCGCCAGTCCACTCGTAATTGCCACTCCACGCGTCGCCAGAGTTTCGCAAAATAATCGTCCCCATCGCTCGCCATACCATGAATGCGCCCGACTCCCCACTACGGCAAGCTGAAATGAATGCAGcgcg              <  1:22862/133‑1 (MQ=255)                                                                                                           cACGCGCCAGTCCACTCGTAATTGCCACTCCACGCGTCGCCAGAGTTTCGCAAAATAATCGTCCCCATCGCTCGCCATACCATGAATGCGCCCGACTCCCCACTACGGCAAGCTGAAATGAATGCAGCGCGTg            <  1:409299/133‑1 (MQ=255)                                                                                                            aCGCGCCAGTCCACTCGTAATTGCCACTCCACGCGTCGCCAGAGTTTCGCAAAATAATCGTCCCCATCGCTCGCCATACCATGAATGCGCCCGACTCCCCACTACGGCAAGCTGAAATGAATGCAGCGCGTGc           >  1:417480/1‑133 (MQ=255)                                                                                                               cgcCAGTCCACTCGTAATTGCCACTCCACGCGTCGCCAGAGTTTCGCAAAATAATCGTCCCCATCGCTCGCCATACCATGAATGCGCCCGACTCCCCACTACGGCAAGCTGAAATGAATGCAGCGCGTGCAGt        <  1:105839/133‑1 (MQ=255)                                                                                                                 ccAGTCCACTCGTAATTGCCACTCCACGCGTCGCCAGAGTTTCGCAAAATAATCGTCCCCATCGCTCGCCATACCATGAATGCGCCCGACTCCCCACTACGGCAAGCTGAAATGAATGCAGCGCGTGCAGTtc      <  2:64998/133‑1 (MQ=255)                                                                                                                   aGTCCACTCGTAATTGCCACTCCACGCGTCGCCAGAGTTTCGCAAAATAATCGTCCCCATCGCTCGCCATACCATGAATGCGCCCGACTCCCCACTACGGCAAGCTGAAATGAATGCAGCGCGTGCAGTtctc    <  1:42884/133‑1 (MQ=255)                                                                                                                     tCCACTCGTAATTGCCACACCACGCGTCGTTAGAGTTTCGCAAAATAATCGTCCCCATCGCTCGCCATACCATGAATGCGCCCGACTCCCCACTACGGCAAGCTGAAATGAATGCAGCGCGTGCAGTTCTCCt  <  1:145034/133‑1 (MQ=255)                                                                                                                                   |                                                                                                                       AGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATGCGCTACACCATCGATTCCACGCGCCAGTCCACTCGTAATTGTCACTCCACGCGTCGCCAGAGTTTCGCAAAATAATCGTCCCCATCGCTCGCCATACCATGAATGCGCCCGACTCCCCACTACGGCAAGCTGAAATGAATGCAGCGCGTGCAGTTCTCCT  >  NC_000913/3433016‑3433264

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