Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,984,411	C→T	*394* (TAG→TAA)	yqeF ←	short chain acyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,984,411	0	C	T	100.0%	165.0 / NA	52	*394* (TAG→TAA)	yqeF	short chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/25);  total (27/25)

TGTCGATCACGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAG  >  NC_000913/2984282‑2984538                                                                                                                                  |                                                                                                                                tgtCGATCACGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCa                                                                                                                        >  2:1222918/1‑139 (MQ=255)                  tCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGc                                                                                                       <  1:1194433/139‑1 (MQ=255)                            ggAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTgg                                                                                             >  2:537336/1‑139 (MQ=255)                                   tGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccgc                                                                                      >  2:828755/1‑139 (MQ=255)                                       ggAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTc                                                                                                                         >  1:1394601/1‑100 (MQ=255)                                       ggAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTc                                                                                                                         <  2:1394601/100‑1 (MQ=255)                                       ggAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTc                                                                                                            <  1:1482035/113‑1 (MQ=255)                                       ggAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTc                                                                                                            >  2:1482035/1‑113 (MQ=255)                                           cccAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACaa                                                                              <  1:537336/139‑1 (MQ=255)                                              aGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGc                                                                           <  1:828755/139‑1 (MQ=255)                                               gAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGcc                                                                                           <  1:763855/122‑1 (MQ=255)                                               gAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGcc                                                                                           >  2:763855/1‑122 (MQ=255)                                                            gACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTGAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGt                                                                                                             >  2:847749/1‑91 (MQ=255)                                                            gACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGt                                                                                                             <  1:847749/91‑1 (MQ=255)                                                              cTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGcc                                                                                           <  2:1012879/107‑1 (MQ=255)                                                              cTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGcc                                                                                           >  1:1012879/1‑107 (MQ=255)                                                              cTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATACAGAGTTATTCg                                                                                                                           <  2:1193827/75‑1 (MQ=255)                                                              cTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATACAGAGTTATTCg                                                                                                                           >  1:1193827/1‑75 (MQ=255)                                                               ttGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGcc                                                                                           >  2:1292788/1‑106 (MQ=255)                                                               ttGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGcc                                                                                           <  1:1292788/106‑1 (MQ=255)                                                                tGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCgg                                                         >  2:705264/1‑139 (MQ=255)                                                                                       aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGcc                                                                                                      >  2:103516/1‑71 (MQ=255)                                                                                       aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGcc                                                                                                      <  1:103516/71‑1 (MQ=255)                                                                                       aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                  >  2:616727/1‑139 (MQ=255)                                                                                         tAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAg                                                                            <  1:1137530/95‑1 (MQ=255)                                                                                         tAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAg                                                                            >  2:1137530/1‑95 (MQ=255)                                                                                         tAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGt                                                <  1:1485343/123‑1 (MQ=255)                                                                                         tAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGt                                                >  2:1485343/1‑123 (MQ=255)                                                                                           aTGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGtaa                                                                                                                        >  2:927626/1‑47 (MQ=255)                                                                                           aTGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCa                                                                                                                        <  1:927626/49‑1 (MQ=255)                                                                                            tGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCa                             >  2:1122358/1‑139 (MQ=255)                                                                                                 gCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGcc                                                                                                      >  2:874617/1‑61 (MQ=255)                                                                                                 gCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGcc                                                                                                      <  1:874617/61‑1 (MQ=255)                                                                                                 gCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGaa                        >  2:183786/1‑139 (MQ=255)                                                                                                   gagacagtATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAc                      <  1:46715/132‑1 (MQ=255)                                                                                                   aGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGcc                                                                                           >  1:764379/1‑70 (MQ=255)                                                                                                   aGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGcc                                                                                           <  2:764379/70‑1 (MQ=255)                                                                                                   aGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAc                      <  1:705264/139‑1 (MQ=255)                                                                                                     atatTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCaga                           <  2:794699/132‑1 (MQ=255)                                                                                                     atatTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCaga                           >  1:794699/1‑132 (MQ=255)                                                                                                          tATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGc                                                                      <  1:731682/84‑1 (MQ=255)                                                                                                          tATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGc                                                                      >  2:731682/1‑84 (MQ=255)                                                                                                                tAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCaa         >  1:857483/1‑139 (MQ=255)                                                                                                                  aggaggAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGcccc                                                                                         <  1:1370377/57‑1 (MQ=255)                                                                                                                  aggaggAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGcccc                                                                                         >  2:1370377/1‑57 (MQ=255)                                                                                                                     aggaAATGCAGAGTTATTCGTCACGTTCAATGGTCAATACAACACCCTGGACCCCGCCGATAAAAAACATTTCAAGTCGTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGcctgtctc    >  2:46715/1‑132 (MQ=255)                                                                                                                       gaAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAg  <  1:1222918/139‑1 (MQ=255)                                                                                                                         aaTGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCa                                                                    <  1:1292403/71‑1 (MQ=255)                                                                                                                         aaTGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCa                                                                    >  2:1292403/1‑71 (MQ=255)                                                                                                                          aTGCAGAGTTATTCGTCACGTTCAATGTTAAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAtt               <  2:257059/123‑1 (MQ=255)                                                                                                                          aTGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAtt               >  1:257059/1‑123 (MQ=255)                                                                                                                               gagTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCttt                                                             <  1:381941/72‑1 (MQ=255)                                                                                                                                  |                                                                                                                                TGTCGATCACGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAG  >  NC_000913/2984282‑2984538

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