Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,012,376	G→T	E126* (GAA→TAA)	pqiA →	intermembrane transport protein PqiA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,012,376	0	G	T	100.0%	144.6 / NA	43	E126* (GAA→TAA)	pqiA	intermembrane transport protein PqiA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (25/18);  total (25/18)

TTACATTACTGGAAATTCCCGGCGTGCTTTTTTCTGAGGACTACGCCAGCCTCGGCACCTTTTTCCTGTTGTTTGTGCAACTGGTTCCCGCGTTTTGTCTGATAACCATTCTGTTACTGGTGAATCGCGCGGAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTGGGGAATGGCGGAGATTTTCCTCGCGGGTGTGCTGGTCAGTTTCGTTAAACTGATGGCTTACGGCAGCATTGGG  >  NC_000913/1012245‑1012507                                                                                                                                    |                                                                                                                                    ttacattacTGGAAATTCCCGGCGTGCTTTTTTCTGAGGACTACGCCAGCCTCGGCACCTTTTTCCTGTTGTTTGTGCAACTGGTTCCCGCGTTTTGTCTGATAACCATTCTGTTACTGGTGAATCGCGCGTAATTACCg                                                                                                                             >  1:842491/1‑140 (MQ=255)           ggAAATTCCCGGCGTGCTTTTTTCTGAGGACTACGCCAGCCTCGGCACCTTTTTCCTGTTGTTTGTGCAACTGGTTCCCGCGTTTTGTCTGATAACCATTCTGTTACTGGTGAATCGCGCGTAATTACCGGTCCGTTTaa                                                                                                                   >  2:620233/1‑140 (MQ=255)              aaTTCCCGGCGTGCTTTTTTCTGAGGACTACGCCAGCCTCGGCACCTTTTTCCTGTTGTTTGTGCAACTGGTTCCCGCGTTTTGTCTGATAACCATTCTGTTACTGGTGAATCGCGCGTAATTACCGGTCCGTTTAAAag                                                                                                                >  2:816603/1‑140 (MQ=255)                  cccGGCGTGCTTTTTTCTGAGGACTACGCCAGCCTCGGCACCTTTTTCCTGTTGTTTGTGCAACTGGTTCCCGCGTTTTGTCTGATAACCATTCTGTTACTGGTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCa                                                                                                            <  2:76650/140‑1 (MQ=255)                         tGCTTTTTTCTGAGGACTACGCCAGCCTCGGCACCTTTTTCCTGTTGTTTGTGCAACTGGTTCCCGCGTTTTGTCTGATAACCATTCTGTTACTGGTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCa                                                                                                     >  2:759678/1‑140 (MQ=255)                           cTTTTTTCTGAGGACTACGCCAGCCTCGGCACCTTTTTCCTGTTGTTTGTGCAACTGGTTCCCGCGTTTTGTCTGATAACCATTCTGTTACTGGTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACg                                                                                                   >  2:742825/1‑140 (MQ=255)                               tttCTGAGGACTACGCCAGCCTCGGCACCTTTTTCCTGTTGTTTGTGCAACTGGTTCCCGCGTTTTGTCTGATAACCATTCTGTTACTGGTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTg                                                                                               <  2:808278/140‑1 (MQ=255)                                            cGCCAGCCTCGGCACCTTTTTCCTGTTGTTTGTGCAACTGGTTCCCGCGTTTTGTCTGATAACCATTCTGTTACTGGTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCa                                                                                  >  1:806673/1‑140 (MQ=255)                                                      ggCACCTTTTTCCTGTTGTTTGTGCAACTGGTTCCCGCGTTTTGTCTGATAACCATTCTGTTACTGGTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTgggg                                                                        >  2:282228/1‑140 (MQ=255)                                                      ggCACCTTTTTCCTGTTGTTTGTGCAACTGGTTCCCGCGTTTTGTCTGATAACCATTCTGTTACTGGTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTgggg                                                                        >  1:653761/1‑140 (MQ=255)                                                      ggCACCTTTTTCCTGTTGTTTGTGCAACTGGTTCCCGCGTTTTGTCTGATAACCATTCTGTTACTGGTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTgggg                                                                        >  1:600110/1‑140 (MQ=255)                                                                     ttgtttgtGCAACTGGTTCCCGCGTTTTGTCTGATAACCATTCTGTTACTGGTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTGGGGAATGGCGGAGAtttt                                                         >  2:722565/1‑140 (MQ=255)                                                                     ttgtttgtGCAACTGGTTCCCGCGTTTTGTCTGATAACCATTCTGTTACTGGTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTGGGGAATGGCGGAGAtttt                                                         >  1:359169/1‑140 (MQ=255)                                                                          tgtgCAACTGGTTCCCGCGTTTTGTCTGATAACCATTCTGTTACTGGTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTGGGGAATGGCGGAGATTTTCCTcg                                                    >  1:357001/1‑140 (MQ=255)                                                                          tgtgCAACTGGTTCCCGCGTTTTGTCTGATAACCATTCTGTTACTGGTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTGGGGAATGGCGGAGATTTTCCTcg                                                    >  2:61016/1‑140 (MQ=255)                                                                            tgCAACTGGTTCCCGCGTTTTGTCTGATAACCATTCTGTTACTGGTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTGGGGAATGGCGGAGATTTTCCTcgcg                                                  >  2:541998/1‑140 (MQ=255)                                                                             gCAACTGGTTCCCGCGTTTTGTCTGATAACCATTCTGTTACTGGTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTGGGGAATGGCGGAGATTTTCCTCGCgg                                                 <  1:722565/140‑1 (MQ=255)                                                                             gCAACTGGTTCCCGCGTTTTGTCTGATAACCATTCTGTTACTGGTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTGGGGAATGGCGGAGATTTTCCTCGCgg                                                 <  1:282228/140‑1 (MQ=255)                                                                             gCAACTGGTTCCCGCGTTTTGTCTGATAACCATTCTGTTACTGGTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTGGGGAATGGCGGAGATTTTCCTCGCgg                                                 <  2:600110/140‑1 (MQ=255)                                                                                  tGGTTCCCGCGTTTTGTCTGATAACCATTCTGTTACTGGTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTGGGGAATGGCGGAGATTTTCCTCGCGGGtgtg                                            <  1:742825/140‑1 (MQ=255)                                                                                   ggTTCCCGCGTTTTGTCTGATAACCATTCTGTTACTGGTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTGGGGAATGGCGGAGATTTTCCTCGCGGGTGTGc                                           <  1:632344/140‑1 (MQ=255)                                                                                      tCCCGCGTTTTGTCTGATAACCATTCTGTTACTGGTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTGGGGAATGGCGGAGATTTTCCTCGCGGGTGTGCTgg                                        >  2:319428/1‑140 (MQ=255)                                                                                       cccGCGTTTTGTCTGATAACCATTCTGTTACTGGTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTGGGGAATGGCGGAGATTTTCCTCGCGGGTGTGCTGGt                                       <  1:541998/140‑1 (MQ=255)                                                                                        ccGCGTTTTGTCTGATAACCATTCTGTTACTGGTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTGGGGAATGGCGGAGATTTTCCTCGCGGGTGTGCTGGTc                                      >  1:575134/1‑140 (MQ=255)                                                                                           cgTTTTGTCTGATAACCATTCTGTTACTGGTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTGGGGAATGGCGGAGATTTTCCTCGCGGGTGTGCTGGTCAGt                                   >  2:519146/1‑140 (MQ=255)                                                                                            gTTTTGTCTGATAACCATTCTGTTACTGGTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTGGGGAATGGCGGAGATTTTCCTCGCGGGTGTGCTGGTCAGtt                                  <  1:519146/140‑1 (MQ=255)                                                                                             ttttGTCTGATAACCATTCTGTTACTGGTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTGGGGAATGGCGGAGATTTTCCTCGCGGGTGTGCTGGTCAGttt                                 >  1:342335/1‑140 (MQ=255)                                                                                               ttGTCTGATAACCATTCTGTTACTGGTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTGGGGAATGGCGGAGATTTTCCTCGCGGGTGTGCTGGTCAGTTTCg                               >  1:93190/1‑140 (MQ=255)                                                                                                      aTAACCATTCTGTTACTGGTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTGGGGAATGGCGGAGATTTTCCTCGCGGGTGTGCTGGTCAGTTTCGTTAAACt                        <  1:305401/140‑1 (MQ=255)                                                                                                       tAACCATTCTGTTACTGGTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTGGGGAATGGCGGAGATTTTCCTCGCGGGTGTGCTGGTCAGTTTCGTTAAACTg                       >  1:260176/1‑140 (MQ=255)                                                                                                              tCTGTTACTGGTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTGGGGAATGGCGGAGATTTTCCTCGCGGGTGTGCTGGTCAGTTTCGTTAAACTGATGGCtt                >  1:509377/1‑140 (MQ=255)                                                                                                              tCTGTTACTGGTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTGGGGAATGGCGGAGATTTTCCTCGCGGGAGTGCTGGTCAGTTTCGTTAAACTGATGGCtt                >  2:364941/1‑140 (MQ=255)                                                                                                                  ttACTGGTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTGGGGAATGGCGGAGATTTTCCTCGCGGGTGTGCTGGTCAGTTTCGTTAAACTGATGGCTTACgg            <  2:575134/140‑1 (MQ=255)                                                                                                                  ttACTGGTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTGGGGAATGGCGGAGATTTTCCTCGCGGGTGTGCTGGTCAGTTTCGTTAAACTGATGGCTTACgg            <  2:260176/140‑1 (MQ=255)                                                                                                                  ttACTGGTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTGGGGAATGGCGGAGATTTTCCTCGCGGGTGTGCTGGTCAGTTTCGTTAAACTGATGGCTTACgg            <  1:759678/140‑1 (MQ=255)                                                                                                                  ttACTGGTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTGGGGAATGGCGGAGATTTTCCTCGCGGGTGTGCTGGTCAGTTTCGTTAAACTGATGGCTTACgg            <  1:319428/140‑1 (MQ=255)                                                                                                                     cTGGTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTGGGGAATGGCGGAGATTTTCCTCGCGGGTGTGCTGGTCAGTTTCGTTAAACTGATGGCTTACGgcag         >  1:619485/1‑140 (MQ=255)                                                                                                                       ggTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTGGGGAATGGCGGAGATTTTCCTCGCGGGTGTGCTGGTCAGTTTCGTTAAACTGATGGCTTACGgcagca       <  2:359169/140‑1 (MQ=255)                                                                                                                       ggTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTGGGGAATGGCGGAGATTTTCCTCGCGGGTGTGCTGGTCAGTTTCGTTAAACTGATGGCTTACGgcagca       <  2:357001/140‑1 (MQ=255)                                                                                                                       ggTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTGGGGAATGGCGGAGATTTTCCTCGCGGGTGTGCTGGTCAGTTTCGTTAAACTGATGGCTTACGgcagca       <  1:816603/140‑1 (MQ=255)                                                                                                                       ggTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTGGGGAATGGCGGAGATTTTCCTCGCGGGTGTGCTGGTCAGTTTCGTTAAACTGATGGCTTACGgcagca       <  2:653761/140‑1 (MQ=255)                                                                                                                        gTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTGGGGAATGGCGGAGATTTTCCTCGCGGGTGTGCTGGTCAGTTTCGTTAAACTGATGGCTTACGGCAGCAt      >  1:368244/1‑140 (MQ=255)                                                                                                                        gTGAATCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTGGGGAATGGCGGAGATTTTCCTCGCGGGTGTGCTGGTCAGTTTCGTTAAACTGATGGCTTACGGCAGCAt      >  1:209193/1‑140 (MQ=255)                                                                                                                            aTCGCGCGTAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTGGGGAATGGCGGAGATTTTCCTCGCGGGTGTGCTGGTCAGTTTCGTTAAACTGATGGCTTACGGCAGCATTggg  >  1:579711/1‑140 (MQ=255)                                                                                                                                    |                                                                                                                                    TTACATTACTGGAAATTCCCGGCGTGCTTTTTTCTGAGGACTACGCCAGCCTCGGCACCTTTTTCCTGTTGTTTGTGCAACTGGTTCCCGCGTTTTGTCTGATAACCATTCTGTTACTGGTGAATCGCGCGGAATTACCGGTCCGTTTAAAAGAGCAACTGGCACGGGTGCTTTTTCAACTCAAAACCTGGGGAATGGCGGAGATTTTCCTCGCGGGTGTGCTGGTCAGTTTCGTTAAACTGATGGCTTACGGCAGCATTGGG  >  NC_000913/1012245‑1012507

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