Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	220,239	C→T	L230L (CTG→CTA)	metQ ←	DL‑methionine transporter subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	220,239	0	C	T	100.0%	221.3 / NA	68	L230L (CTG→CTA)	metQ	DL‑methionine transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (33/35);  total (33/35)

CCAGCCTTTAACAGCTCCGCCGTTAAACACTTTGTTTGCTGCTTCGTAAACTTCGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGCGTCGTCCAGAGAAC  >  NC_000913/220116‑220363                                                                                                                            |                                                                                                                             ccAGCCTTTAACAGCTCCGCCGTTATACACTTTGTTTGCTGCTTCGTAAACTTCGTCAGACTGATAAGCCTAGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCCCGCGTCACGATTAGGTTTACGTACggg                                                                                                               <  2:1185905/139‑1 (MQ=255) ccAGCCTTTAACAGCTCCGCCGTTAAACACTTTGTTTGCTGCTTCGTAAACTTCGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACggg                                                                                                               <  2:933130/139‑1 (MQ=255) ccAGCCTTTAACAGCTCCGCCGTTAAACACTTTGTTTGCTGCTTCGTAAACTTCGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACggg                                                                                                               <  2:633652/139‑1 (MQ=255)         tAACAGCTCCGCCGTTAAACACTTTGTTTGCTGCTTCGTAAACTTCGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGACTCttt                                                                                                       <  2:271841/139‑1 (MQ=255)            cAGCTCCGCCGTTAAACACTTTGTTTGCTGCTTCGTAAACTTCGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATc                                                                                                    <  2:195911/139‑1 (MQ=255)               cTCCGCCGTTAAACACTTTGTTTGCTGCTTCGTAAACTTCGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGACTCttt                                                                                                       <  1:189811/133‑1 (MQ=255)               cTCCGCCGTTAAACACTTTGTTTGCTGCTTCGTAAACTTCGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGACTCttt                                                                                                       >  2:189811/1‑133 (MQ=255)               cTCCGCCGTTAAACACTTTGTTTGCTGCTTCGTAAACTTCGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTc                                                                                                 >  2:1114038/1‑139 (MQ=255)                  cgccgTTAAACACTTTGTTTGCTGCTTCGTAAACTTCGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTcaac                                                                                              >  1:1654868/1‑139 (MQ=255)                       ttAAACACTTTGTTTGCTGCTTCGTAAACTTCGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGa                                                                                         >  1:571863/1‑139 (MQ=255)                                 tgtttgCTGCTTCGTAAACTTCGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTta                                                                                                                       >  2:1465554/1‑99 (MQ=255)                                 tgtttgCTGCTTCGTAAACTTCGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTta                                                                                                                       <  1:1465554/99‑1 (MQ=255)                                 tgtttgCTGCTTCGTAAACTTCGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGAc                                                                                                            >  1:1791234/1‑110 (MQ=255)                                 tgtttgCTGCTTCGTAAACTTCGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGAc                                                                                                            <  2:1791234/110‑1 (MQ=255)                                           ttCGTAAACTTCGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTc                                                                     >  2:336307/1‑139 (MQ=255)                                            tCGTAAACTTCGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCa                                                                    >  1:1259910/1‑139 (MQ=255)                                                    ttCGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTcaaca                                                                                             >  2:433933/1‑106 (MQ=255)                                                    ttCGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTcaaca                                                                                             <  1:433933/106‑1 (MQ=255)                                                      cGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTGACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATctg                                                          >  1:1313935/1‑139 (MQ=255)                                                        tCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATctggc                                                        <  1:1114038/139‑1 (MQ=255)                                                                 tAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAgg                                               <  2:1755133/139‑1 (MQ=255)                                                                 tAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAgg                                               <  2:571863/139‑1 (MQ=255)                                                                      cTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTtacg                                                                                                                     >  2:830256/1‑64 (MQ=255)                                                                      cTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTtacg                                                                                                                     <  1:830256/64‑1 (MQ=255)                                                                        ggACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTAtt                                        >  2:997276/1‑139 (MQ=255)                                                                        ggACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTAtt                                        >  1:1521249/1‑139 (MQ=255)                                                                         gACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAgg                                               <  1:1636016/131‑1 (MQ=255)                                                                         gACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAgg                                               >  2:1636016/1‑131 (MQ=255)                                                                          aCGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTtacg                                                                                                                     >  1:227446/1‑60 (MQ=255)                                                                          aCGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTtacg                                                                                                                     <  2:227446/60‑1 (MQ=255)                                                                          aCGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATc                                                                                                    <  2:1830421/77‑1 (MQ=255)                                                                          aCGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATc                                                                                                    >  1:1830421/1‑77 (MQ=255)                                                                          aCGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGa                                      >  1:1668376/1‑139 (MQ=255)                                                                           cGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGAt                                     <  2:1259910/139‑1 (MQ=255)                                                                            gAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATa                                    <  2:1313935/139‑1 (MQ=255)                                                                              aTTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGAc                                                                                                            >  1:1610976/1‑65 (MQ=255)                                                                              aTTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGAc                                                                                                            >  1:1544481/1‑65 (MQ=255)                                                                              aTTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGAc                                                                                                            <  2:1544481/65‑1 (MQ=255)                                                                              aTTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGAc                                                                                                            <  2:1610976/65‑1 (MQ=255)                                                                               tttCTTCACGTTCTCGGCGTCTTTGTTATCTTCATGCGTCACGATTAGGTTTACGTACggg                                                                                                               <  1:593750/61‑1 (MQ=255)                                                                               tttCTTCACGTTCTCGGCGTCTTTGTTATCTTCATGCGTCACGATTAGGTTTACGTACggg                                                                                                               >  2:593750/1‑61 (MQ=255)                                                                               tttCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACgggg                                                                                                              >  2:271353/1‑62 (MQ=255)                                                                               tttCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACgggg                                                                                                              <  1:271353/62‑1 (MQ=255)                                                                                          tctcGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCtggctggc                                                    >  2:805396/1‑109 (MQ=255)                                                                                          tctcGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCtggctggc                                                    <  1:805396/109‑1 (MQ=255)                                                                                            tcGGCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTATTTCGCCGGCGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGAttt                    >  1:1562461/1‑139 (MQ=255)                                                                                              ggCGTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTgc                  <  2:349795/139‑1 (MQ=255)                                                                                                 gTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGcgcg               <  2:1654868/139‑1 (MQ=255)                                                                                                 gTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGcgcg               <  2:1521249/139‑1 (MQ=255)                                                                                                 gTCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGcgcg               <  1:997276/139‑1 (MQ=255)                                                                                                  tCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATctgg                                                         <  1:87158/96‑1 (MQ=255)                                                                                                  tCTTTGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATctgg                                                         >  2:87158/1‑96 (MQ=255)                                                                                                    tttGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAggtggt                                           <  1:801115/108‑1 (MQ=255)                                                                                                    tttGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAggtggt                                           >  2:801115/1‑108 (MQ=255)                                                                                                    tttGTTATCTTCACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGcgtcg            >  1:951560/1‑139 (MQ=255)                                                                                                              tCACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGATGGCATAGGTGGTCTTGATAACTGCCAGAGCGATTTGCGCGTAGTCCAGAGAAc  >  1:1766763/1‑139 (MQ=255)                                                                                                               cACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTcaaca                                                                                             <  2:1509180/47‑1 (MQ=255)                                                                                                               cACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTcaaca                                                                                             >  1:1509180/1‑47 (MQ=255)                                                                                                               cACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATctgg                                                         <  2:1434851/83‑1 (MQ=255)                                                                                                               cACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATctgg                                                         >  1:1434851/1‑83 (MQ=255)                                                                                                               cACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAg                                                <  2:261040/92‑1 (MQ=255)                                                                                                               cACGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAg                                                >  1:261040/1‑92 (MQ=255)                                                                                                                aCGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAgg                                                                  <  2:340101/73‑1 (MQ=255)                                                                                                                aCGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAgg                                                                  >  1:340101/1‑73 (MQ=255)                                                                                                                aCGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAAc                                  <  2:1146695/105‑1 (MQ=255)                                                                                                                aCGCGTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAAc                                  >  1:1146695/1‑105 (MQ=255)                                                                                                                  gcgTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCtggctggc                                                    >  2:783518/1‑85 (MQ=255)                                                                                                                  gcgTCACGATTAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCtggctggc                                                    <  1:783518/85‑1 (MQ=255)                                                                                                                            |                                                                                                                             CCAGCCTTTAACAGCTCCGCCGTTAAACACTTTGTTTGCTGCTTCGTAAACTTCGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGCGTCGTCCAGAGAAC  >  NC_000913/220116‑220363

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