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

ATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAAT  >  NC_000913/3015614‑3015853                                                                                                                             |                                                                                                                    atatGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTaaa                                                                                                       >  1:837541/1‑139 (MQ=255)     gCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGAt                                                                                                   >  2:914918/1‑139 (MQ=255)         cAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCg                                                                                               >  2:115122/1‑139 (MQ=255)          aGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGt                                                                                              >  1:509751/1‑139 (MQ=255)                tCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAat                                                                                        >  2:873876/1‑139 (MQ=255)                       gggAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAAt                                                                                 <  2:649908/139‑1 (MQ=255)                       gggAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAAt                                                                                 <  1:931787/139‑1 (MQ=255)                            aTCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAAc                                                                                                                >  2:557877/1‑103 (MQ=255)                            aTCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAAc                                                                                                                <  1:557877/103‑1 (MQ=255)                            aTCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTc                                                                            >  2:158972/1‑139 (MQ=255)                                     ccGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCa                                                                   <  1:115122/139‑1 (MQ=255)                                      cGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCAt                                                                  <  2:109859/139‑1 (MQ=255)                                        tCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAGACGTAGTTAAACGATTGCGTTCAAATCTTAAATCCTTCCGGCACTCATTa                                                                >  2:722599/1‑139 (MQ=255)                                        tCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTTAAATATTTAATCCTTCCGGCACTCATTa                                                                >  2:279074/1‑139 (MQ=255)                                        tCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTa                                                                >  1:211897/1‑139 (MQ=255)                                        tCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTa                                                                >  2:881451/1‑139 (MQ=255)                                           aaaTAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCa                                                             <  1:914918/139‑1 (MQ=255)                                                     tGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCa                                                                                           >  1:388387/1‑99 (MQ=255)                                                     tGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCa                                                                                           <  2:388387/99‑1 (MQ=255)                                                           aataatTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCa                                             <  1:254905/139‑1 (MQ=255)                                                                 tttAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACATAGTTACACGATTGCGTTCCAATATTTAATCCTTCCGGCACTCATTATAATAATGAAAAAGAATCACTGAtt                                       >  2:192357/1‑139 (MQ=255)                                                                    aGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGt                                                                                                              >  1:497713/1‑65 (MQ=255)                                                                    aGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGt                                                                                                              <  2:497713/65‑1 (MQ=255)                                                                        atTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTgagaa                                >  1:531976/1‑139 (MQ=255)                                                                                  cGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTgagaaaagagaaaat                      >  1:419329/1‑138 (MQ=255)                                                                                    gacgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTaaaaaa                                                                                                                 <  1:707746/46‑1 (MQ=255)                                                                                    gacgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTaaaaaa                                                                                                                 >  2:707746/1‑46 (MQ=255)                                                                                    gacgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAAc                                                                                                                <  2:60203/47‑1 (MQ=255)                                                                                    gacgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAAc                                                                                                                >  1:60203/1‑47 (MQ=255)                                                                                    gacgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCa                                                                       <  1:524567/88‑1 (MQ=255)                                                                                    gacgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCa                                                                       >  2:524567/1‑88 (MQ=255)                                                                                    gacgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCa                                                             >  1:509197/1‑98 (MQ=255)                                                                                    gacgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCa                                                             <  2:509197/98‑1 (MQ=255)                                                                                    gacgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTgagaa                                >  1:723402/1‑127 (MQ=255)                                                                                    gacgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTgagaa                                <  2:723402/127‑1 (MQ=255)                                                                                    gacgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAAtt                    >  1:653792/1‑139 (MQ=255)                                                                                          tttATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATaaa              >  1:50316/1‑139 (MQ=255)                                                                                              ttaCAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATggg          <  2:893907/139‑1 (MQ=255)                                                                                                aCAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCt        <  2:653792/139‑1 (MQ=255)                                                                                                      aaGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACt                                                                     >  1:765322/1‑72 (MQ=255)                                                                                                      aaGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACt                                                                     <  2:765322/72‑1 (MQ=255)                                                                                                      aaGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCaat  >  2:376389/1‑139 (MQ=255)                                                                                                                   cTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCAtt                                                                 <  2:441119/63‑1 (MQ=255)                                                                                                                   cTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCAtt                                                                 >  1:441119/1‑63 (MQ=255)                                                                                                                        aTTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCt        <  2:255165/115‑1 (MQ=255)                                                                                                                        aTTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCt        >  1:255165/1‑115 (MQ=255)                                                                                                                         tttAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCt        >  1:237903/1‑114 (MQ=255)                                                                                                                         tttAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCt        <  2:237903/114‑1 (MQ=255)                                                                                                                             |                                                                                                                    ATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAAT  >  NC_000913/3015614‑3015853

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: ‑