Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,015,738	G→A	*193* (TAG→TAA)	mocA →	CTP:molybdopterin cytidylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,015,738	0	G	A	100.0%	230.8 / NA	71	*193* (TAG→TAA)	mocA	CTP:molybdopterin cytidylyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (37/34);  total (37/34)

GTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGA  >  NC_000913/3015607‑3015871                                                                                                                                    |                                                                                                                                      gTCACCAATATACGTCATGAACTTCTTATGGGAGATCATTGTTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTTATACTCCGTACGACTTTATAACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGt                                                                                                                                <  1:1108079/139‑1 (MQ=255)                cAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCg                                                                                                                 >  1:1348427/1‑139 (MQ=255)                 aGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGt                                                                                                                >  1:1559375/1‑139 (MQ=255)                     aCTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCaa                                                                                                            <  1:1351731/139‑1 (MQ=255)                         cTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAatat                                                                                                        >  2:1282340/1‑139 (MQ=255)                                   aTCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAAc                                                                                                                                  >  1:517214/1‑103 (MQ=255)                                   aTCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAAc                                                                                                                                  <  2:517214/103‑1 (MQ=255)                                   aTCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGt                                                                                                                                <  1:26931/105‑1 (MQ=255)                                   aTCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGt                                                                                                                                >  2:26931/1‑105 (MQ=255)                                   aTCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGAtt                                                                                                                    >  2:1559781/1‑117 (MQ=255)                                   aTCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGAtt                                                                                                                    <  1:1559781/117‑1 (MQ=255)                                   aTCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTc                                                                                                              <  1:251842/123‑1 (MQ=255)                                   aTCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTc                                                                                                              >  2:251842/1‑123 (MQ=255)                                   aTCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTc                                                                                              >  1:608594/1‑139 (MQ=255)                                   aTCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTc                                                                                              >  2:800624/1‑139 (MQ=255)                                   aTCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATCTAATCCTTc                                                                                              >  2:975660/1‑139 (MQ=255)                                        tattCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCa                                                                                         <  1:1205278/139‑1 (MQ=255)                                             cGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCAt                                                                                    <  2:1559375/139‑1 (MQ=255)                                             cGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCAt                                                                                    <  1:1282340/139‑1 (MQ=255)                                               tCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCa                                                                                         <  2:848374/132‑1 (MQ=255)                                               tCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCa                                                                                         >  1:848374/1‑132 (MQ=255)                                               tCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACt                                                                                       >  2:1408860/1‑134 (MQ=255)                                               tCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACt                                                                                       <  1:1408860/134‑1 (MQ=255)                                               tCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTa                                                                                  >  1:70321/1‑139 (MQ=255)                                               acGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCTGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTa                                                                                  >  2:1434420/2‑139 (MQ=255)                                                          aaTGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATa                                                                       <  2:608594/139‑1 (MQ=255)                                                               aGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCa                                                                                         >  1:1133096/1‑116 (MQ=255)                                                               aGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCa                                                                                         <  2:1133096/116‑1 (MQ=255)                                                                      atTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCa                                                                                                             <  1:121121/89‑1 (MQ=255)                                                                      atTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCa                                                                                                             >  2:121121/1‑89 (MQ=255)                                                                       ttttAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGAt                                                          <  2:219834/139‑1 (MQ=255)                                                                                 tGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACt                                                                                       >  1:1518363/1‑100 (MQ=255)                                                                                 tGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACt                                                                                       <  2:1518363/100‑1 (MQ=255)                                                                                    tACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGGAATTTAAAAAACGTAGTTAAACGGTTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAAGCACTGATTTTgagaaaagag                                             <  2:894739/139‑1 (MQ=255)                                                                                        ccGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCaa                                                                                                            >  2:490101/1‑72 (MQ=255)                                                                                        ccGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCaa                                                                                                            <  1:490101/72‑1 (MQ=255)                                                                                         cGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCa                                                                                         <  1:579294/90‑1 (MQ=255)                                                                                         cGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCa                                                                                         >  2:579294/1‑90 (MQ=255)                                                                                         cGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGaaaaa                                        >  1:454412/1‑139 (MQ=255)                                                                                           gacgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCa                                                                                                             <  1:557035/68‑1 (MQ=255)                                                                                           gacgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCa                                                                                                             >  2:557035/1‑68 (MQ=255)                                                                                           gacgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCtt                                                                                               >  1:1614391/1‑82 (MQ=255)                                                                                           gacgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCtt                                                                                               <  2:1614391/82‑1 (MQ=255)                                                                                           gacgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTc                                                                                              <  2:1535942/83‑1 (MQ=255)                                                                                           gacgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTc                                                                                              >  1:1535942/1‑83 (MQ=255)                                                                                           gacgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCgg                                                                                           <  2:525778/86‑1 (MQ=255)                                                                                           gacgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCgg                                                                                           >  1:525778/1‑86 (MQ=255)                                                                                           gacgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCa                                                                                         >  2:1298469/1‑88 (MQ=255)                                                                                           gacgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCa                                                                                         <  1:1298469/88‑1 (MQ=255)                                                                                           gacgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTa                                                                                  >  2:87839/1‑95 (MQ=255)                                                                                           gacgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTa                                                                                  <  1:87839/95‑1 (MQ=255)                                                                                           gacgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAAtt                                      >  2:1486930/1‑139 (MQ=255)                                                                                               acTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATa                                  <  2:1348427/139‑1 (MQ=255)                                                                                                 tttATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATaaa                                <  2:454412/139‑1 (MQ=255)                                                                                                   tattaCAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCTTTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATg                              >  1:882892/1‑139 (MQ=255)                                                                                                   tattaCAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCAt                                                                                    >  2:941414/1‑85 (MQ=255)                                                                                                   tattaCAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCAt                                                                                    <  1:941414/85‑1 (MQ=255)                                                                                                     ttaCAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATggg                            <  1:800624/139‑1 (MQ=255)                                                                                                          gCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGc                       <  1:975660/139‑1 (MQ=255)                                                                                                          gCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGc                       >  1:1298399/1‑139 (MQ=255)                                                                                                          gCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGc                       >  1:543028/1‑139 (MQ=255)                                                                                                          gCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGc                       >  2:91584/1‑139 (MQ=255)                                                                                                            gAAGGAAAGGTATAATGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCaa                     <  2:501453/139‑1 (MQ=255)                                                                                                             aaGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATa                                                                             <  1:1340302/82‑1 (MQ=255)                                                                                                             aaGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATa                                                                             >  2:1340302/1‑82 (MQ=255)                                                                                                             aaGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCaat                    >  2:628823/1‑139 (MQ=255)                                                                                                                 aaaGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCaataata                <  1:1486930/139‑1 (MQ=255)                                                                                                                        taCTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCa         >  1:1255236/1‑139 (MQ=255)                                                                                                                               aTTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCa                                                               <  1:1526984/78‑1 (MQ=255)                                                                                                                               aTTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCa                                                               >  2:1526984/1‑78 (MQ=255)                                                                                                                               aTTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGa  >  2:168425/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                      GTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGA  >  NC_000913/3015607‑3015871

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