Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,572,722	A→G	C827C (TGT→TGC)	pqqL ←	putative zinc peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,572,722	0	A	G	100.0%	152.1 / NA	42	C827C (TGT→TGC)	pqqL	putative zinc peptidase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/21);  total (21/21)

TACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGACAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCG  >  NC_000913/1572586‑1572845                                                                                                                                         |                                                                                                                            tacgctGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGGCaa                                                                                                                          >  2:154746/1‑140 (MQ=255) tacgctGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGGCaa                                                                                                                          >  1:389553/1‑140 (MQ=255)         tGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGGCAAGTAAAAGc                                                                                                                  <  1:779651/140‑1 (MQ=255)                  aGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGGCAAGTAAAAGCCAGCAAATg                                                                                                         >  1:739260/1‑140 (MQ=255)                     cTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGGCAAGTAAAAGCCAGCAAATGACt                                                                                                      >  1:521054/1‑140 (MQ=255)                       gcgcTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGGCAAGTAAAAGCCAGCAAATGACTGa                                                                                                    <  2:286773/140‑1 (MQ=255)                              aCGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGGCAAGTAAAAGCCAGCAAATGACTGATATCttt                                                                                             <  2:739260/140‑1 (MQ=255)                                                     ttCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGGCAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGc                                                                      >  1:647345/1‑140 (MQ=255)                                                      tCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGGCAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGAt                                                                                  <  2:435781/127‑1 (MQ=255)                                                      tCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGGCAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGAt                                                                                  >  1:435781/1‑127 (MQ=255)                                                       cTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGGCAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCga                                                                    >  2:374426/1‑140 (MQ=255)                                                          gCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGGCAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGagaa                                                                 >  2:893059/1‑140 (MQ=255)                                                           cTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGGCAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGagaag                                                                <  1:341138/140‑1 (MQ=255)                                                           cTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGGCAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGagaag                                                                <  1:893059/140‑1 (MQ=255)                                                               cTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGGCAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAAc                                                            >  1:203665/1‑140 (MQ=255)                                                                  aTCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGGCAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCt                                                         <  2:942511/140‑1 (MQ=255)                                                                  aTCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGCTGGCAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCt                                                         >  2:951158/1‑140 (MQ=255)                                                                   tCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGGCAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTg                                                        <  2:521054/140‑1 (MQ=255)                                                                          aGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTAGGCAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTc                                                 >  2:931915/1‑140 (MQ=255)                                                                                 cGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGGCAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGc                                          <  2:203665/140‑1 (MQ=255)                                                                                      aaCCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGGCAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGtt                                     <  1:931915/140‑1 (MQ=255)                                                                                         cATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGGCAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCAc                                  <  1:692163/140‑1 (MQ=255)                                                                                            cACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGGCAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGaa                               >  1:315353/1‑140 (MQ=255)                                                                                              cTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGGCAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAata                             >  1:392926/1‑140 (MQ=255)                                                                                              cTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGGCAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAata                             >  2:281307/1‑140 (MQ=255)                                                                                              cTTCATTCGCTAACGTTAACAGTTCATCATGTCGATCTGGTTGGCAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGGATGCTCCAGATGCCTGTTCAGGAata                             >  2:923252/1‑140 (MQ=255)                                                                                                tCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGGCAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATAtt                           <  2:839035/140‑1 (MQ=255)                                                                                                 cATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGGCAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTa                          <  2:315353/140‑1 (MQ=255)                                                                                                  aTTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGGCAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTaa                         <  2:392926/140‑1 (MQ=255)                                                                                                   ttCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGGCAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAac                        >  1:1751/1‑140 (MQ=255)                                                                                                        tAACGTTAACAGTTCATCATGTCGTTCTGGTTGGCAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTa                   <  2:474883/140‑1 (MQ=255)                                                                                                        cagcGTTAACAGTTCATCATGTCGTTCTGGTTGGCAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTa                   <  2:590554/137‑1 (MQ=255)                                                                                                         aaCGTTAACAGTTCATCATGTCGTTCTGGTTGGCAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAg                  >  1:1010689/1‑140 (MQ=255)                                                                                                          aCGTTAACAGTTCATCATGTCGTTCTGGTTGGCAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGa                 <  1:281307/140‑1 (MQ=255)                                                                                                                    ttCATCATGTCGTTCTGGTTGGCAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCtt             >  2:677947/1‑134 (MQ=255)                                                                                                                    ttCATCATGTCGTTCTGGTTGGCAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCtt             <  1:677947/134‑1 (MQ=255)                                                                                                                    ttCATCATGTCGTTCTGGTTGGCAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCa       >  1:1021299/1‑140 (MQ=255)                                                                                                                      catcatGTCGTTCTGGTTGGCAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGc         >  2:740955/1‑136 (MQ=255)                                                                                                                      catcatGTCGTTCTGGTTGGCAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGc         <  1:740955/136‑1 (MQ=255)                                                                                                                         catGTCGTTCTGGTTGGCAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCg  <  1:374426/140‑1 (MQ=255)                                                                                                                              cGTTCTGGTTGGCAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCg  >  1:309130/1‑135 (MQ=255)                                                                                                                              cGTTCTGGTTGGCAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCg  <  2:309130/135‑1 (MQ=255)                                                                                                                                         |                                                                                                                            TACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGACAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCG  >  NC_000913/1572586‑1572845

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