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%	166.1 / NA	52	*193* (TAG→TAA)	mocA	CTP:molybdopterin cytidylyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/25);  total (27/25)

TGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGA  >  NC_000913/3015606‑3015871                                                                                                                                     |                                                                                                                                      tGTCACAAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACg                                                                                                                                 <  2:1028379/139‑1 (MQ=255)              cGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGAtt                                                                                                                    >  1:591765/1‑139 (MQ=255)                 cAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCg                                                                                                                 <  2:591765/139‑1 (MQ=255)                   ggCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGtt                                                                                                               >  2:728139/1‑139 (MQ=255)                    gCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTc                                                                                                              >  1:568556/1‑139 (MQ=255)                        tCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCg                                                                                                                 >  1:1052331/1‑132 (MQ=255)                        tCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCg                                                                                                                 <  2:1052331/132‑1 (MQ=255)                                    aTCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTc                                                                                              >  1:876101/1‑139 (MQ=255)                                    aTCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTc                                                                                              >  1:1218416/1‑139 (MQ=255)                                    aTCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTc                                                                                              >  1:613363/1‑139 (MQ=255)                                    aTCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTACTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTc                                                                                              >  1:972638/1‑139 (MQ=255)                                      cATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCg                                                                                            >  1:970337/1‑139 (MQ=255)                                       attattCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCgg                                                                                           >  2:1283882/1‑139 (MQ=255)                                       attattCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCgg                                                                                           <  2:568556/139‑1 (MQ=255)                                       attattCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCgg                                                                                           <  2:613363/139‑1 (MQ=255)                                         tattCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCa                                                                                         <  1:418165/139‑1 (MQ=255)                                         tattCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCa                                                                                         <  2:468092/139‑1 (MQ=255)                                         tattCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCa                                                                                         <  2:970337/139‑1 (MQ=255)                                         tattCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCa                                                                                         <  1:728139/139‑1 (MQ=255)                                         tattCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTACTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCa                                                                                         <  2:972638/139‑1 (MQ=255)                                            tCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACt                                                                                       >  1:527316/1‑138 (MQ=255)                                            tCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACt                                                                                       <  2:527316/138‑1 (MQ=255)                                              cGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCAt                                                                                    <  1:1204582/139‑1 (MQ=255)                                                tCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCAt                                                                                    <  2:402349/137‑1 (MQ=255)                                                tCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCAt                                                                                    >  1:402349/1‑137 (MQ=255)                                                       aGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTGAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAAt                                                                           <  1:810625/139‑1 (MQ=255)                                                           aaTGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCACAATGATa                                                                       <  1:7260/139‑1 (MQ=255)                                                                                          cGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGaaaaa                                        >  1:803617/1‑139 (MQ=255)                                                                                          cGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGaaaaa                                        >  1:771910/1‑139 (MQ=255)                                                                                            gacgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCa                                                                                         >  1:855939/1‑88 (MQ=255)                                                                                            gacgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCa                                                                                         <  2:855939/88‑1 (MQ=255)                                                                                              cgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTc                                                                                              >  2:621923/1‑81 (MQ=255)                                                                                              cgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTc                                                                                              <  1:621923/81‑1 (MQ=255)                                                                                                           gCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCaa                                                                                                            >  2:629613/1‑54 (MQ=255)                                                                                                           gCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCaa                                                                                                            <  1:629613/54‑1 (MQ=255)                                                                                                           gCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCa                                                                                     >  2:1176885/1‑77 (MQ=255)                                                                                                           gCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCa                                                                                     <  1:1176885/77‑1 (MQ=255)                                                                                                            cGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGAtt                                                                                                                    <  1:788263/45‑1 (MQ=255)                                                                                                            cGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGAtt                                                                                                                    >  2:788263/1‑45 (MQ=255)                                                                                                                         taCTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGa                                                                   >  2:1287291/1‑81 (MQ=255)                                                                                                                         taCTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGa                                                                   <  1:1287291/81‑1 (MQ=255)                                                                                                                           cTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCa                                                                                         <  2:363358/57‑1 (MQ=255)                                                                                                                           cTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCa                                                                                         >  1:363358/1‑57 (MQ=255)                                                                                                                              aaaTTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGt    >  2:244631/1‑139 (MQ=255)                                                                                                                                aTTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATa                                                                             >  1:1145281/1‑64 (MQ=255)                                                                                                                                aTTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATa                                                                             <  2:1145281/64‑1 (MQ=255)                                                                                                                                aTTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAAt                                                                           <  2:725433/66‑1 (MQ=255)                                                                                                                                aTTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAAt                                                                           >  1:725433/1‑66 (MQ=255)                                                                                                                                aTTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGa  >  1:948858/1‑139 (MQ=255)                                                                                                                                aTTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGa  >  1:729090/1‑139 (MQ=255)                                                                                                                                 tttAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATa                                                                       <  2:160272/69‑1 (MQ=255)                                                                                                                                 tttAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATa                                                                       >  1:160272/1‑69 (MQ=255)                                                                                                                                     |                                                                                                                                      TGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGA  >  NC_000913/3015606‑3015871

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