Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,219,088	G→A	intergenic (‑12/‑406)	aer ← / → patA	fused signal transducer for aerotaxis sensory component/methyl accepting chemotaxis component/putrescine:2‑oxoglutaric acid aminotransferase, PLP‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,219,088	0	G	A	100.0%	226.1 / NA	69	intergenic (‑12/‑406)	aer/patA	fused signal transducer for aerotaxis sensory component/methyl accepting chemotaxis component/putrescine:2‑oxoglutaric acid aminotransferase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (37/32);  total (37/32)

TCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTCCTTGATTTAAAATACTTT  >  NC_000913/3218956‑3219223                                                                                                                                     |                                                                                                                                        tCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAgg                                                                                                                                   >  2:1235681/1‑139 (MQ=255)  cACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGt                                                                                                                                  >  2:1627660/1‑139 (MQ=255)      tGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTtata                                                                                                                              <  2:1019830/139‑1 (MQ=255)                         cATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTc                                                                                                           <  2:610085/139‑1 (MQ=255)                             aGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCa                                                                                                                                     >  1:516222/1‑109 (MQ=255)                             aGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCa                                                                                                                                     <  2:516222/109‑1 (MQ=255)                                 atatAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTa                                                                                                   >  1:1214575/1‑139 (MQ=255)                                           gCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTc                                                                                                           <  2:1506327/121‑1 (MQ=255)                                           gCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTc                                                                                                           >  1:1506327/1‑121 (MQ=255)                                           gCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGcccc                                                                                         >  2:1122778/1‑139 (MQ=255)                                            cAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCgg                                                                                             <  1:909820/134‑1 (MQ=255)                                            cAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCgg                                                                                             >  2:909820/1‑134 (MQ=255)                                            cAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCg                                                                                        <  1:163402/139‑1 (MQ=255)                                            cAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCg                                                                                        >  1:1649908/1‑139 (MQ=255)                                               aTCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTg                                                                                     >  1:137872/1‑139 (MQ=255)                                                 cGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACACACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTCCg                                                                                   >  1:1106025/1‑139 (MQ=255)                                                  ggTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCgg                                                                                  <  2:1649908/139‑1 (MQ=255)                                                  ggTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCgg                                                                                  <  1:1235681/139‑1 (MQ=255)                                                           cATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTg                                                                                     <  1:1016113/127‑1 (MQ=255)                                                           cATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTATTGTCAGAGTTTATGTCGGCCCCGCTg                                                                                     >  2:1016113/1‑127 (MQ=255)                                                              cAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAActt                                                                                                                 <  2:524058/96‑1 (MQ=255)                                                              cAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAActt                                                                                                                 >  1:524058/1‑96 (MQ=255)                                                              cAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCGGTTATCTTTAAcc                                                                      >  1:1346044/1‑139 (MQ=255)                                                              cAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCGGTTATCTTTAAcc                                                                      <  2:163228/139‑1 (MQ=255)                                                                gagTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCGGTTATCTTTAACCGa                                                                    <  1:699106/139‑1 (MQ=255)                                                                     gTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCGGTTATCTTTAACCGATTAAt                                                               <  2:1106025/139‑1 (MQ=255)                                                                      tATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTc                                                                                                           >  1:1324204/1‑94 (MQ=255)                                                                      tATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTc                                                                                                           <  2:1324204/94‑1 (MQ=255)                                                                           tCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCgg                                                                                  >  1:1628145/1‑114 (MQ=255)                                                                           tCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCgg                                                                                  <  2:1628145/114‑1 (MQ=255)                                                                           tCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGAtt                                                         >  1:1076593/1‑139 (MQ=255)                                                                           tCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTACGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCGGTTATCTTTaa                                                                        <  2:1613068/124‑1 (MQ=255)                                                                           tCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTACGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCGGTTATCTTTaa                                                                        >  1:1613068/1‑124 (MQ=255)                                                                               ccAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTtata                                                                                                                              <  1:1243146/66‑1 (MQ=255)                                                                               ccAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTtata                                                                                                                              >  2:1243146/1‑66 (MQ=255)                                                                               ccAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGa                                                     >  1:1292040/1‑139 (MQ=255)                                                                               ccAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGa                                                     >  2:1177401/1‑139 (MQ=255)                                                                                cAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGAt                                                    >  2:1475915/1‑139 (MQ=255)                                                                                      ggTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGccc                                                                                          <  1:907368/95‑1 (MQ=255)                                                                                      ggTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGccc                                                                                          >  2:907368/1‑95 (MQ=255)                                                                                         aTTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTg                                           <  2:1346044/139‑1 (MQ=255)                                                                                               cTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGAttt                                     <  1:1475915/139‑1 (MQ=255)                                                                                                  ggTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGt                                                                                                      >  1:1310284/1‑71 (MQ=255)                                                                                                  ggTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGt                                                                                                      <  2:1310284/71‑1 (MQ=255)                                                                                                  ggTGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTaaa                                  >  1:32007/1‑139 (MQ=255)                                                                                                    tGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCgg                                                                                  <  1:823581/89‑1 (MQ=255)                                                                                                    tGACATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCgg                                                                                  >  2:823581/1‑89 (MQ=255)                                                                                                      aCATACGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAAcaca                              <  2:137872/139‑1 (MQ=255)                                                                                                           cGGATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCt                         <  1:722482/139‑1 (MQ=255)                                                                                                             gATGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTaa                       >  2:388831/1‑139 (MQ=255)                                                                                                              aTGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCGGt                                                                                 >  2:548847/1‑80 (MQ=255)                                                                                                              aTGAGAAGACATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCGGt                                                                                 <  1:548847/80‑1 (MQ=255)                                                                                                                 agaagaCATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCagag                                                                                                       >  2:1104260/1‑55 (MQ=255)                                                                                                                 agaagaCATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCagag                                                                                                       <  1:1104260/55‑1 (MQ=255)                                                                                                                 agaagaCATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCGGt                                                                                 >  2:931833/1‑77 (MQ=255)                                                                                                                 agaagaCATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCGGt                                                                                 <  1:931833/77‑1 (MQ=255)                                                                                                                 agaagaCATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTcc                   >  1:205449/1‑139 (MQ=255)                                                                                                                 agaagaCATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTcc                   >  2:1158549/1‑139 (MQ=255)                                                                                                                 agaagaCATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTcc                   >  1:1136599/1‑139 (MQ=255)                                                                                                                   aagaCATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCg                                                                                        <  1:1319950/68‑1 (MQ=255)                                                                                                                   aagaCATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCg                                                                                        >  2:1319950/1‑68 (MQ=255)                                                                                                                   aagaCATGTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCTATTTAAACACAAATCTAATTCCtt                 >  2:932888/1‑139 (MQ=255)                                                                                                                          gTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCGGTTATCTTTAACCGAt                                                                   >  1:181225/1‑82 (MQ=255)                                                                                                                          gTTCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCGGTTATCTTTAACCGAt                                                                   <  2:181225/82‑1 (MQ=255)                                                                                                                           ttCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGAttt                                     <  1:1269291/111‑1 (MQ=255)                                                                                                                           ttCGCGTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGAttt                                     >  2:1269291/1‑111 (MQ=255)                                                                                                                             cgcgTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCGGTTATCTTTAACCGAt                                                                   >  2:461864/1‑79 (MQ=255)                                                                                                                             cgcgTCCTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCGGTTATCTTTAACCGAt                                                                   <  1:461864/79‑1 (MQ=255)                                                                                                                                  ccTATGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTCCTTGATTTAAAATACttt  <  2:205449/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                                        TCACCTGCACAAAAGTGTCATTAGCATGAGTGATATAGCTTTGCAGATCGGTAGTGGACATCAGAGTGGTATCGTCCGCCAGCGGGGTATTTTGCTGGGTGACATACGGATGAGAAGACATGTTCGCGTCCTGTGCAGGTTATATGGTTGTTAACTTCTTGTCAGAGTTTATGTCGGCCCCGCTGCGGTTATCTTTAACCGATTAATTTGATTTAGATCGCAATTTGCGATTTAAACACAAATCTAATTCCTTGATTTAAAATACTTT  >  NC_000913/3218956‑3219223

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