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

TGCAGTGTCGATCACGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAG  >  NC_000913/2984277‑2984538                                                                                                                                       |                                                                                                                                tGCAGTGTCGATCACGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTAtt                                                                                                                             >  2:814680/1‑139 (MQ=255)           aTCACGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTc                                                                                                                         <  2:73397/133‑1 (MQ=255)           aTCACGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTc                                                                                                                         >  1:73397/1‑133 (MQ=255)           aTCACGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTc                                                                                                                   >  2:482542/1‑139 (MQ=255)                 aaTATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGt                                                                                                             <  1:814680/139‑1 (MQ=255)                  atatGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTc                                                                                                            <  2:232908/139‑1 (MQ=255)                  atatGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTc                                                                                                            >  1:232908/1‑139 (MQ=255)                          tGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCcac                                                                                                    >  1:713009/1‑139 (MQ=255)                                tGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCATCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGcc                                                                                                      >  2:1074184/1‑131 (MQ=255)                                tGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGcc                                                                                                      <  1:1074184/131‑1 (MQ=255)                                        tGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccgc                                                                                      >  2:522494/1‑139 (MQ=255)                                            ggAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGt                                                                                                             <  1:419446/112‑1 (MQ=255)                                            ggAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGt                                                                                                             >  2:419446/1‑112 (MQ=255)                                                        atgtatgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAATTGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGc                                                                      <  1:522494/139‑1 (MQ=255)                                                          gtatgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCATACCCTGGCCCCCGCCGATACAAAGCGTTGCCa                                                                    >  2:542807/1‑139 (MQ=255)                                                            atgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGt                                                                  <  1:482542/139‑1 (MQ=255)                                                              gtCGACTTGCCAGTCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTa                                                                                                                               >  2:1075464/1‑76 (MQ=255)                                                                  aCTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTg                                                            >  1:1389791/1‑139 (MQ=255)                                                                   cTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGcccc                                                                                         <  1:481505/109‑1 (MQ=255)                                                                   cTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGcccc                                                                                         >  2:481505/1‑109 (MQ=255)                                                                   cTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGc                                                           >  2:409190/1‑139 (MQ=255)                                                                             atcatTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACg                                                 <  1:393297/139‑1 (MQ=255)                                                                                attattCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGttt                                              >  1:651603/1‑139 (MQ=255)                                                                                       aaCAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAcc                                          >  1:1076741/1‑136 (MQ=255)                                                                                       aaCAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAcc                                          <  2:1076741/136‑1 (MQ=255)                                                                                         cAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc                                     >  2:761310/1‑139 (MQ=255)                                                                                         cAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc                                     >  2:1388313/1‑139 (MQ=255)                                                                                         cAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc                                     <  2:395843/139‑1 (MQ=255)                                                                                         cAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc                                     >  1:471206/1‑139 (MQ=255)                                                                                            aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGcc                                                                     <  2:1003052/104‑1 (MQ=255)                                                                                            aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGcc                                                                     >  1:1003052/1‑104 (MQ=255)                                                                                            aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGttt                                              <  1:143559/127‑1 (MQ=255)                                                                                            aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGttt                                              >  2:143559/1‑127 (MQ=255)                                                                                            aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                  >  1:601588/1‑139 (MQ=255)                                                                                              tAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGt                                                                                                             >  1:1188022/1‑62 (MQ=255)                                                                                              tAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGt                                                                                                             <  2:1188022/62‑1 (MQ=255)                                                                                              tAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGaa                                >  1:392615/1‑139 (MQ=255)                                                                                                 tGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGt                                                                                                             <  2:924309/59‑1 (MQ=255)                                                                                                 tGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGt                                                                                                             >  1:924309/1‑59 (MQ=255)                                                                                                      gCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCa                                                                                                           >  2:1297342/1‑56 (MQ=255)                                                                                                      gCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCa                                                                                                           <  1:1297342/56‑1 (MQ=255)                                                                                                      gCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGaa                        >  2:821793/1‑139 (MQ=255)                                                                                                      gCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCAGGCCCCCGCCGATACAAAGc                                                                           <  1:602749/88‑1 (MQ=255)                                                                                                      gCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCAGGCCCCCGCCGATACAAAGc                                                                           >  2:602749/1‑88 (MQ=255)                                                                                                        aGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAc                      <  1:1388313/139‑1 (MQ=255)                                                                                                               tATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACAGCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAtt               <  1:761310/139‑1 (MQ=255)                                                                                                                    gTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAt                                                     <  2:1413669/96‑1 (MQ=255)                                                                                                                    gTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAt                                                     >  1:1413669/1‑96 (MQ=255)                                                                                                                     tAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTg                                                            >  2:126015/1‑88 (MQ=255)                                                                                                                     tAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTg                                                            <  1:126015/88‑1 (MQ=255)                                                                                                                       aggaggAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGcc                                                                                                      <  2:466503/44‑1 (MQ=255)                                                                                                                       aggaggAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGcc                                                                                                      >  1:466503/1‑44 (MQ=255)                                                                                                                        ggaggaAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCg      >  1:113487/1‑139 (MQ=255)                                                                                                                         gaggaAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCgg     <  2:471206/139‑1 (MQ=255)                                                                                                                            gaAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAg  <  1:409190/139‑1 (MQ=255)                                                                                                                            gaAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAg  <  2:1389791/139‑1 (MQ=255)                                                                                                                               aTGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAg  >  1:1161172/1‑136 (MQ=255)                                                                                                                               aTGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAg  <  2:1161172/136‑1 (MQ=255)                                                                                                                                       |                                                                                                                                TGCAGTGTCGATCACGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAG  >  NC_000913/2984277‑2984538

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