Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,696,493	T→C	G12G (GGT→GGC)	bcsE →	c‑di‑GMP‑binding protein BcsE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,696,493	0	T	C	100.0%	218.4 / NA	61	G12G (GGT→GGC)	bcsE	c‑di‑GMP‑binding protein BcsE
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (28/33);  total (28/33)

TTCTAAATTGATTCACATGCTAAATCTGATAAGTTTTAATTTCAATGGTAGGTTTATTTCTTAGCTTTCGCTAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACACGTCGC  >  NC_000913/3696356‑3696627                                                                                                                                          |                                                                                                                                       ttCTAAATTGATTCACATGCTAAATCTGATAAGTTTTAATTTCAATGGTAGGTTTATTTCTTAGCTTTCGCTAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCAt                                                                                                                                      >  1:800053/1‑140 (MQ=255) ttCTAAATTGATTCACATGCTAAATCTGATAAGTTTTAATTTCAATGGTAGGTTTATTTCTTAGCTTTCGCTAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCAt                                                                                                                                      >  1:3873/1‑140 (MQ=255)  tCTAAATTGATTCACATGCTAAATCTGATAAGTTTTAATTTCAATGGTAGGTTTATTTCTTAGCTTTCGCTAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCAtc                                                                                                                                     >  1:805005/1‑140 (MQ=255)  tCTAAATTGATTCACATGCTAAATCTGATAAGTTTTAATTTCAATGGTAGGTTTATTTCTTAGCTTTCGCTAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCAtc                                                                                                                                     >  2:670556/1‑140 (MQ=255)   cTAAATTGATTCACATGCTAAATCTGATAAGTTTTAATTTCAATGGTAGGTTTATTTCTTAGCTTTCGCTAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCAtct                                                                                                                                    <  2:504124/140‑1 (MQ=255)   cTAAATTGATTCACATGCTAAATCTGATAAGTTTTAATTTCAATGGTAGGTTTATTTCTTAGCTTTCGCTAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCAtct                                                                                                                                    <  1:360275/140‑1 (MQ=255)     aaaTTGATTCACATGCTAAATCTGATAAGTTTTAATTTCAATGGTAGGTTTATTTCTTAGCTTTCGCTAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCtca                                                                                                                                  <  1:365762/140‑1 (MQ=255)             tCAGATGCTAAATCTGATAAGTTTTAATTTCAATGGTAGGTTTATTTCTTAGCTTTCGCTAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATg                                                                                                                          >  2:606504/1‑140 (MQ=255)               acaTGCTAAATCTGATAAGTTTTAATTTCAATGGTAGGTTTATTTCTTAGCTTTCGCTAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATggg                                                                                                                        <  1:168256/140‑1 (MQ=255)                caTGCTAAATCTGATAAGTTTTAATTTCAATGGTAGGTTTATTTCTTAGCTTTCGCTAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGa                                                                                                                       >  1:341341/1‑140 (MQ=255)                 aTGCTAAATCTGATAAGTTTTAATTTCAATGGTAGGTTTATTTCTTAGCTTTCGCTAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGAt                                                                                                                      <  2:174346/140‑1 (MQ=255)                 aTGCTAAATCTGATAAGTTTTAATTTCAATGGTAGGTTTATTTCTTAGCTTTCGCTAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGAt                                                                                                                      <  2:341341/140‑1 (MQ=255)                      aaaTCTGATAAGTTTTAATTTCAATGGTAGGTTTATTTCTTAGCTTTCGCTAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCt                                                                                                                 <  1:789743/140‑1 (MQ=255)                       aaTCTGATAAGTTTTAATTTCAATGGTAGGTTTATTTCTTAGCTTTCGCTAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCt                                                                                                                 <  2:381825/139‑1 (MQ=255)                       aaTCTGATAAGTTTTAATTTCAATGGTAGGTTTATTTCTTAGCTTTCGCTAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCt                                                                                                                 >  1:381825/1‑139 (MQ=255)                       aaTCTGATAAGTTTTAATTTCAATGGTAGGTTTATTTCTTAGCTTTCGCTAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTg                                                                                                                <  1:485454/140‑1 (MQ=255)                                  ttttAATTTCAATGGTAGGTTTATTTCTTAGCTTTCGCTAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGcc                                                                                                     <  2:760276/140‑1 (MQ=255)                                    ttAATTTCAATGGTAGGTTTATTTCTTAGCTTTCGCTAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCcag                                                                                                   <  2:805005/140‑1 (MQ=255)                                            aaTGGTAGGTTTATTTCTTAGCTTTCGCTAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAggcggc                                                                                           >  1:279341/1‑140 (MQ=255)                                            aaTGGTAGGTTTATTTCTTAGCTTTCGCTAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAggcg                                                                                             <  2:279896/138‑1 (MQ=255)                                            aaTGGTAGGTTTATTTCTTAGCTTTCGCTAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAggcg                                                                                             >  1:279896/1‑138 (MQ=255)                                                  aGGTTTATTTCTTAGCTTTCGCTAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCtgg                                                                                     <  2:3873/140‑1 (MQ=255)                                                  aGGTTTATTTCTTAGCTTTCGCTAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCtgg                                                                                     <  2:800053/140‑1 (MQ=255)                                                      ttATTTCTTAGCTTTCGCTAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATTGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTggtggt                                                                                 >  2:72953/1‑140 (MQ=255)                                                      ttATTTCTTAGCTTTCGCTAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTggtggt                                                                                 >  2:602472/1‑140 (MQ=255)                                                        aTTTCTTAGCTTTCGCTAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGttt                                                                               >  2:281605/1‑140 (MQ=255)                                                              tAGCTTTCGCTAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTc                                                                         >  1:723281/1‑140 (MQ=255)                                                              tAGCTTTCGCTAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTc                                                                         >  1:15491/1‑140 (MQ=255)                                                              tAGCTTTCGCTAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTc                                                                         >  2:293661/1‑140 (MQ=255)                                                              tAGCTTTCGCTAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTc                                                                         >  1:175448/1‑140 (MQ=255)                                                               aGCTTTCGCTAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGCCATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTggtggt                                                                                 >  1:639087/1‑131 (MQ=255)                                                               aGCTTTCGCTAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTggtggt                                                                                 <  2:639087/131‑1 (MQ=255)                                                                     cGCTAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGcc                                                                  <  2:279341/140‑1 (MQ=255)                                                                      gCTAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCa                                                                 >  2:39554/1‑140 (MQ=255)                                                                        tAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATg                                                               <  2:137653/140‑1 (MQ=255)                                                                        tAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATg                                                               <  1:602472/140‑1 (MQ=255)                                                                          gTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGaa                                                             <  2:15491/140‑1 (MQ=255)                                                                          gTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGaa                                                             <  2:773032/140‑1 (MQ=255)                                                                          gTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGaa                                                             <  1:281605/140‑1 (MQ=255)                                                                          gTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGaa                                                             <  2:723281/140‑1 (MQ=255)                                                                            aaaCTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGa                                                           <  1:293661/140‑1 (MQ=255)                                                                            aaaCTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGa                                                           <  1:674525/140‑1 (MQ=255)                                                                             aaCTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGAt                                                          <  1:670556/140‑1 (MQ=255)                                                                                  aTAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATg                                                               <  2:388310/130‑1 (MQ=255)                                                                                  aTAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATg                                                               >  1:388310/1‑130 (MQ=255)                                                                                     aTCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAg                                                  >  2:906688/1‑140 (MQ=255)                                                                                            aaaaTGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCg                                           >  2:761892/1‑140 (MQ=255)                                                                                                   tGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCaaa                                    <  1:906688/140‑1 (MQ=255)                                                                                                    gACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAAc                                   >  1:252571/1‑140 (MQ=255)                                                                                                        tGAGGGACATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAAtt                               <  1:599737/140‑1 (MQ=255)                                                                                                              aCATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATcc                         <  1:39554/140‑1 (MQ=255)                                                                                                              aCATTGTGGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATcc                         <  2:699896/140‑1 (MQ=255)                                                                                                                  tgtgGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAgg                     >  2:739992/1‑140 (MQ=255)                                                                                                                  tgtgGACCCTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTAGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAgg                     >  1:147546/1‑140 (MQ=255)                                                                                                                         ccTGTATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCGTGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAAc              >  2:169754/1‑140 (MQ=255)                                                                                                                             tATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATcc                         >  1:742339/1‑125 (MQ=255)                                                                                                                             tATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATcc                         <  2:742339/125‑1 (MQ=255)                                                                                                                             tATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGca          <  2:659632/140‑1 (MQ=255)                                                                                                                             tATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGca          >  1:659632/1‑140 (MQ=255)                                                                                                                             tATTCTCTATCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGca          >  2:840298/1‑140 (MQ=255)                                                                                                                                     aTCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCGTGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACACGTcgc  <  1:169754/140‑1 (MQ=255)                                                                                                                                     aTCGGCATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGACTGAAACCGCACACGTcgc  >  1:628375/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                       TTCTAAATTGATTCACATGCTAAATCTGATAAGTTTTAATTTCAATGGTAGGTTTATTTCTTAGCTTTCGCTAGTAAACTGATAATCAGACAAAATGGTGACATGAGGGACATTGTGGACCCTGTATTCTCTATCGGTATCTCATCATTATGGGATGAGCTGCGACATATGCCAGCAGGCGGCGTCTGGTGGTTTAACGTCGATCGCCATGAAGATGCTATCAGTCTGGCGAATCAAACAATTGCATCCCAGGCTGAAACCGCACACGTCGC  >  NC_000913/3696356‑3696627

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