Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,325,794	A→G	Q350R (CAG→CGG)	yciQ →	DUF2207 domain‑containing protein YciQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,325,794	0	A	G	100.0%	196.7 / NA	55	Q350R (CAG→CGG)	yciQ	DUF2207 domain‑containing protein YciQ
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (29/26);  total (29/26)

TGTGCCGATATCGTAAATCTGATTGTGAAAGGAAAAATTCATCTGGAAGATCAGTCTGACAAGAACCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGG  >  NC_000913/1325661‑1325927                                                                                                                                      |                                                                                                                                      tgtgCCGATATCGTAAATCTGATTGTGAAAGGAAAAATTCATCTGGAAGATCAGTCTGACAAGAACCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGttatt                                                                                                                                 >  1:657616/1‑140 (MQ=255)   tgCCGATATCGTAAATCTGATTGTGAAAGGAAAAATTCATCTGGAAGATCAGTCTGACAAGAACCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTAc                                                                                                                               <  1:440236/140‑1 (MQ=255)      cGATATCGTAAATCTGATTGTGAAAGGAAAAATTCATCTGGAAGATCAGTCTGACAAGAACCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTgg                                                                                                                            >  1:180827/1‑140 (MQ=255)          atCGTAAATCTGATTGTGAAAGGAAAAATTCATCTGGAAGATCAGTCTGACAAGAACCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGc                                                                                                                        >  1:485800/1‑140 (MQ=255)           tCGTAAATCTGATTGTGAAAGGAAAAATTCATCTGGAAGATCAGTCTGACAAGAACCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCg                                                                                                                       <  1:406227/140‑1 (MQ=255)             gTAAATCTGATTGTGAAAGGAAAAATTCATCTGGAAGATCAGTCTGACAAGAACCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAg                                                                                                                         <  2:412944/136‑1 (MQ=255)             gTAAATCTGATTGTGAAAGGAAAAATTCATCTGGAAGATCAGTCTGACAAGAACCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAg                                                                                                                         >  1:412944/1‑136 (MQ=255)                     gATTGTGAAAGGAAAAATTCATCTGGAAGATCAGTCTGACAAGAACCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTa                                                                                                             >  1:1029267/1‑140 (MQ=255)                      aTTGTGAAAGGAAAAATTCATCTGGAAGATCAGTCTGACAAGAACCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTaa                                                                                                            >  1:202413/1‑140 (MQ=255)                      aTTGTGAAAGGAAAAATTCATCTGGAAGATCAGTCTGACAAGAACCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTaa                                                                                                            <  2:858990/140‑1 (MQ=255)                      aTTGTGAAAGGAAAAATTCATCTGGAAGATCAGTCTGACAAGAACCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTaa                                                                                                            >  1:88936/1‑140 (MQ=255)                       ttGTGAAAGGAAAAATTCATCTGGAAGATCAGTCTGACAAGAACCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTaaa                                                                                                           <  2:180827/140‑1 (MQ=255)                       ttGTGAAAGGAAAAATTCATCTGGAAGATCAGTCTGACAAGAACCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTaaa                                                                                                           <  2:1003743/140‑1 (MQ=255)                                gAAAAATTCATCTGGAAGATCAGTCTGACAAGAACCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAg                                                                                                  >  2:537919/1‑140 (MQ=255)                                gAAAAATTCATCTGGAAGATCAGTCTGACAAGAACCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAg                                                                                                  <  1:537919/140‑1 (MQ=255)                                    aaTTCATCTGGAAGATCAGTCTGACAAGAACCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCg                                                                                              <  2:395236/140‑1 (MQ=255)                                          tCTGGAAGATCAGTCTGACAAGAACCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTa                                                                                        >  1:178803/1‑140 (MQ=255)                                           cTGGAAGATCAGTCTGACAAGAACCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTAc                                                                                       >  2:448504/1‑140 (MQ=255)                                             ggAAGATCAGTCTGACAAGAACCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACgg                                                                                     <  1:1069163/140‑1 (MQ=255)                                             ggAAGATCAGTCTGACAAGAACCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACgg                                                                                     <  2:1059994/140‑1 (MQ=255)                                              gAAGATCAGTCTGACAAGAACCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACggg                                                                                    >  1:311391/1‑140 (MQ=255)                                               aaGATCAGTCTGACAAGAACCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACgggg                                                                                   <  2:82856/140‑1 (MQ=255)                                               aaGATCAGTCTGACAAGAACCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACgggg                                                                                   <  1:182901/140‑1 (MQ=255)                                               aaGATCAGTCTGACAAGAACCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACgggg                                                                                   <  1:258198/140‑1 (MQ=255)                                                  aTCAGTCTGACAAGAACCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGgaga                                                                                <  1:976622/140‑1 (MQ=255)                                                  aTCAGTCTGACAAGAACCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGgaga                                                                                <  2:178803/140‑1 (MQ=255)                                                  aTCAGTCTGACAAGAACCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGgaga                                                                                <  2:313670/140‑1 (MQ=255)                                                   tCAGTCTGACAAGAACCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGAc                                                                               >  2:370518/1‑140 (MQ=255)                                                   tCAGTCTGACAAGAACCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGAc                                                                               >  2:339107/1‑140 (MQ=255)                                                       tCTGACAAGAACCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCaa                                                                           >  2:5097/1‑140 (MQ=255)                                                        cTGACAAGAACCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAAc                                                                          >  2:119258/1‑140 (MQ=255)                                                         tGACAAGAACCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACa                                                                         <  2:935276/140‑1 (MQ=255)                                                         tGACAAGAACCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACa                                                                         <  2:763327/140‑1 (MQ=255)                                                         tGACAAGAACCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACa                                                                         <  2:657616/140‑1 (MQ=255)                                                         tGACAAGAACCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACa                                                                         >  1:419081/1‑140 (MQ=255)                                                                 aCCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTa                                                                 >  2:329192/1‑140 (MQ=255)                                                                 aCCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTa                                                                 <  2:235051/140‑1 (MQ=255)                                                                 aCCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTa                                                                 >  1:803668/1‑140 (MQ=255)                                                                  ccAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTAc                                                                >  1:704133/1‑140 (MQ=255)                                                                   cAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACg                                                               <  1:167526/140‑1 (MQ=255)                                                                     gCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCaa                                                                           <  1:869825/126‑1 (MQ=255)                                                                     gCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCaa                                                                           >  2:869825/1‑126 (MQ=255)                                                                                    tgtTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGa                                              >  1:479915/1‑140 (MQ=255)                                                                                        aaTGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCa                                          <  2:1029267/140‑1 (MQ=255)                                                                                        aaTGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCa                                          >  2:590028/1‑140 (MQ=255)                                                                                              ggCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAAtt                                    >  1:517916/1‑140 (MQ=255)                                                                                                              atGCGGTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTaaa                    <  2:1030844/140‑1 (MQ=255)                                                                                                                   gTATTACTGCCCGCAGAGCGGTTATTACTTGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTc               >  1:851695/1‑140 (MQ=255)                                                                                                                   gTATTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTc               >  2:142041/1‑140 (MQ=255)                                                                                                                     aTTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCtt             >  1:306170/1‑140 (MQ=255)                                                                                                                     aTTACTGCCCGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCtt             <  1:590028/140‑1 (MQ=255)                                                                                                                       tACTGCCCGCAGAGCGGTTATTACTTGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCg           <  2:851695/140‑1 (MQ=255)                                                                                                                            cccGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTa      >  1:823827/1‑140 (MQ=255)                                                                                                                              cGCAGAGCGGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTAcc    >  2:783617/1‑140 (MQ=255)                                                                                                                                cAGAGCGGTTATTACTGGAAGCATTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCgg  >  2:425156/1‑140 (MQ=255)                                                                                                                                      |                                                                                                                                      TGTGCCGATATCGTAAATCTGATTGTGAAAGGAAAAATTCATCTGGAAGATCAGTCTGACAAGAACCAGCAAATCCTGATTCGTGTTAATGAAGGCGCGACCAGAAATAATGCGGTATTACTGCCCGCAGAGCAGTTATTACTGGAAGCGTTATTTCGTAAAGGCGATAAGGTCGTTCTTACGGGGAGACGCAACAGAGTCTTACGCAGGGCATTTTTACGGATGCAGAAATTTTATCTGCCGCGTAAAAAGTCTTCGTTTTACCGG  >  NC_000913/1325661‑1325927

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