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

ATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATA  >  NC_000913/3015605‑3015857                                                                                                                                      |                                                                                                                        aTGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAAc                                                                                                                    <  1:563018/139‑1 (MQ=255)                 tCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTaaaaaa                                                                                                                     >  2:24137/1‑122 (MQ=255)                 tCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTaaaaaa                                                                                                                     <  1:24137/122‑1 (MQ=255)                     gCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTc                                                                                                >  1:1120018/1‑139 (MQ=255)                               tGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCa                                                                                               <  1:665458/130‑1 (MQ=255)                               tGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCa                                                                                               >  2:665458/1‑130 (MQ=255)                                     aTCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTc                                                                                >  2:1388930/1‑139 (MQ=255)                                     aTCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTc                                                                                >  1:286551/1‑139 (MQ=255)                                     aTCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTc                                                                                >  1:1349271/1‑139 (MQ=255)                                          tattCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCa                                                                           >  1:1067892/1‑139 (MQ=255)                                          tattCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCa                                                                           <  2:1067892/139‑1 (MQ=255)                                              ccGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCa                                                                       <  1:1073929/139‑1 (MQ=255)                                               cGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCAt                                                                      <  1:1388930/139‑1 (MQ=255)                                                 tCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTaaaaaa                                                                                                                     >  2:637263/1‑90 (MQ=255)                                                 tCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTaaaaaa                                                                                                                     <  1:637263/90‑1 (MQ=255)                                                 tCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTa                                                                    >  1:16929/1‑139 (MQ=255)                                                 tCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACCCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTa                                                                    >  1:872319/1‑139 (MQ=255)                                                       tAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATaa                                                              <  2:1046973/139‑1 (MQ=255)                                                                                     aTACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGAtt                                                                                                      <  1:889198/69‑1 (MQ=255)                                                                                     aTACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGAtt                                                                                                      >  2:889198/1‑69 (MQ=255)                                                                                        cTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTa                                                                                                             >  2:821684/1‑59 (MQ=255)                                                                                        cTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTa                                                                                                             <  1:821684/59‑1 (MQ=255)                                                                                          ccGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCa                                                                           >  1:1170301/1‑91 (MQ=255)                                                                                          ccGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCa                                                                           <  2:1170301/91‑1 (MQ=255)                                                                                          ccGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTgagaaaagagaaaa                           >  1:351986/1‑139 (MQ=255)                                                                                          ccGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTgagaaaagagaaaa                           >  1:1132364/1‑139 (MQ=255)                                                                                             gacgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAACAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCa                                                                           <  2:1010400/88‑1 (MQ=255)                                                                                             gacgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCg                                                                              >  2:686555/1‑85 (MQ=255)                                                                                             gacgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCg                                                                              <  1:686555/85‑1 (MQ=255)                                                                                             gacgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCa                                                                           >  1:1010400/1‑88 (MQ=255)                                                                                             gacgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGAt                                                          <  1:182365/105‑1 (MQ=255)                                                                                             gacgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGAt                                                          >  2:182365/1‑105 (MQ=255)                                                                                             gacgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAAtt                        >  2:732619/1‑139 (MQ=255)                                                                                                  cTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATaa                   <  2:872319/139‑1 (MQ=255)                                                                                                   tttATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTc                                                                                <  2:1365892/77‑1 (MQ=255)                                                                                                   tttATTACAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTc                                                                                >  1:1365892/1‑77 (MQ=255)                                                                                                     tattaCAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTCGTTAAACGATTGCGTCCAAATATTTAATCCTTCCGGCACt                                                                         <  1:872935/82‑1 (MQ=255)                                                                                                     tattaCAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACt                                                                         >  2:872935/1‑82 (MQ=255)                                                                                                       ttaCAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATggg              <  2:351986/139‑1 (MQ=255)                                                                                                       ttaCAGCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATggg              <  2:1132364/139‑1 (MQ=255)                                                                                                            gCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAAt                                                                                     <  1:1133075/63‑1 (MQ=255)                                                                                                            gCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAAt                                                                                     >  2:1133075/1‑63 (MQ=255)                                                                                                            gCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCt            <  2:115918/136‑1 (MQ=255)                                                                                                            gCGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCt            >  1:115918/1‑136 (MQ=255)                                                                                                             cGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACt                                                                         <  2:246980/74‑1 (MQ=255)                                                                                                             cGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACt                                                                         >  2:967068/1‑74 (MQ=255)                                                                                                             cGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACt                                                                         >  1:246980/1‑74 (MQ=255)                                                                                                             cGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCa                                                                 <  1:1073299/82‑1 (MQ=255)                                                                                                             cGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCa                                                                 >  2:1073299/1‑82 (MQ=255)                                                                                                             cGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCa                                                 <  1:94014/98‑1 (MQ=255)                                                                                                             cGAAGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCa                                                 >  2:94014/1‑98 (MQ=255)                                                                                                             cGAAGGAAAGGTATACTGAAAATTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACt                                                                         <  1:967068/74‑1 (MQ=255)                                                                                                               aaGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATa                                                               >  2:363254/1‑82 (MQ=255)                                                                                                               aaGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATa                                                               <  1:363254/82‑1 (MQ=255)                                                                                                               aaGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATaa                   <  1:1095223/126‑1 (MQ=255)                                                                                                               aaGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATaa                   >  2:1095223/1‑126 (MQ=255)                                                                                                               aaGGAAAGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCaat      >  2:1270676/1‑139 (MQ=255)                                                                                                                   aaaGGTATACTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCaataata  <  2:286551/139‑1 (MQ=255)                                                                                                                         ataCTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCa                                                                       >  2:437701/1‑64 (MQ=255)                                                                                                                         ataCTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCa                                                                       <  1:437701/64‑1 (MQ=255)                                                                                                                            cTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATa                                                               <  1:1242854/69‑1 (MQ=255)                                                                                                                            cTGAAATTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATa                                                               >  2:1242854/1‑69 (MQ=255)                                                                                                                               aaaTTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCa                                                                           <  2:552998/54‑1 (MQ=255)                                                                                                                               aaaTTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCa                                                                           >  1:552998/1‑54 (MQ=255)                                                                                                                                 aTTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCa                                                                           >  2:243630/1‑52 (MQ=255)                                                                                                                                 aTTTAAAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCa                                                                           <  1:243630/52‑1 (MQ=255)                                                                                                                                      |                                                                                                                        ATGTCACCAATATGCGTCAGGCACTCCTTATGGGAGATCATTATTCCGTCGAAATAGAAAATGCAGAAATAATTTTAGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATA  >  NC_000913/3015605‑3015857

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