Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,572,854	A→G	A783A (GCT→GCC)	pqqL ←	putative zinc peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,572,854	0	A	G	100.0%	217.0 / NA	62	A783A (GCT→GCC)	pqqL	putative zinc peptidase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (37/25);  total (37/25)

GACAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAA  >  NC_000913/1572721‑1572985                                                                                                                                      |                                                                                                                                    gACAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCg                                                                                                                               >  2:399802/1‑140 (MQ=255)   cAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAg                                                                                                                             >  2:135540/1‑140 (MQ=255)     aGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAgcg                                                                                                                           <  2:124044/140‑1 (MQ=255)            gCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATAcg                                                                                                                    >  1:524640/1‑140 (MQ=255)              cAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATAcgcg                                                                                                                  >  2:926665/1‑140 (MQ=255)                 cAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCg                                                                                                               >  2:1124913/1‑140 (MQ=255)                      gACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGa                                                                                                          >  1:866750/1‑140 (MQ=255)                              tCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAgg                                                                                                  >  1:303436/1‑140 (MQ=255)                                tttGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGc                                                                                                             >  2:48235/1‑127 (MQ=255)                                tttGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGc                                                                                                             <  1:48235/127‑1 (MQ=255)                                    gCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCg                                                                                            >  2:339492/1‑140 (MQ=255)                                    gCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCg                                                                                            <  2:28242/140‑1 (MQ=255)                                      cTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCggg                                                                                          <  2:357085/140‑1 (MQ=255)                                      cTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCggg                                                                                          <  1:1124913/140‑1 (MQ=255)                                         aGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACa                                                                                                    >  2:253266/1‑127 (MQ=255)                                         aGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACa                                                                                                    <  1:253266/127‑1 (MQ=255)                                         aGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGgtg                                                                                                >  1:427988/1‑131 (MQ=255)                                         aGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGgtg                                                                                                <  2:427988/131‑1 (MQ=255)                                         aGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAAt                                                                                       >  1:755705/1‑140 (MQ=255)                                         aGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAAt                                                                                       >  1:764674/1‑140 (MQ=255)                                              cAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATaa                                                                                  >  2:38000/1‑140 (MQ=255)                                                   gagaGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCtt                                                                             >  2:571292/1‑140 (MQ=255)                                                   gagaGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCtt                                                                             <  1:926665/140‑1 (MQ=255)                                                    agagGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTc                                                                            >  1:860693/1‑140 (MQ=255)                                                     gagGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTcc                                                                           <  2:805467/140‑1 (MQ=255)                                                      agGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTcc                                                                           >  1:243404/1‑139 (MQ=255)                                                      agGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTcc                                                                           <  2:243404/139‑1 (MQ=255)                                                        gCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCtt                                                                             >  2:946018/1‑135 (MQ=255)                                                        gCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCtt                                                                             <  1:946018/135‑1 (MQ=255)                                                         cGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACtg                                                                       >  1:747203/1‑140 (MQ=255)                                                           agaagaAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACtgtg                                                                     >  2:86839/1‑140 (MQ=255)                                                           agaagaAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACtgtg                                                                     <  1:339492/140‑1 (MQ=255)                                                           agaagaAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACtgtg                                                                     >  2:1059478/1‑140 (MQ=255)                                                           agaagaAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACtgtg                                                                     <  1:418752/140‑1 (MQ=255)                                                             aagaaACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGaa                                                                   <  1:571292/140‑1 (MQ=255)                                                             aagaaACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGaa                                                                   <  2:866750/140‑1 (MQ=255)                                                              agaaACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGaaa                                                                  <  1:806456/140‑1 (MQ=255)                                                               gaaACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAAc                                                                 <  1:873293/140‑1 (MQ=255)                                                                       gTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCca                                                         <  1:38000/140‑1 (MQ=255)                                                                            cTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGt                                                    <  1:855725/140‑1 (MQ=255)                                                                              ccAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGAAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTc                                                  >  2:347441/1‑140 (MQ=255)                                                                                aGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACA‑GTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCAtt                                               >  2:482762/1‑140 (MQ=255)                                                                                  aTGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCAttt                                              >  2:1124877/1‑140 (MQ=255)                                                                                  aTGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCAttt                                              >  2:684591/1‑140 (MQ=255)                                                                                    gCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTg                                            <  2:755705/140‑1 (MQ=255)                                                                                       tGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTc                                         >  1:970465/1‑140 (MQ=255)                                                                                                          gTAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCg                      >  1:1019168/1‑140 (MQ=255)                                                                                                           tAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTTTGCCACAGGTTCATTTTGTTCTTTTACAGTAAcc                            <  2:1018157/133‑1 (MQ=255)                                                                                                           tAGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAAcc                            >  1:1018157/1‑133 (MQ=255)                                                                                                            aGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCg                      >  1:357009/1‑138 (MQ=255)                                                                                                            aGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGtt                    >  1:1107649/1‑140 (MQ=255)                                                                                                            aGATCTTTTGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTGTACAGTAACCGATGCg                      <  2:357009/138‑1 (MQ=255)                                                                                                                  tttGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGt              >  1:345153/1‑140 (MQ=255)                                                                                                                   ttGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTc             >  2:286275/1‑140 (MQ=255)                                                                                                                    tGCCAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTcg            <  2:764674/140‑1 (MQ=255)                                                                                                                       cAGTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGCAACCGATGCGTTGTCCGTcgcgc         <  1:1124877/140‑1 (MQ=255)                                                                                                                         gTGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGa       >  2:521368/1‑140 (MQ=255)                                                                                                                          tGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAg      >  2:315885/1‑140 (MQ=255)                                                                                                                          tGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAg      >  2:431669/1‑140 (MQ=255)                                                                                                                          tGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAg      >  2:1014259/1‑140 (MQ=255)                                                                                                                          tGCGACGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAg      >  1:889222/1‑140 (MQ=255)                                                                                                                              aCGTTAAAGGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGAAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTaa  <  1:347441/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                    GACAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAA  >  NC_000913/1572721‑1572985

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