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%	195.2 / NA	62	*394* (TAG→TAA)	yqeF	short chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (32/30);  total (32/30)

TGTCGATCACGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTA  >  NC_000913/2984282‑2984544                                                                                                                                  |                                                                                                                                      tgtCGATCACGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCa                                                                                                                              >  1:234636/1‑139 (MQ=255)   tCGATCACGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACg                                                                                                                            >  1:473478/1‑139 (MQ=255)              tatGTAATGAGAATGGAAAAATGAGGGAGCCGAGAAAATGTATGTCGACTTGCCAGCCATAATTATTCAACAGATGAATGACGACAGATATTTAACGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCa                                                                                                                 <  2:680233/139‑1 (MQ=255)                 gTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATg                                                                                                              >  1:764924/1‑139 (MQ=255)                     tGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTc                                                                                                                               >  1:1320888/1‑118 (MQ=255)                     tGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTc                                                                                                                               <  2:1320888/118‑1 (MQ=255)                     tGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCcac                                                                                                          >  2:269188/1‑139 (MQ=255)                     tGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCcaa                                                                                                          >  1:204275/1‑138 (MQ=255)                       agaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCca                                                                                                           <  2:1308367/136‑1 (MQ=255)                       agaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCca                                                                                                           >  1:1308367/1‑136 (MQ=255)                           tGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTg                                                                                                    >  2:869685/1‑139 (MQ=255)                            ggAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTCTTCGTCACGTTCAATGGTCAATGCCACACCCTgg                                                                                                   >  1:863707/1‑139 (MQ=255)                            ggAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTgg                                                                                                   >  1:1137309/1‑139 (MQ=255)                              aaaaaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGcc                                                                                                 <  2:234636/139‑1 (MQ=255)                              aaaaaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGcc                                                                                                 <  1:269188/139‑1 (MQ=255)                              aaaaaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGcc                                                                                                 <  2:1137309/139‑1 (MQ=255)                                        gAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATa                                                                                       >  1:663758/1‑139 (MQ=255)                                                     gtatgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGcc                                                                                                 >  2:461054/1‑116 (MQ=255)                                                     gtatgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGcc                                                                                                 <  1:461054/116‑1 (MQ=255)                                                       atgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGt                                                                        <  1:869685/139‑1 (MQ=255)                                                       atgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGACGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGt                                                                        <  2:863707/139‑1 (MQ=255)                                                          tCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCt                                                                     >  1:195209/1‑139 (MQ=255)                                                              cTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGcc                                                                                                            >  2:1056458/1‑96 (MQ=255)                                                              cTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGcc                                                                                                            <  1:1056458/96‑1 (MQ=255)                                                               ttGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCg                                                                <  2:195209/139‑1 (MQ=255)                                                                        atcatTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACg                                                       <  2:663758/139‑1 (MQ=255)                                                                          catTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACATGCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCGGTCCTTTGCGGGCATTACGtt                                                     >  2:785378/1‑139 (MQ=255)                                                                                 cAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTc                                                                                                                  <  1:1147665/71‑1 (MQ=255)                                                                                 cAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTc                                                                                                                  >  2:1147665/1‑71 (MQ=255)                                                                                    cAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc                                           <  2:764924/139‑1 (MQ=255)                                                                                    cAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc                                           >  1:1204407/1‑139 (MQ=255)                                                                                    cAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc                                           >  1:1313529/1‑139 (MQ=255)                                                                                     aGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCa                                          <  2:337892/139‑1 (MQ=255)                                                                                       aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                        <  2:473478/139‑1 (MQ=255)                                                                                       aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                        <  2:1204407/139‑1 (MQ=255)                                                                                          aaTGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGc                                                                            <  1:1229320/100‑1 (MQ=255)                                                                                          aaTGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGc                                                                            >  2:1229320/1‑100 (MQ=255)                                                                                              aCGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCaga                                 >  1:444261/1‑139 (MQ=255)                                                                                                 gCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGt                                                                                                                                <  1:1050683/41‑1 (MQ=255)                                                                                                 gCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGt                                                                                                                                >  2:1050683/1‑41 (MQ=255)                                                                                                 gCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCa                                                                                                                        >  1:373185/1‑49 (MQ=255)                                                                                                 gCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCa                                                                                                                        <  2:373185/49‑1 (MQ=255)                                                                                                 gCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGcc                                                                                                            <  1:474698/61‑1 (MQ=255)                                                                                                 gCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGcc                                                                                                            >  2:474698/1‑61 (MQ=255)                                                                                                 gCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTg                                                                  <  1:56768/103‑1 (MQ=255)                                                                                                 gCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTg                                                                  >  2:56768/1‑103 (MQ=255)                                                                                                       atTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCg                                                                <  1:704644/99‑1 (MQ=255)                                                                                                       atTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCg                                                                >  2:704644/1‑99 (MQ=255)                                                                                                               gTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCaga                                 >  1:575637/1‑122 (MQ=255)                                                                                                               gTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCaga                                 <  2:575637/122‑1 (MQ=255)                                                                                                                tAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCaa               >  1:1328932/1‑139 (MQ=255)                                                                                                                      ggaAATGCAGAGTTATTCGTCACGTTCAATGGTCa                                                                                                                 <  1:620855/35‑1 (MQ=255)                                                                                                                      ggaAATGCAGAGTTATTCGTCACGTTCAATGGTCa                                                                                                                 >  2:620855/1‑35 (MQ=255)                                                                                                                        aaaTGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgc       >  1:428845/1‑139 (MQ=255)                                                                                                                          aTGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAt                                                           <  1:752010/85‑1 (MQ=255)                                                                                                                          aTGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAt                                                           >  2:752010/1‑85 (MQ=255)                                                                                                                             cAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTa  <  2:1328932/139‑1 (MQ=255)                                                                                                                              agagTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCaga                                 >  2:1134716/1‑107 (MQ=255)                                                                                                                              agagTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCaga                                 <  1:1134716/107‑1 (MQ=255)                                                                                                                               gagTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAc                            <  2:348150/111‑1 (MQ=255)                                                                                                                               gagTTATTCGTCACGTTAAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAGTACGTTTCACCATTTCATGAACCAGAGAAAc                            <  2:234057/111‑1 (MQ=255)                                                                                                                               gagTTATTCGTCACGATCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAc                            >  1:234057/1‑111 (MQ=255)                                                                                                                                  |                                                                                                                                      TGTCGATCACGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTA  >  NC_000913/2984282‑2984544

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