Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,572,223	T→C	intergenic (‑178/+184)	gadB ← / ← pqqL	glutamate decarboxylase B/putative zinc peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,572,223	0	T	C	100.0%	206.2 / NA	53	intergenic (‑178/+184)	gadB/pqqL	glutamate decarboxylase B/putative zinc peptidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (27/26);  total (27/26)

GATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCAGACCTTAAAATATCACTATTAGTACTTGATTATTATTTTGAACGCATTTATAAAATTATTACATAAAAATA  >  NC_000913/1572087‑1572358                                                                                                                                         |                                                                                                                                        gATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAg                                                                                                                                      <  2:685668/140‑1 (MQ=255) gATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAg                                                                                                                                      <  1:543869/140‑1 (MQ=255)                        tGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTAt                                                                                                                             <  2:581692/126‑1 (MQ=255)                        tGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTAt                                                                                                                             >  1:581692/1‑126 (MQ=255)                        tGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTAt                                                                                                               >  2:858995/1‑140 (MQ=255)                        tGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTAt                                                                                                               >  2:395042/1‑140 (MQ=255)                           aaaGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATgcg                                                                                                            <  2:642262/140‑1 (MQ=255)                             aGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCa                                                                                                          >  2:801869/1‑140 (MQ=255)                             aGGACTCGTGGTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATgcgct                                                                                                          >  2:409207/1‑139 (MQ=255)                                  tCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCtg                                                                                                     >  1:843361/1‑140 (MQ=255)                                  tCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCtg                                                                                                     >  2:572932/1‑140 (MQ=255)                                   cGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCtgt                                                                                                    <  1:813545/140‑1 (MQ=255)                                        ttAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGtctc                                                                                               <  2:590550/140‑1 (MQ=255)                                         tAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTcc                                                                                              >  1:587266/1‑140 (MQ=255)                                               aCAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAAtt                                                                                        >  2:309927/1‑140 (MQ=255)                                                cAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTa                                                                                       >  2:558049/1‑140 (MQ=255)                                                  aaaTCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATc                                                                                     <  1:858995/140‑1 (MQ=255)                                                  aaaTCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATc                                                                                     <  1:558049/140‑1 (MQ=255)                                                  aaaTCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATc                                                                                     <  1:309927/140‑1 (MQ=255)                                                           tGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTaaaaa                                                                            >  2:524502/1‑140 (MQ=255)                                                              tatttatCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCa                                                                         <  2:587266/140‑1 (MQ=255)                                                              tatttatCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCa                                                                         <  1:500822/140‑1 (MQ=255)                                                               atttatCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCAg                                                                        >  1:763223/1‑140 (MQ=255)                                                                       tGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCAGACCt                                                                    >  1:812120/1‑136 (MQ=255)                                                                       tGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCAGACCt                                                                    <  2:812120/136‑1 (MQ=255)                                                                         agaTATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCAGACCTTAAAat                                                              >  2:287331/1‑140 (MQ=255)                                                                          gaTATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCa                                                                         >  1:510652/1‑128 (MQ=255)                                                                          gaTATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCa                                                                         <  2:510652/128‑1 (MQ=255)                                                                           atatTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCAGACCTTAAAatat                                                            <  1:395042/140‑1 (MQ=255)                                                                           atatTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCAGACCTTAAAatat                                                            <  2:743074/140‑1 (MQ=255)                                                                            tatTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCAGACCTTAAAATATc                                                           <  1:572932/140‑1 (MQ=255)                                                                                 cgcgAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCAGACCTTAAAATATCACTAt                                                      >  2:604451/1‑140 (MQ=255)                                                                                     aataatATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCAGACCTTa                                                                  >  1:855074/1‑124 (MQ=255)                                                                                     aataatATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCAGACCTTa                                                                  <  2:855074/124‑1 (MQ=255)                                                                                     aataatATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCAGACCTTAAAata                                                             >  2:186643/1‑129 (MQ=255)                                                                                     aataatATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCAGACCTTAAAata                                                             <  1:186643/129‑1 (MQ=255)                                                                                      ataataTTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCAGACCTTAAAatat                                                            >  1:768472/1‑129 (MQ=255)                                                                                      ataataTTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCAGACCTTAAAatat                                                            <  2:768472/129‑1 (MQ=255)                                                                                      ataataTTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCAGACCTTAAAATATCACTAt                                                      >  1:891950/1‑135 (MQ=255)                                                                                      ataataTTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCAGACCTTAAAATATCACTAt                                                      <  2:891950/135‑1 (MQ=255)                                                                                      ataataTTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCAGACCTTAAAATATCACTATTAGTa                                                 <  1:342104/140‑1 (MQ=255)                                                                                       taataTTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCAGACCTTAAAATATCACTATTAGTAc                                                <  2:843361/140‑1 (MQ=255)                                                                                         atatTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCAGACCTTAAAATATCACTATTAGTACtt                                              <  1:524502/140‑1 (MQ=255)                                                                                               tttCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCAGACCTTAAAATATCACTATTAGTACTTGattat                                        <  1:768519/140‑1 (MQ=255)                                                                                                 tCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCAGACCTTAAAATATCACTATTAGTACTTGattatta                                      >  2:115466/1‑140 (MQ=255)                                                                                                          aaCAATACAATATGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCAGACCTTAAAATATCACTATTAGTACTTGATTATTATTTTGAACg                             <  1:287331/140‑1 (MQ=255)                                                                                                                   atatGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCAGACCTTAAAATATCACTATTAGTACTTGATTATTATTTTGAACGCATTTATaa                    <  1:115466/140‑1 (MQ=255)                                                                                                                      tGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCAGACCTTAAAATATCACTATTAGTACTTGATTATTATTTTGAACGCATTTATAAAat                 >  2:929822/1‑140 (MQ=255)                                                                                                                      tGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCAGACCTTAAAATATCACTATTAGTACTTGATTATTATTTTGAACGCATTTATAAAat                 >  2:883154/1‑140 (MQ=255)                                                                                                                      tGAAATTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCAGACCTTAAAATATCACTATTAGTACTTGATTATTATTTTGAACGCATTTATAAAat                 >  1:754991/1‑140 (MQ=255)                                                                                                                      tGAAATTCTTGGGTGGTGGCAAGGAGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCAGACCTTAAAATATCACTATTAGTACTTGATTATTATTTTGAACGCATTTATAAAat                 >  2:913639/1‑140 (MQ=255)                                                                                                                         aaTTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCAGACCTTAAAATATCACTATTAGTACTTGATTATTATTTTGAACGCATTTATAAAattat              >  2:151276/1‑140 (MQ=255)                                                                                                                         aaTTCTTGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCAGACCTTAAAATATCACTATTAGTACTTGATTATTATTTTGAACGCATTTATAAAattat              >  1:781297/1‑140 (MQ=255)                                                                                                                              ttGGGTGGTGGCAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCAGACCTTAAAATATCACTATTAGTACTTGATTATTATTTTGAACGCATTTATAAAATTATTACAt         >  1:94396/1‑140 (MQ=255)                                                                                                                                   tggtggcAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCAGACCTTAAAATATCACTATTAGTACTTGATTATTATTTTGAACGCATTTATAAAATTATTACATaaaaa    >  1:201719/1‑140 (MQ=255)                                                                                                                                    ggtggcAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCAGACCTTAAAATATCACTATTAGTACTTGATTATTATTTTGAACGCATTTATAAAATTATTACATAAAAAt   >  1:133779/1‑140 (MQ=255)                                                                                                                                     gtggcAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCAGACCTTAAAATATCACTATTAGTACTTGATTATTATTTTGAACGCATTTATAAAATTATTACATAAAAATa  >  2:177024/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                        GATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAATCAGACCTTAAAATATCACTATTAGTACTTGATTATTATTTTGAACGCATTTATAAAATTATTACATAAAAATA  >  NC_000913/1572087‑1572358

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