Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,088,257	T→C	intergenic (+397/‑27)	metK → / → galP	S‑adenosylmethionine synthetase/D‑galactose transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,088,257	0	T	C	100.0%	242.8 / NA	67	intergenic (+397/‑27)	metK/galP	S‑adenosylmethionine synthetase/D‑galactose transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (37/30);  total (37/30)

ACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTG  >  NC_000913/3088140‑3088392                                                                                                                      |                                                                                                                                        aCTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCa                                                                                                                             >  1:86143/1‑130 (MQ=255) aCTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCa                                                                                                                             <  2:86143/130‑1 (MQ=255)      attattTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGcatc                                                                                                               >  1:504084/1‑139 (MQ=255)       ttattTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGcatca                                                                                                              >  1:1016233/1‑139 (MQ=255)               aTGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGAc                                                                                                      >  2:1252761/1‑139 (MQ=255)                   aCTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTa                                                                                                  <  2:1016233/139‑1 (MQ=255)                   aCTTAAATCATAACTAAGATAAATGTTAGTGTAACCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTa                                                                                                  <  1:967563/139‑1 (MQ=255)                              aaCTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAgggg                                                                                       >  1:1166420/1‑139 (MQ=255)                              aaCTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAgggg                                                                                       <  2:504084/139‑1 (MQ=255)                               aCTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGc                                                                                      >  2:396026/1‑139 (MQ=255)                               aCTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGc                                                                                      >  2:1323092/1‑139 (MQ=255)                               aCTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGc                                                                                      >  2:68562/1‑139 (MQ=255)                                   aGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAAt                                                                                                                                  >  2:611440/1‑91 (MQ=255)                                   aGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAAt                                                                                                                                  <  1:611440/91‑1 (MQ=255)                                   aGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAgg                                                                                                                    >  2:466409/1‑105 (MQ=255)                                   aGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAgg                                                                                                                    <  1:466409/105‑1 (MQ=255)                                                    aGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAg                                                                                                                     >  2:1049035/1‑87 (MQ=255)                                                    aGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAg                                                                                                                     <  1:1049035/87‑1 (MQ=255)                                                          tACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAg                                                                                                                     <  1:1495563/81‑1 (MQ=255)                                                          tACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAg                                                                                                                     >  2:1495563/1‑81 (MQ=255)                                                          tACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAgggg                                                                                       >  2:1046627/1‑111 (MQ=255)                                                          tACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAgggg                                                                                       <  1:1046627/111‑1 (MQ=255)                                                          tACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAgg                                                                                         >  1:1132851/1‑109 (MQ=255)                                                          tACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAgg                                                                                         <  2:1132851/109‑1 (MQ=255)                                                            cacTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTc                                                         <  1:1323092/139‑1 (MQ=255)                                                                 aTGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTa                                                                                                  <  2:835452/93‑1 (MQ=255)                                                                 aTGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTa                                                                                                  >  1:835452/1‑93 (MQ=255)                                                                 aTGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAggg                                                                                        >  1:420359/1‑103 (MQ=255)                                                                 aTGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAggg                                                                                        <  2:420359/103‑1 (MQ=255)                                                                 aTGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCtt                                                    >  1:870245/1‑139 (MQ=255)                                                                  tgtgATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTc                                                   >  2:547565/1‑139 (MQ=255)                                                                   gtgATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTcc                                                  >  1:858021/1‑139 (MQ=255)                                                                      aTTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTg                                               >  1:972573/1‑139 (MQ=255)                                                                      aTTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTg                                               >  1:149243/1‑139 (MQ=255)                                                                       tttGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATaa                                                                                                                                >  1:718966/1‑57 (MQ=255)                                                                       tttGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATaa                                                                                                                                <  2:718966/57‑1 (MQ=255)                                                                           cTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACACAAAATAACCATATTGGAGCGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGct                                          >  2:786558/1‑139 (MQ=255)                                                                           cTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCa                                                                        <  1:644730/109‑1 (MQ=255)                                                                           cTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCa                                                                        >  2:644730/1‑109 (MQ=255)                                                                           cTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGct                                          <  2:1166420/139‑1 (MQ=255)                                                                                caTCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCt                                                 <  2:1467100/127‑1 (MQ=255)                                                                                caTCTTTTTACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCt                                                 >  1:1467100/1‑127 (MQ=255)                                                                                     ttttACGTCGTACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGAt                                >  1:124125/1‑139 (MQ=255)                                                                                           gtcgtACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTActct                          <  2:124125/139‑1 (MQ=255)                                                                                             cgtACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGcat                                                                                                                >  1:1412349/1‑51 (MQ=255)                                                                                             cgtACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGcat                                                                                                                <  2:1412349/51‑1 (MQ=255)                                                                                             cgtACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGa                                                                                                       >  2:1276711/1‑60 (MQ=255)                                                                                             cgtACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGa                                                                                                       <  1:1276711/60‑1 (MQ=255)                                                                                             cgtACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCttt                        >  2:1441103/1‑139 (MQ=255)                                                                                             cgtACTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCttt                        >  1:328463/1‑139 (MQ=255)                                                                                                aCTCACCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGc                     <  1:1441103/139‑1 (MQ=255)                                                                                                    aCCTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTgg                 >  1:281418/1‑139 (MQ=255)                                                                                                      cTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAAt                                                                      <  2:1446654/84‑1 (MQ=255)                                                                                                      cTATCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAAt                                                                      >  1:1446654/1‑84 (MQ=255)                                                                                                        aTCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGatat             <  2:281418/139‑1 (MQ=255)                                                                                                        aTCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGatat             <  2:681621/139‑1 (MQ=255)                                                                                                        aTCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGCCTGCTTCCTTGCCGCACTGGCGGGATTACTCTTTGGCCTGGatat             <  2:870245/139‑1 (MQ=255)                                                                                                         tCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCg                                            >  1:1157493/1‑107 (MQ=255)                                                                                                         tCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCg                                            <  2:1157493/107‑1 (MQ=255)                                                                                                         tCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATc            <  1:786558/139‑1 (MQ=255)                                                                                                         tCTTAATTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATc            >  1:564722/1‑139 (MQ=255)                                                                                                              aTTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAgggg                                                                                       <  2:543986/59‑1 (MQ=255)                                                                                                              aTTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAgggg                                                                                       >  1:543986/1‑59 (MQ=255)                                                                                                              aTTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAggg                                                                                        >  2:1442610/1‑58 (MQ=255)                                                                                                              aTTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAggg                                                                                        <  1:1442610/58‑1 (MQ=255)                                                                                                              aTTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGAt                                <  1:1254630/114‑1 (MQ=255)                                                                                                              aTTCACAACAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGAt                                >  2:1254630/1‑114 (MQ=255)                                                                                                                   caacaaAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTgg                                      <  2:841891/103‑1 (MQ=255)                                                                                                                   caacaaAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTg  >  2:1194001/1‑139 (MQ=255)                                                                                                                      |                                                                                                                                        ACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTG  >  NC_000913/3088140‑3088392

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