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

GTGCAGTGTCGATCACGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAG  >  NC_000913/2984276‑2984538                                                                                                                                        |                                                                                                                                gTGCAGTGTCGATCACGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTAt                                                                                                                              >  1:1092229/1‑139 (MQ=255)            aTCACGACTATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTc                                                                                                                   >  2:1005961/1‑139 (MQ=255)            aTCACGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTc                                                                                                                   >  1:790941/1‑139 (MQ=255)            aTCACGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTc                                                                                                                   >  2:1189025/1‑139 (MQ=255)                  aaTATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGt                                                                                                             <  1:1189025/139‑1 (MQ=255)                  aaTATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGt                                                                                                             <  2:383560/139‑1 (MQ=255)                              gaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACAcc                                                                                                 <  1:460732/139‑1 (MQ=255)                                 tGGAAAAATGAGGGAGCCCAGAAAATGTGTGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGAGGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACCCCCTg                                                                                              >  1:952607/1‑139 (MQ=255)                                    aaaaaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGcc                                                                                           <  2:790941/139‑1 (MQ=255)                                       aaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCa                                                                                                                  <  2:1318877/113‑1 (MQ=255)                                       aaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCa                                                                                                                  >  1:1318877/1‑113 (MQ=255)                                             ggAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGAt                                                                                  >  1:584841/1‑139 (MQ=255)                                                 cccAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACaa                                                                              >  2:992909/1‑139 (MQ=255)                                                    aGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGc                                                                           <  2:1092229/139‑1 (MQ=255)                                                    aGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGc                                                                           >  2:417205/1‑139 (MQ=255)                                                                tCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGt                                                                                                             >  2:484098/1‑93 (MQ=255)                                                                tCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGt                                                                                                             <  1:484098/93‑1 (MQ=255)                                                                     ttGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTTGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCg                                                          <  2:1187083/139‑1 (MQ=255)                                                                       gccagccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCa                                                                                                                  >  2:1050090/1‑81 (MQ=255)                                                                       gccagccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCa                                                                                                                  <  1:1050090/81‑1 (MQ=255)                                                                        ccagccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGcc                                                                     <  1:895504/125‑1 (MQ=255)                                                                        ccagccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTAGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGcc                                                                     >  2:895504/1‑125 (MQ=255)                                                                                          cAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAtt                                                    >  2:347020/1‑124 (MQ=255)                                                                                          cAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAtt                                                    <  1:347020/124‑1 (MQ=255)                                                                                             aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCCCGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTg                                                            <  2:1358932/113‑1 (MQ=255)                                                                                             aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCa                                                                                                                  <  1:1495740/59‑1 (MQ=255)                                                                                             aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCa                                                                                                                  >  2:1495740/1‑59 (MQ=255)                                                                                             aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTg                                                            >  1:1358932/1‑113 (MQ=255)                                                                                             aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGtt                                               >  1:1077495/1‑126 (MQ=255)                                                                                             aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGtt                                               >  1:1118058/1‑126 (MQ=255)                                                                                             aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGtt                                               <  2:1077495/126‑1 (MQ=255)                                                                                             aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGtt                                               <  2:1118058/126‑1 (MQ=255)                                                                                             aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCa                                         <  2:1174930/132‑1 (MQ=255)                                                                                             aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCa                                         >  1:1174930/1‑132 (MQ=255)                                                                                             aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                  >  2:637494/1‑139 (MQ=255)                                                                                             aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                  >  2:1295278/1‑139 (MQ=255)                                                                                               tAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTACAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGAGGGCATTACGTTTCACAATTTCATGaa                                >  2:790358/1‑139 (MQ=255)                                                                                                    aCGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGt                                                                                                             <  1:460164/57‑1 (MQ=255)                                                                                                    aCGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGt                                                                                                             >  2:460164/1‑57 (MQ=255)                                                                                                    aCGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCaga                           >  2:803772/1‑139 (MQ=255)                                                                                                    aCGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCaga                           <  1:1031678/139‑1 (MQ=255)                                                                                                     cGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                  >  1:1315718/1‑131 (MQ=255)                                                                                                     cGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                  <  2:1315718/131‑1 (MQ=255)                                                                                                       gCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccgc                                                                                      <  1:544618/77‑1 (MQ=255)                                                                                                       gCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccgc                                                                                      >  2:544618/1‑77 (MQ=255)                                                                                                         aGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAc                      <  1:1427968/139‑1 (MQ=255)                                                                                                            tatTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAg                   >  2:1235473/1‑139 (MQ=255)                                                                                                             atTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAgg                  >  2:748902/1‑139 (MQ=255)                                                                                                              tttATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGa                 <  1:748902/139‑1 (MQ=255)                                                                                                                        aggaggAAATGCAGAGTTATTCGTCACGTTCAATGGt                                                                                                             >  2:1429440/1‑37 (MQ=255)                                                                                                                        aggaggAAATGCAGAGTTATTCGTCACGTTCAATGGt                                                                                                             <  1:1429440/37‑1 (MQ=255)                                                                                                                           aggaAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGa    <  2:952607/139‑1 (MQ=255)                                                                                                                             gaAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAg  >  1:771142/1‑139 (MQ=255)                                                                                                                             gaAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAg  <  1:1235473/139‑1 (MQ=255)                                                                                                                                aTGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccg                                                                                       >  1:298076/1‑51 (MQ=255)                                                                                                                                aTGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccg                                                                                       <  2:298076/51‑1 (MQ=255)                                                                                                                                aTGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACg                                                 >  2:591511/1‑89 (MQ=255)                                                                                                                                aTGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACg                                                 <  1:591511/89‑1 (MQ=255)                                                                                                                                    agagTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccgc                                                                                      >  1:95943/1‑48 (MQ=255)                                                                                                                                    agagTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccgc                                                                                      <  2:95943/48‑1 (MQ=255)                                                                                                                                    agagTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGt                                                                         <  2:707736/61‑1 (MQ=255)                                                                                                                                    agagTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGt                                                                         >  1:707736/1‑61 (MQ=255)                                                                                                                                        |                                                                                                                                GTGCAGTGTCGATCACGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAG  >  NC_000913/2984276‑2984538

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